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ABSTRACT 

This  report  is  a user’s  guide  for  a computerized  mathematical 
model  development  in  support  of  the  Aircraft  Undersea  Sound  Experiment 
(AUSEX)  Program.  The  AUSEX  Program  is  funded  and  directed  by  the 
Defense  Advanced  Research  Projects  Agency  (DARPA). 

The  subject  mathematical  model  is  termed  the  AUSEX  Acoustic 
Propagation  Model  and  it  describes  the  acoustic  propagation  from 
a moving  sound  source  in  air,  to  and  across  a rough  air-water  inter- 
face, and  subsequently  through  the  water  to  an  arbitrarily  located 
point-acoustic  receiver. 

This  report  describes  the  computer  program  architecture  and 
the  input  and  output  data  associated  with  the  program’s  use.  Two 
examples  of  the  program’s  application  are  included.  In  addition, 
the  complete  program  listing  for  the  Portron  IV  coding  is  included. 

A description  of  the  mathematical  basis  for  the  Acoustic 
Propagation  Model  Is  contained  in  a companion  report,  BBN  TM  W311, 
’’AUSEX  Air-Water  Acoustic  Program  Model’’. 
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SECTION  1 
INTRODUCTION 

1.1  Overview  and  Objective  of  the  Report 

The  Defense  Advanced  Research  Projects  Agency  (DARPA)  has 
been  sponsoring  a research  and  development  program  termed  AUSEX 
(Aircraft  Undersea  Sound  Experiments).  Briefly,  the  objective 
of  th2  AUSEX  Program  is  the  development  of  generic  algorithms 
for  the  detection,  classification,  and  tracking  of  air  vehicles 
by  their  underwater  acoustic  signatures. 

The  principal  elements  and  their  functional  relationship 
are  indicated  in  Fig.  1-1.  One  of  the  primary  inputs  t:>  the 
detection,  classification,  tracking  algorithm  design  is  the 
outputs  from  the  air-water  acoustic  propagation  model.1*  This 
model  was  designed  to  provide  the  algorithm  designer  with  i 
detailed  description  of  the  signal  field  at  a point  in  the 
ocean  as  a function  of  time,  the  signal  being  the  radiated 
acoustic  energy  arsociated  with  an  air  vehicle  (i.e.  fixed  and 
rotary  wing  aircraft  and  cruise  missiles). 

Other  AUSEX  Program  elements  have  been  concerned  with  the 
development  of  detailed  acoustic  source  characterizations  of 

p 

air  vehicles  of  interest.  And  another  dealt  with  the  quanti- 
fication of  the  effects  of  real  ocean  surfaces  on  the  propaga- 
tion of  sound  from  air  to  water. ^ 

The  air-water  propagation  model  was  computer  programmed  as 
a tool  for  the  detection  algorithm  designers.  This  technical 
memorandum  is  intended  to  serve  as  a user's  guide  and  as  such 


•Superscript  numbers  identify  references  as  listed  on  page  R-l. 
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FIGURE  1-1  PRINCIPAL  ELEMENTS  OF  THE  AUSEX  PROGRAM 
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it  describes  the  computer  program,  its  required  input,  and  sub- 
sequent output.  Typical  examples  of  its  utilization  are  given. 
The  program  listing  in  Fortran  IV  for  the  BBN  TENEX  System  is 
also  included.  A description  of  the  mathematical  basis  for  the 
Propagation  Model  is  given  in  BBN  Technical  Memorandum  TM  W311, 
"AUSEX  Air-Water  Program  Model".  This  report  is  in  publication. 

1.2  Report  Organization 

Section  2 is  a summary  of  the  overall  air-water  acoustic 
propagation  program.  The  program  is  designed  as  an  executive 
with  several  subroutines  which  concern  such  things  as  the  air- 
water  interface  model,  the  several  distinct  underwater  acoustic 
propagation  modes,  the  atmospheric  acoustic  propagation,  and  the 
dynamics  of  the  air-source  underwater-receiver  encounter. 

Section  3 displays  some  typical  results  in  manually  plotted 
form.  In  addition,  this  section  includes  some  discussion  of  the 
individual  results  for  the  example  cases  chosen. 
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SECTION  2 
PROGRAM  SUMMARY 

2.1  Intended  Utilization 

The  ultimate  objective  of  the  AUSEX  Program  is  the  develop- 
ment of  generic  algorithms  for  the  detection,  classification, 
and  tracking  of  air  vehicles  by  their  underwater  acoustic  signa- 
tures. The  AUSEX  air-water  acoustic  propagation  model  is 
intended  to  provide  the  detection  algorithm  designer  with  a 
detailed,  time  dependent  description  of  the  acoustic  field 
produced  at  a (moving)  point  in  the  ocean  by  an  air  vehicle. 

The  model  output  provides  a primary  information  base  which  the 
designer  will  utilize  to  construct  and  optimize  detection/classi- 
fication/tracking algorithm  schemes. 

The  user  specifies  atmospheric  and  ocean  environmental 
parameters  and  air  vehicle  and  receiver  tra'’k  parameters.  The 
model  code  then  marches  through  time  as  the  encounter  unfolds, 
calculating  time  histories  of  the  following  quantities: 

• Transmission  loss  from  air  vehicle  to  the  receiver 
for  each  underwater  mode  of  propagation  (i.e.,  direct 
path,  bottom  bounce,  surface  duct  and  convergence 
zone) . 

• Depression/elevation  and  azimuthal  arrival  angles 
for  each  propagation  mode. 

• Received  frequency  by  propagation  mode . 

• Range  variables. 

In  obtaining  these  variables,  the  model  accounts  for  atmospheric 
propagation,  air-water  interface  transmission,  underwater  propa- 
gation and  source/receiver  geometry  dynamics . 
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The  parameter  space  c.  vered  by  the  model  Is  given  in 
Table  2-1 


The  family  of  model  outputs  for  anticipated  scenarios 
provides  a hypothesis  space  for  the  algorithm  designer  to  use 
in  the  synthesis  of  detection  schemes. 

2.2  Program  Architecture 

Figure  2-1  shows  the  basic  program  structure.  There  are 
two  independent  programs,  one  of  which  calculates  direct  path 
and  bottom  bounce  outputs,  while  the  other  computes  convergence 
zone  and  surface  duct  quantities. 

For  each  of  the  programs,  the  input  data  is  inserted  and 
the  program  begins  at  a user-specified  initial  time  and  calls 
upon  each  of  a pair  of  propagation  subroutines,  in  turn.  The 
subroutines  are  almost  completely  self-contained  and  mutually 
independent.  (The  subroutines  are  preceded  by  a sub-program 
which  calculates  a small  number  of  common  quantities.)  Each 
of  the  propagation  subroutines  calculates  all  the  pertinent 
output  for  only  one  type  of  underwater  propagation  mode. 
Incorporated  in  the  subroutines  are  atmospheric  propagation, 
air-water  interface  transmission,  underwater  propagation  and 
source-receiver  geometry  dynamics.  For  a given  time  instant, 
each  of  the  subroutines  outputs  time  variables  (emission  time, 
arrival  time) , transmission  loss,  depression/elevation  and 
azimuthal  arrival  angles,  received  frequency  and  range  variables 
(range  at  emission  time,  range  at  arrival  time).  The  time  is 
then  increased  by  a user-specified  increment  and  the  calculation 
process  repeated  until  the  incremented  time  variable  exceeds  a 
user-specified  final  time. 


-5- 
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TABLE  2-1 

PROGRAM  PARAMETER  SPACE 


Air  source  parameters 
Character : 

Speed: 

Heights : 

Lateral  offsets: 
Radiation  frequency: 
Track: 

Receiver  parameters 
Character: 

Speed : 

Depths : 

Track: 

Range  Interval 

Closest  approach  to  Inner 
zone . 

Sea  3tate 

Character: 


point,  omnidirectional 

any  speed  not  large  compared  to 
the  sound  speed. in  air 

any 

> 0 yards 

> 0 - 10,000  Hz 

straight  line  motion  at  constant 
height  in  any  direction 


point,  omnidirectional 

small  compared  with  water  sound 
speed 

any  small  fraction  of  ocean 
depth  (for  instance,  < 600  ft) 

straight  line  motion  at 
constant  depth 


edge  of  second  convergence 


fully  developed  or  smooth 
(partially  arisen  seas  and 
swells  not  included) 

any 


Wind  speed: 

Wind  direction: 


any 
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Table  2-1  continued. 


Atmosphere 

Air  temperature: 

Ocean 

Bottom: 

Bottom  types: 

Sound  speed  profiles : 


-10°F  to  90°F 


flat 

LRAPP  types  3 and  5^ 

any  two  or  three  linearly  seg 
mented  profile  with  a depth 
excess  approximating  a deep 
ocean  winter  profile  (three 
segments)  or  summer  profile 
(two  segments) 


Depth:  deep  ocean  depths 

shallow  depths  out  to  horizontal 
ranges  where  triple  bottom  bounce 
dominates 
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INPUT 


1.  ENVIRONMENTAL  PARAMETERS 

2.  SOURCE  AND  RECEIVER  TRACK 

PARAMETERS 

.3.  SOURCE  FREQUENCY 
4.  TIME  STEP  SIZE 


t E t 


initial 


PROPAGATION  SUBROUTINES 
DIRECT  PATH  AND  BOTTOM  BOUNCE 


CONVERGENCE  ZONE  AND 
SURFACE  DUCT 


ERlilT 

1.  TIME  VARIABLES 

2.  TRANSMISSION  LOSS  BY  MODE 

3.  D/E  AND  AZIMUTHAL  ARRIVAL 

ANGLE  BY  MODE 

4.  RECEIVED  FREQUENCY  BY  MODE 
'5.  RANGE 


t ■ t ♦ at 


l> ‘final S-HL 


FIGURE  2-1  Schematic  of  Propagation  Model  Program 
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Initial  and  final  time  and  time  increments  are  selected 
by  the  user  on  the  basis  of  the  range  interval  to  be  covered 
and  the  temporal  detail  required. 

No  computer  plotting  capabilities  are  included  in  the 
current  model  code;  all  results  must  be  hand  plotted. 

The  program  is  written  in  FORTRAN  IV.  It  is  operational 
on  the  Interactive  Sciences  Corporation  system  and  the  BBN 
T.5NEX  system. 

Detailed  descriptions  of  the  component  programs  are  given 
in  Appendices  A and  B. 

2.3  Input  Data 

The  program  accepts  and  inputs  the  quantities  listed  in 
Table  2-2.  Refer  to  Figs.  2-2  and  2-3  for  a pictorial  defini- 
tion of  the  variables. 


The  air  source/receiver/time  inputs  define  the  encounter 
scenario  and  radiation  frequency  of  the  source.  The  atmospheric 
inputs  quantify  atmospheric  propagation  effects  and  the  character 
of  the  air-water  interface.  Sound  velocity  profiles  and  bottom 
reflection  losses  are  specified  by  the  ocean  medium  input  group. 

Appendices  A and  B contain  detailed  descriptions  of  the 
input  data  format. 

2.4  Output  Data 

Table  2-3  lists  the  major  program  outputs.  See  Fig.  2-2 
for  a schematic  definition  of  some  of  the  output  variables. 

For  each  value  of  absolute  time,  the  output  is  grouped  by 
propagation  mode. 
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TABLE  2-2 

PROPAGATION  MODEL  PROGRAM  INPUTS 


Times 


starting  time  of  encounter  (in  secs) 

final  time  (in  secs) 

step  size  in  time  (in  secs) 

Air  Source 


v , v 
ax*  ay 


cia*  yia 


h 

f 

Receiver 

D 

vs 

x 

•is 


cartesian  components,  in  kts,  of  air 
source  velocity 

vector  position  of  air  source,  in  ft,  at 
time  t = 0 ("initial  coordinates") 

height,  in  ft,  of  air  source 

frequency  of  air  source,  in  Hz 


depth  (ft) 

receiver  speec'  along  x-axis  (in  kts) 

x-position  of  receiver  (in  ft)  at  time 
t * 0 ("initial  coordinates") 


Atmosphere 

air  temp. — either  gre?"r  or  less  than  50°F 
— wind  speec*  (kts 


sx*nsy 


Ocean  Medium 


cmin*^min 


cartesian  direction  cosines  of  wind 
velocity  vector  (dimensionless) 


— surface  sound  speed  (ft/sec) 

sound  speed  at  sound  speed  profile 
minimum  (ft/sec),  depth  of  profile 
minimum  (ft) 


-10- 
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Table  2-2  continued. 


cb»Dg  — sound  speed  at  ocean  bottom  (ft/sec), 

depth  to  bottom  (ft) 

cmax,(*s  — sound  speed  at  bottom  of  surface  duct  (if 

one  exists)  (ft/sec),  depth  of  surface 
duct  (ft) 

Bottom  type  — either  LRAPP  type  3 or  5 


-11- 


Bolt  Beranek  and  Newman  Inc. 


TM  W307 
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FIGURE  2-2  Schematic  of  Surface  and  Source  / Receiver  Variables. 
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TABLE  2-3 
PROGRAM  OUTPUTS 


Propagation  Mode  Characterization 

- Propagation  mode  type  (i.e.,  single  bottom  bounce, 
direct  path,  etc.) 


Time  Variables 

- absolute  time 

- arrival  time  when  ray  reaches  receiver 

- time  difference  between  arrival  time  and  time  of 
arrival  of  direct  ray  from  CPA  point 


Range  Variables 

- horizontal  source/receiver  range  at  time  ray  is 
emitted  by  air  source  (R’(t)  in  Fig.  2) 

- horizontal  -source/receiver  range  at  time  ray  is 
received 


Arrival  Angle  Variables 

- depression/elevation  angle (s)  of  arriving  ray(s) 

at  arrival  time  of  ray(s)  for  each  propagation  mode 
("D/E"  in  Fig.  2) 

- azimuthal  arrival  angle  of  arriving  ray(s)  at 
arrival  time  of  ray(s)  for  each  propagation  mode 
(♦  in  Fig.  2) 


Received  Frequency 

- ratio  of  received  frequency  to  source  frequency  of 
arrival (s)  at  arrival  time  for  each  propagation 
mode 
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Table  2-3  continued. 


Transmi. oion  Loss  Variables 


transmission  loss  of  received  pressure  level  for  each 
serrate  propagation  mode,  at  appropriate  arrival  time, 
referred  to  air  source  level  at  1 yd  from  the  source, 
including  water  volumetric  and  bottom  losses  (where 
appropriate)  and  air  volumetric  losses,  for  both  smooth 
and  wind-driven  ocean  surfaces. 
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Following  the  mode  identification,  a number  of  time 
variables  are  given.  The  first  is  the  time  when  the  source 
radiates.  Since  the  source  is  assumed  to  be  continuously 
radiating,  this  time  is  the  same  as  the  absolute  time.  The 
arrival  time  when  the  ray  reaches  the  receiver  is  given  next. 
Finally,  the  time  difference  between  the  arrival  time  and  an 
arbitrarily  defined  reference  time  is  given.  This  reference 
time,  which  only  has  physical  significance  if  the  distance  of 
closest  approach  is  less  than  several  kiloyards,  is  the  time 
of  arrival  of  the  direct  path  ray.  The  program  calculates 
direct  path  information  as  though  the  ocean  were  isospeed, 
regardless  of  the  actual  sound  velocity  profile.  Physically, 
this  is  a good  approximation  for  source/receiver  nnges  of 
several  kiloyards  or  less.  For  greater  ranges,  a direct  path 
connecting  the  effective  surface  source  and  the  receiver  will 
not  exist  because  of  refraction.  The  program,  however,  will 
still  calculate  a direct  path  arrival  tine.  The  user  Is 
consequently  advised  to  ignore  the  time  difference  output 
variable  when  the  encounter  CPA  is  more  than  one  or  two  kilo- 
yards . 

Two  range  variables  are  calculated.  The  first  is  the 
horizontal  source  /receiver  range  at  the  time  the  source  radiates 
(same  as  absolute  time).  The  second  is  the  horizontal  source/ 
receiver  range  when  the  signal  is  received  (i.e.,  at  the  arrival 
time) . 

Angle  information  includes  the  depression/elevation  and 
azimuthal  arrival  angle (s)  at  the  appropriate  arrival  time. 

The  ratio  of  received  frequency  to  source  frequency 
(Doppler  shift)  for  each  arrival  is  given. 
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The  transmission  loss  for  each  arrival  is  calculated  for 
both  a smooth  and  wind-driven  air/water  interface,  at  a single 
source  frequency.  The  transmission  loss  is  defined  as  10  log 
of  the  ratio  of  the  mean  squared  received  pressure  to  the 
square  of  the  air  source  pressure  at  one  yard  from  the  source. 

Broadband  noise  transmission  loss  may  be  treated  by  sub- 
dividing the  noise  spectrum  into  a number  of  frequency  bands 
and  using  the  single  frequency  results  as  an  estimate  of  the 
transmission  loss  for  each  sub-band. 
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SECTION  3 
TYPICAL  RESULTS 


Sample  < ltputs  of  the  programmed  algorithm* are  presented 
and  discussed  in  this  section. 

3.1  Representative  North  Atlantic  Environmental  Conditions 

The  input  parameters  selected  for  the  first  example  are 
summarized  in  Table  3-1. 

The  sound  speed  profiles  are  chose-'  to  be  representative 
of  North  Atlantic  mid-latitude  profiles,  and  represent  a com- 
promise between  the  FNWC  profiles  of  area  50^  and  the  profiles 
for  Marsden  square  078,  as  given  in  the  LRAPP  volumes.5*  The 

bottom  type  is  a low  loss  type  which  occurs  frequently  in  the 
North  Atlantic  area. 

Sea  State  3 (average  wind  speed  8.8  kts ) occurs  with  W 
probability  in  the  summer  and  51%  probability  in  the  winter  in 
Marsden  square  078  and  is  taken  to  be  representative  of  the 
area’s  general  wind  conditions.  There  is  no  reported  pre- 
dominant wind  direction;  a direction  of  *15°  with  respect  to 
the  air  source  track  is  assumed  arbitrarily  and  should  provide 
sufficient  generality. 


Air  temperatures  are  taken  to  be  less  than  50°F  for  the 
winter  profile  and  greater  than  50°F  for  the  summer  profile. 
(The  program  actually  uses  air  loss  curves  representative  of 
0°F  for  the  case  when  the  temperature  is  less  than  50°F;  when 
the  air  temperature  is  greater  than  50°F,  it  uses  air  loss 
curves  representative  of  70°F.) 


* Fleet  Numerical  Weather  Central  (FNWC)  and  Long  Range  Acoustic 
Propagation  Project  (LRAPP)  have  surveyed  the  deep  ocean 
archival  data  and  established  representative  sound  speed 
profiles  and  bottom  loss  curves . 
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TABLE  3-1 

EXAMPLE  INPUT  PARAMETERS 


SOURCE/RECEIVER 

Air  Source: 

Speed  ■ 220  kts 

Height  « 10,000  ft 

Frequency  = 150  Hz 


Receiver; 

Speed  = 7 kts 

Depth  * 400  ft 

Source/Receiver  Tracks : 

Track  directions:  parallel 

Closest  point  of  approach:  4 kyds 

Range:  extends  from  CPA  to  beyond  first  convergence 

zone  120  kyd) 

Time  of  CPA:  t = 0 


OCEAN 


Sound  Speed  Profiles: 


Winter Summer 


Depth 

Sound  Speed 

Depth 

Sound  Speed 

(ft) 

(ft/sec ) 

(ft) 

(ft/sec ) 

0 

4990.4 

0 

5052 

420 

4998.5 

3440 

4875 

3720 

4886.5 

15660 

5053 

15540 

5050 

Bottom  Type:  LRAPP  bottom  type  3 
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Table  3-1  continued. 


ATMOSPHERE 

Wind; 

Wind  speed:  8.8  kts 

Direction:  45°  with  respect  to  source  (receiver)  track 

Air  Temperature: 

Winter:  Less  than  50°P 

Summer:  Greater  than  50°P 
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3.1.1  Summer 

The  program  output  for  the  summer  profile  is  shown  in 
Pigs.  3-1  and  3-2  (CPA = 4 kyds).  Rough  surface  effects  are, 
in  general,  insignificant  for  the  parameters  of* this  example 
and  will  not  be  discussed.  (Section  3«2  presents  an  example 

for  which  effects  are  substantial  and  are  elaborated  upon 

in  detail. ) 

All  results  are  plotted  against  time  on  the  abscissa. 

The  time  shown  in  the  figure  is  the  difference,  in  secs, 
between  the  arrival  time  of  the  signal  and  the  arrival  time 
of  the  direct  ray  from  the  source  at  the  CPA  point.  This 
latter  arrival  time  is  merely  an  arbitrarily  selected  time 
origin.  It  is  generally  very  close  to  the  time  of  CPA,  if 
the  CPA  distance  is  not  exceedingly  large.  With  the  source/ 
receiver  track  configuration  assumed,  negative  time  differences 
occur  when  the  source  is  approaching  the  receiver  and  positive 

^^erences  occur  when  the  source  is  moving  away  from  the 
receiver . 

Also  shown  on  the  abscissa  are  two  ranges,  both  in 
kyds  ■ (1)  R’(t),  the  horizontal  source/receiver  range  at  the 
time  the  sound  energy  leaves  the  air  source  and  (2)  R»(tA), 
the  horizontal  sc uroe/receiver  range  at  the  time  the  signal 
arrives  at  the  receiver.  Since  a signal  leaving  the  source  at 
a given  time  may  travel  to  the  receiver  by  different  propaga- 
tion modes,  with  different  propagation  times,  R'(tA)  may  differ 
between  the  modes. 

The  upper  portion  of  Fig.  3-1  displays  the  .transmission 
loss  time  history,  where  the  transmission  loss  is  referred  to 
the  air  source  level  at  1 yd  from  the  source.  To  obtain 
absolute  intensity  levels,  the  air  source  level,  at  1 yd,  must 
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be  added  to  the  TL  values  shown.  The  TL  for  each  propagation 
mode  is  shown  separately,  both  for  the  smooth  surface  conditions 
(wind  speed  = 0 lets)  and  for  the  most  probable  sea  state  (wind 
speed  » 8.8  kts).  Three  propagation  mode  types,  are  identifiable 

• single  bottom  bounce 

• convergence  zone 

• double  bottom  bounce. 

The  direct  path,  if  it  occurs,  will  be  important  only  near 
the  time  origin  (near  CPA  time)  and,  although  always  calculated 
by  the  program,  should  be  ignored  for  horizontal  ranges  beyond 
several  kyds  since  for  ranges  greater  than  this,  the  direct 
pa.th  mode  cannot  exist  physically.  When  the  direct  path  does 
not  exist,  bottom  bounce  and/or  surface  duct  modes  (if  a duct 
is  present)  must  be  considered.  Note  that  although  the  direct 
path  mode  does  not  actually  exist,  we  have  still  used  the 
arrival  time  of  the  direct  path  from  the  CPA  point  as  an 
arbitrary  time  origin.  One  could  just  as  well  plot  all  results 
versus  the  actual  arrival  time. 

For  the  summer  profile,  a surface  duct  is  not  present, 
and  bottom  bounce  modes  will  be  the  only  available  propagation 
paths  from  times  near  CPA  time  out  to  the  times  corresponding 
to  the  convergence  zone  arrival.  In  mostvcases,  the  single 
bottom  bounce  mode  will  dominate  until  the  convergence  zone 
arrival  appears.  However,  for  this  example,  the  double  bottom 
bounce  arrival  is  significant  even  within  the  convergence  zone 
because  the  bottom  loss  is  almost  negligible  for  the  grazing 
angles  associated  with  this  arrival.  While  the  bottom  loss 
for  the  single  bottom  bounce  arrival  is  also  small  for  arrivals 
originating  near  the  convergence  zone,  the  spreading  loss 
-ncreases  rapidly  as  the  range  approaches  the  convergence  zone 
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range.  Thus,  in  this  example,  the  single  bottom  bounce  arrival 
becomes  weaker  than  the  double  bottom  bounce  arrival  as  the 
inner  edge  of  the  convergence  zone  is  approached.  The  single 
bounce  signal  eventually  disappears  completely  within  the 
convergence  zone. 


In  this  example,  the  sound  speed  excess  is  very  small 
(c^  almost  equal  to  cg).  Consequently,  the  convergence  zone 
width  is  very  narrow  and  the  associated  transmission  loss 
relatively  large,  as  a result  of  the  fact  that  only  a narrow 
tube  of  rays  emitted  from  the  virtual  surface  source  c in  reach 
the  zone.  If  the  sound  speed  excess  w eve  greater,  the  zone 
width  would  be  broader  and  the  TL  would  be  less  (in  magnitude). 

Beyond  the  convergence  zone,  the  single  bottom  bounce 
mode  does  not  exist  and  double  bottom  bounce  arrivals  will 
generally  dominate.  In  the  current  example,  the  signal  level 
of  the  double  bounce  mode  is  rather  high  because  the  bottom 
losses  are  negligible. 

It  may  be  observed  that  there  is  an  asymmetry  in  the  TL 
time  histories  about  the  time  origin,  which  becomes  more 
pronounced  for  larger  times.  This  effect  is  a consequence  of 
the  ray  travel  time. 

The  time  history  of  the  Doppler  shift  (ratio  of  received 
frequency  to  source  frequency)  is  shown  below  the  TL  time 
history.  The  shift  will,  in  theory,  be  different  for  each 
propagation  mode.  For  large  times,  however,  the  shifts  for 
all  modes  will  be  essentially  the  same  because  of  the  large 
distances  involved,  as  the  example  illustrates.  As  the  CPA 
point  is  approached,  the  shifts  usually  differ.  Generally, 
the  bottom  bounce  shift  will  be  more  gradual  than  the  direct 
path  shift,  if  a direct  path  exists.  This  is  a consequence 
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of  the  fact  that  the  effective  source  for  the  bottom  bounce 
arrivals  lies  at  one  ocean  depth  below  the  ocean  bottom  for 
the  single  bounce  mode  (three  ocean  depths  below  the  bottom 
for  the  double  bounce  arrivals). 

The  D/E  arrival  angle  time  history  (Pig.  3-2)  displays 
the  angles  mode  by  mode . Positive  angles  correspond  to  surface- 
reflected  arrivals,  negative  to  direct  arrivals. 

Azimuthal  arrival  angle  time  histories  (bottom  of  Pig.  3-2) 
generally  display  very  gradual  changes  for  large  times  and  a 
rather  rapid  transition  near  the  CPA  tine.  Because  of  ray  travel 
time  differences,  the  time  histories  near  the  time  origin  may 
be  mode-dependent . 

3.1.2  Winter 

Results  for  the  ’winter  profile  are  displayed  in  Pigs.  3-3 
and  3-J». 

The  winter  profile  has  a surface  duct  and  can  therefore 
permit  an  additional  propagation  mode,  the  surface  duct  mode, 
which  may  be  important  if  the  receiver  depth  is  less  than  the 
surface  duct  depth.  In  the  present  example,  the  receiver  is 
in  the  surface  duct,  but  the  diffractive  losses  due  to  below— 
layer  leakage  are  extremely  large.  Consequently,  the  surface 
duct  TL  is  greater  than  200  dB  (in  magnitude)  and  this  mode 
may  be  ignored.  In  general,  surface  duct  propagation  will 
be  negligible  for  low  frequencies,  as  illustrated  by  the 
example. 

The  only  essential  difference  between  the  winter  and  summer 
profile  results  arises  from  the  fact  that  the  sound  speed  excess 
Is  considerably  larger  for  the  winter  profile.  Consequently, 
the  convergence  zone  is  broader  and  the  TL  less  (in  magnitude). 
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3.2  Rough  Air-Water  Interface 

Figure  3-5  presents  the  results  of  the  algorithm  for  a 
case  in  which  rough  surface  effects  are  significant.  (These 
results  were  computed  by  hand  at  any  early  stage  of  the  program 
development.  All  the  equations  used  are  the  same  as  in  the 
programmed  algorithm;  however,  the  results  may  differ  in  some 
minor  details  from  the  program  output  because  of  calculational 
errors  in  the  hand  calculation.  A limited  set  of  spot  checks 
indicates  that,  on  the  whole,  the  results  shown  in  Fig.  3-5 
and  the  program  output  will  essentially  agree.)  In  this 
example,  the  direct  path  does  exist  close  to  the  receiver  (with- 
in several  kyd.s)  since  the  CPA  is  1/3  kyd. 

Arrival  angle  and  Doppler  time  histories  are  unchanged  by 
rough  surface  effects.  In  contrast,  the  TL  levels  may  be 
substantially  different,  although  the  general  shape  of  the  TL 
history  Is  unchanged  except  in  the  direct  path  region  (ace,  for 
example,  the  1000  ft  source  height  results).  A number  of 
effects  are  operating  simultaneously  to  produce  the  differences 
in  the  TL  curves.  These  are  summarized  in  T?ble  3-2.  (See 
Ref.  1 for  a more  detailed  discussion  of  the  table  parameters.) 
In  the  following,  each  of  the  TL  curves  will  be  considered  In 
turn,  and  the  effects  of  the  parameters  in  Table  3-2  will  be 
discussed.  The  virtual  source  plots  given  in  Ref.  1 should  be 
kept  in  mind. 

Consider  first  the  50  ft  source  height  curve,  with  the 
wind  equal  2*4.5  kts.  In  the  double  bottom  bounce  regions,  the 
emission  angles  (at  the  virtual  surface  source)  are  moderately 
small  (15°-l8°)  and  the  effects  of  the  most  probable  slope 
are  still  Important.  With  the  assumed  wind  and  source/receiver 
track  configuration,  the  receiver  sees  an  arrival  which  has 
been  emitted  parallel  to  the  wind  for  negative  times  (source 
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FREQUENCY,  HIGH  SEA  STATE 
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TABLE  3-2 


FOR  GIVEN  SEA  STATE, TRANSMISSION  GAIN  RELATIVE  TO 
A SMOOTH  INTERFACE  DEPENDS  ON: 

0 - Depression  angle  of  ray  leaving  virtual  source 

^ “ Angle  between  plane  of  acoustic  path  and 
direction  of  waves 

f - Source  frequency 

h - Source  height 

c - RMS  slope  (effective) 

L - Correlation  length  of  surface  slopes  (effective) 
IT  - Average  number  of  specular  paths 


i|>  - Most  probable  slope 
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approaching  receiver)  while  the  arrival  for  positive  times 
(source  moving  away  from  the  receiver)  was  emitted  anti-parallel 
to  the  wind.  The  intensity  of  the  outgoing  ray  at  the  virtual 
source  is  consequently  somewhat  greater  than  the  flat  surface 
intensity  for  negative  times  and  somewhat  smaller  for  positive 
times.  The  average  number  of  specular  paths  is  about  unity  for 
this  low  source  height.  Thus  the  TL  lies  slightly  above  the 
flat  surface  TL  for  negative  times  and  slightly , below  for  positive 
times . 

In  the  single  bottom  bounce  domain,  the  emission  angles 
are  moderately  smal^.  .for  large  negative  or  positive  times  and 
increase  as  the  source/reneiver  range  decreases  (smaller  negative/ 
positive  times).  For  large  negative/positive  times,  the  same 
effects  operate  as  in  the  double  bottom  bounce  case.  As  the 
time  origin  is  approached,  the  emission  angles  get  very  steep 
and  the  rough  surface  effects  diminish.  Thus  the  TL  approaches 
flat  surface  values  as  the  source/receiver  range  diminishes. 

Going  back  to  the  convergence  zone  regions,  it  can  be  seen 
from  the  D/E  angle  time  histories  that  emission  angles  from  0° 
to  13°  are  involved.  All  the  energy  emitted  by  the  source  in 
this  angular  range  appears  in  the  zone.  For  large  negative  times, 
the  received  arrivals  are  emitted  parallel  to  the  wind  and  hence 
a modest  gain,  over  the  smooth  surface,  is  seen;  for  large  posi- 
tive times,  the  arrivals  are  emitted  anti-parallel  to  the  wind 
and  a modest  loss,  relative  to  smooth  surface,  is  the  result. 

As  In  the  bottom  bounce  case,  the  average  number  of  specular 
paths  is  about  unity. 

The  direct  path  region  displays  the  effects  of  the  most 
probable  slope  most  strongly.  Over  most  of  the  direct  path  region, 
only  very  small  emission  angles  are  involved.  Furthermore,  the 
angle  between  the  plane  of  the  acoustic  path  and  the  wind  changes 
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rapialy . At  the  outer  edge  of  the  region  for  negative  times, 
the  receiver  sees  rays  emitted  at  very  small  depression  angles 
in  a direction  almost  parallel  to  the  wind.  Consequently,  there 
is  a large  gain  relative  to  the  smooth  interface  results.  At 
the  other  edge  for  positive  times,  this  gain  is  substantially 
diminished  because  the  arrivals  are  emitted  in  the  direction 
anti-parallel  to  the  wind.  As  the  CPA  point  is  approached, 
the  received  rays  are  emitted  more  and  more  in  a direction 
perpendicular  to  the  wind  and  the  effect  of  the  most  probable 
slope  becomes  increasingly  less  pronounced.  In  addition,  the 
emission  angle  increases.  Hence  only  slight  gains  are  predicted. 

When  the  source  height  is  increased,  the  effective  rms 
tlope  decreases  somewhat,  and  in  the  absence  of  any  other  effects, 
the  gains  would  actually  decrease  relative  to  the  50  ft  case. 
However,  the  average  number  of  specular  paths  goes  up  substan- 
tially and  results  in  generally  large  gains,  relative  to  the 
smooth  surface,  for  all  times*.  This  effect  is  most  pronounced 
for  negative  times,  where  the  received  rays  are  emitted  parallel 
to  the  wind,  and  the  favorable  disposition  of  the  most  probable 
slope  further  enhances  the  gains  due  to  the  large  N.  Fox1  positive 
times,  where  the  received  rays  are  emitted  anti-parallel  to  the 
wind,  the  effect  of  the  multiple  specular  paths  is  diminished 
by  the  unfavorable  orientation  of  the  most  probable  slope  in 
regions  where  the  emission  angles  are  small  (direct  path  and 
convergence  zone),  and  only  slightly  affected  in  the  bottom  bounce 
regions,  where  the  emission  angles  are  larger. 


The  formalism  used  to  calculate  the  average  number  of  specular 
paths  is  probably  not  correct  when  the  Fresnel  zone  size  is  on 

°r  2reater  than  surface  slope  correlation  length, 
ll  this  is  the  case  for  the  parameters  of  Fig.  3-5(a),  then  the 
decrease  in  transmission  loss  for  the  1000  ft  source  height  may 
not  be  correct . The  resolution  of  this  problem  requires  further 
investigation  of  the  model.  A study  is  currently  in  progress. 
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For  very  large  source  heights,  the  effective  rms  slope 
would  become  very  small,  and  the  average  number  of  specular 
paths  again  approaches  unity  since  the  Fresnel  zone  size  in- 
creases faster  than  the  correlation  length  of  surface  slopes. 

The  surface,  in  effect,  again  looks  smooth  and  the  TL  should 
approach  the  smooth  surface  results.*  Physically,  the  most 
probable  slope  should  also  vanish  for  large  heights  (large 
Fresnel  zone  size);  however,  in  the  algorithm,  the  most  probable 
slope  depends  only  on  wind  speed  and  hence  some  residual  effects 
for  small  emission  angles  will  still  be  predicted.  The  correr- 
tion  of  this  artifact  must  await  a suitable  theory  for  the  most 
probable  slope. 


Provided  that  the  acoustic  Rayleigh  parameter  for  the  sea 
surface  is  small.  Further  development  of  the  model  to  account 
underwa^er^n^  as  °^^ose^  ^o  sPe°ular  transmission)  is  currently 
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APPENDIX  A 

SURFACE  DUCT  AND  CONVERGENCE  ZONE  SUBROUTINE 
Introduction  and  Summary 

The  program  described  below  was  written  to  calculate  the 
transmission  loss  (and  some  other  pertinent  variables)  between 
a moving  source  above  the  surface  of  the  water  and  a moving 
receiver  lying  beneath  the  surface.  The  program  calculates  (1) 
the  transmission  loss  for  sound  rays  which  travel  in  a surface 
duct,  if  one  is  present,  when  the  receiver  is  in  the  duct,  and 
(2)  the  transmission  loss  for  sound  rays  which  travel  to  the 
first  convergence  zone  when  the  receiver  is  within  the  conver- 
gence zone.  For  both • cases , the  program  also  computes  arrival 
angles  and  the  Doppler  shift  of  the  arrival. 

The  program  user  supplies  the  parameters  which  describe 
the  positions  and  velocities  of  the  source  and  receiver  at  the 
time  origin  ("initial”)  as  well  as  the  conditions  of  the  air, 
water  and  surface  between  them.  The  program  then  determines 
the  propagation  modes  which  apply  and  calculates  the  appropriate 
transmission  loss  and  other  transmission  variables.  A number 
of  outputs  are  printed  for  the  user.  These  include,  besides 
the  TL,  various  arrival  times,  the  depression/elevation  angle, 
azimuthal  angle,  Doppler  shift,  and  horizontal  and  slant  ranges 
for  the  pertinent  propagation  mode.  If  the  surface  is  rough, 
both  smooth  and  rough  surface  results  are  given.  The  program 
is  driven  by  supplying  time  as  the  continuing  parameter.  It  Is 
possible  to  have  certain  spans  of  time  during  which  some  propa- 
gation modes  are  not  possible.  In  this  case  the  transmission 
loss  must  be  considered  to  be  Infinite  and  no  results  are 
printed.  This  situation  arises  when  (1)  there  is  no  surface 
duct,  or  (2)  there  is  a surface  duct  but  the  range  is  too  close 
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(<5280  feet)  for  the  surface  duct  algorithm  used  to  be  valid, 
and  (3)  the  source  and  receiver  are  too  close  or  too  far  for 
convergence  zone  propagation. 


A.l  Glossary 
Input 


Wind : 

- speed  in  knots  (WS) 

- x,y  direction  cosines;  -1  < 0 < 1 (XWC,YWC) 

Air  source: 

- x,y  components  of  velocity  in  kts  (VAX,VAY) 

- initial  x,y  position  (at  t = 0)  in  ft  (XIA,YIA) 

- source  height  in  feet  (H) 

- frequency  in  Hz  (P) 

Receiver: 

- velocity  along  x-axis  in  kts  (VS) 

- initial  x position  (at  t=0)  in  ft  (XIS) 

- depth  in  feet  (D) 

(y  coordinate  of  receiver  always  assumed  to  be  zero) 

Sound  velocity  profile: 

- depth's  in  ft  (0,  DS,  DMIN,  DB) 

- velocities  in  ft/sec  (CS,  CMAX,  CMIN,  CB) 

Sound  velocities  (ft/sec): 

CS  speed  at  surface 

- CMAX  speed  at  SVP  maximum 

- CMIN  speed  at  SVP  minimum 

- CB  speed  at  ocean  bottom 

Depths  (ft): 

DS  depth  of  surface  duct 

- DMIN  depth  of  SVP  minimum 

DB  depth  to  ocean  bottom 


Time : 

- starting,  ending  and  step  size  in  seconds 
(TI,  TP,  DT) 

Questions : 

- is  all  input  data  already  in  computer?  (DATAIN) 

Y(es)  or  N(o) 
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- is  air  temperature  greater  than  50°P?  (QTEMP) 

Y(es)  or  N(o) 

- is  there  a surface  duct?  (LANS) 

Y(es)  or  N(o) 

- rough  or  smooth  surface?  (RSS) 

R or  S 


Output ; 


Propagation  mode : 

- SD  surface  duct,  smooth  surface 

- SDR  surface  duct,  rough  surface 

- CZS  convergence  zone  with  surface  duct,  smooth 

surface 

- CZSR  convergence  zone  with  surface  duct,  rough 

surface 

CZ  convergence  zone,  no  surface  duct,  smooth 
surface 

CZR  convergence  zone,  no  surface  duct,  rough 
surface 

Times  (in  secs): 

- *»  V V^MAX-  VT'  TA>  TDIP>  T0> 

t * time  of  emission 

t^  = arrival  time  of  propagated  energy 

tA“tAMAX  * dif,f,erence  in  arrival  times  between 
arrival  at  t^  and  arrival  time  of 
direct  path  from  CPA 

tQ  * time  of  CPA 

Ranges : 

- R',  R (in  ft  - RP,  R;  in  kyds  - RP3,  R3) 

- R't_t  (in  kyd  - RPA3) 

A 

- Rq  (in  ft  - RO) 

Doppler  shift: 

- dimensionless  (DOP) 

Arrival  angles : 

- D/E  (depression/elevation,  in  deg  - DE) 

- $ (azimuthal,  in  deg  - PHI) 

Transmission  loss: 

- squared  pressure  ratio  re  air  source  level  (pj) 

at  1 ft,  smooth  or  rough  surfaces  (w/o  atmospheric 
attenuation)  (PRSQ,  PRS) 
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- in  dB,  smooth  or  rough  surfaces  (w/o  atmospheric 
attenuation)  (TL,  TLRS) 

- in  dB  re  air  source  level  at  1 yd  (including 
atmospheric  attenuation)  (TL3) 

A. 2 Description 

The  program  consists  of  a main  section  plus  numerous 
subroutines.  The  main  section  calculates  some  parameters  and 
directs  the  flow  of  the  calculations.  The  subroutines  are  used 
to  input  data,  evaluate  certain  functions,  and  to  calculate 
the  transmission  losses.  Table  A-l  gives  a general 
outline  of  these  subroutines.  Details  of  the  calculations 
performed  in  the  main  program  and  the  subroutines  which  calcu- 
late transmission  loss  follow. 

The  following  sections  give  further  details  on  the  calcu- 
lations performed  by  the  main  program  and  its  subroutines. 

A. 2.1  Main  program 

The  first  event  in  the  main  program  is  inputting  the  data. 
This  may  be  done  within  the  program  or  partially  done  by  using 
a BLOCK  DATA  subroutine.  In  either  case  the  following  values 
will  be  input:  C1,  C2,  n,  wind  speed,  and  n . for  the  wind, 

vax»  vay*  xia»  yia»  h»  vs»  xis»  D»  Cs*  Cb*  V Cmin»  and  dmin* 

The  program  also  demands  some  non-numerical  input.  It 
must  be  told  whether  or  not  the  air  temperature  is  greater 
than  50°P,  if  a surface  duct  exists  and  if  the  surface  of  the 
sea  is  rough  or  smooth.  If  a surface  duct  exists,  values  must 
be  supplied  for  Cmax  and  dg.  Finally,  values  for  the  range  of 
time  to  be  covered,  t±,  , At,  are  input  along  with  the  fre- 

quency f.  The  following  values  are  now  calculated  by  the 
program : 
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TABLE  A -1 


NAME  • DESCRIPTION 

Subroutines  which  calculate  transmission  loss 


SD 


SDR 

CZS 


CZSR 


CZ 


CZR 


Calculates  ta,  ta-tAMAX>  <f>»  D/E,  doppler,  slant  range 
and  horizontal  range, pressure2  ratio;  and  transmission 
loss  for  sound  traveling  in  the  surface  duct  when  the 
sea  surface  is  smooth. 

Calculates  the  corrected  values  of  pressure2  ratio  and 
transmission  loss  for  sound  traveling  in  the  surface 
duct  when  the  sea  surface  is  rough. 

Calculates  ta,  ta-tamax,  <j>,  doppler,  slant  range  and 
horizontal  range  and  location  of  the  convergence  zone 
in  the  presence  of  £ surface  duct  for  a smooth  sur- 
face. If  the  receiver  lies  within  the  convergence 
zone,  the  subroutine  also  calculates  D/E,  the  pressure2 
ratios  and  the  transmission  loss. 

Calculates  the  corrected  values  of  pressure2  ratio  and 
transmission  ratio  when  the  surface  is  rough  and  there 
Is  a surface  duct  and  the  receiver  lies  within  the 
convergence  zone. 

Calculates  ta,  ta-tamax,<}>,  doppler,  slant  range  and 
horizontal  range,  and  the  location  of  the  convergence 
zone  when  the  sea  surface  is  smooth  and  there  Is  no 
surface  duct.  If  the  receiver  lies  within  the  con- 
vergence zone,  the  subroutine  also  calculates  D/E, 
the  pressure2  ratio  and  the  transmission  loss. 

Calculates  the  corrected  values  of  pressure2  ratio 
and  transmission  loss  when  the  surface  js  rough,  there 
is  no  surface  duct  and  the  receiver  lies  in  the  con- 
vergence zone. 


WOUT  Calculates  corrected  values  of  transmission  loss  by 

including  the  effects  of  the  air  path.  This  subroutine 
also  converts  ranges  to  kyds  and  gives-  the  transmission 
loss  re  1 yd. 
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Table  A-i  continued. 


NAME 


DESCRIPTION 


Subroutines  which  evaluate  functions 

INT  Evaluates  the  function  1(0.  a <j>  ) , 

■A  w W 

I ■ - T I> *T  **f 


S VIP 


y.W«r*  5 - VT  j A,  ■ ^ , 

A»sTc,^) 


ERF 


Evaluates  the  error  function. 


Subroutines  which  Input  data 
WIND 


AIRV 

AIRC 

AIRH 

AIRF 

SUBV 


Inputs  wind  speed  In  knots  and  the  dimensionless 
x and  y direction  cosines  of  the  wind. 

Inputs  the  aircraft  velocity  as  vectors  In  the  x 
and  y direction*.  The  velocities  are  input  In 
knots  and  converted  to  ft/sec. 

Inputs  the  Initial  (i.e.,  t=0)  x and  y position  of 
the  air  source. 

Inputs  the  height  of  the  aircraft  In  ft. 

Inputs  the  frequency  of  aircraft  radiation  in  Hz. 

Inputs  the  velocity  of  the  receiver  in  the  x direction 
in  knots  and  converts  to  ft/sec.  (The  y-direction 
velocity  is  assumed  to  be  zero.) 
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Table  A 

SUBC 

SUBD 

ENVS 

ENVC 

ENVD 

ENVA 

TIM 

WATIN 


1 continued. 


Inputs  the  x position  of  the  receiver  at  t*0  in  ft. 
(The  y-position  is  assumed  to  be  zero.) 

Inputs  the  depth  of  the  receiver  in  ft. 

Inputs  the  sound  velocity  at  the  surface  of  the  sea 
in  ft/sec. 

Inputs  the  sound  velocity  at  the  bottom  in  ft/sec. 
Also  inputs  the  depth  of  the  bottom  in  ft. 

Inputs  the  minimum  sound  velocity  (deep  sound 
channel)  in  ft/sec  along  with  the  depth  at  which 
it  orcurs  in  ft . 

Inputs  the  sound  velocity  at  the  bottom  of  the 
surface  duct  in  ft/sec  (maximum  sound  velocity). 

The  depth  of  the  surface  duct  is  also  input  in  ft. 

Inputs  the  starting  and  ending  times  for  the  calcu- 
lations as  well  as  the  step  size  in  time.  All  these 
values  are  In  seconds . 

Subroutine  used  to  change  values  of  the  Input  and 
then  rerun  the  program. 
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If  the  surface  is  rough,  the  following  values  also  are  calculated: 

t • (W)‘ 


where  gp  is  the  acceleration  due  to  gravity  in  ft/sec2  and  Up 
is  the  wind  speed  in  ft/sec. 

<r„  -yb  [ '•*>  ^ VF 


where  Uc  and  gR  are  Up  and  gp  converted  to  cgs  units.  The 
program  also  calculates  -f*  . - 

(1)M  ' V t‘t ' '1 1 
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(7) 


where  h is  the  source  height  in  cgs  units. 

w * 

(6) 


2.K*'0**Uc 


(9) 


where  the  error  functions,  erf,  are  evaluated  by  subroutine 
ERF. 


All  these  parameters  are  not  functions  of  time.  The 
program  now  is  ready  to  march  through  time.  The  parameter 
t is  set  equal  to  the  starting  value  of  time  (nt  necessarily 
equal  to  zero;  may  plus  or  minus)  rnd  the  foil  g 
parameters  are  calculated: 

* %U  (10) 

(ID 

+<fi± 

* l/ ( *4.  ~ *«.)*  * 


(13) 
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v Hjvyr7r? 
fiifL 


4 a 


where  fK  is  the  frequency  in  kHz, 


-■  (£■) 


(1*0 

(15) 


(16) 


Cfc  **  Cih'iH 

"*  d'min 


(17) 


If  there  is  a surface  duct,  the  program  will  check  to  see 
If  the  receiver  Is  within  it.  If  the  receiver  is  not  in  the 
surface  duct  the  program  skips  to  the  convergence  zone  calcula- 
tion. If  the  receiver  does  lie  in  the  surface  duct,  the  SD 
subroutine  will  be  called,  and  the  calculated  values  will  be 
typed  out.  The  subroutine  WOUT  will  be  called  next  to  convert 
the  values  Just  calculated  to  other  units  and  to  include  the 
effects  of  air  losses.  If  the  sea  surface  is  rough,  the  sub- 
routine SDR  will  now  be  called  to  calculate  the  effects  of 
this  roughness  on  the  transmission  loss.  WOUT  wilx  be  called 
again  to  operate  on  this  new  transmission  loss 

The  program  next  will  go  to  the  convergence  zone  calcula- 
tion. If  there  is  a surface  duct,  subroutine  CZS  will  be 
called.  In  the  absence  of  a surface  duct,  CZ  will  be  called. 

In  either  case  the  subroutine  returns  a tag  which  indicates 
whether  or  not  the  receiver  lies  within  the  convergence  zone. 
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If  the  tag  says  no,  the  program  skips  further  calculations  and 
chooses  a new  value  of  t.  If  the  receiver  is  in  the  convergence 
zone  the  calculated  parameter  are  typed  out  and  WOUT  is  then 
called.  If  the  sea  surface  rough  the  subroutine  CZSR  or  CZR 

will  be  called,  depending  on  the  presence  or  absence  of  a 

surface  duct.  The  values  calculated  will  be  typed  out  and  WOUT 
will  again  be  called. 


Next  a new  value  of  t is  calculated  by  adding  At  to  the 
present  value.  This  value  is  checked  against  the  final  value, 
tf,  and  if  it  is  smaller  or  equal  to  it,  the  program  will  start 
calculating  again  at  the  point  where  dependence  on  t begins. 

If  the  new  value  of  t is  beyond  the  final  value  allowed  for  the 
time  parameter,  subroutine  WATIN  is  called.  This  subroutine 
allows  changes  to  be  made  in  many  of  the  input  parameters  at 
which  point  the  program  may  be  continued.  WATIN  also  allows 
the  program  to  end.  If  that  Is  the  desired  option. 


A. 2. 2 SD  subroutine 


The  SD  subroutine  calculates  the  following: 


A-  4.  -t-t.  * 


(19) 


dopp  l< 


‘-Jtihr-I 


(20) 


I - X “W 


A-ll 
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*Wt  (**•*'  " ‘Ml- 

(22) 

d/e  »J  1 

''MW 

(23) 

A.-  ' ““€] 

(24) 

vihux  SSH*  |2.7sj  dipand»nj  on  TX^Cw'in^  »n  kH) ■ 

>!«*• 

(25) 

" Cmin 

(26) 

if ’«*.(  y /o  *• } 

(27) 

TL.loh y(-&) 

(28) 

A- 2 . 3 SDR  subroutine 

This  subroutine  uses  the  values  calculated  by  SD  as  well 
as  the  original  input  to  calculate  the  following: 

*«(*»“  **)  - 

«•»(<&,)’  ^ ’H#_  (29) 
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m 


ft  . (c«ffitU«4  of  f Cf  f [ ^2  L ±2  ) + cr  Jf  / 

1 * \ t/%  CZ 


VT 


( *tn  g S 

£s!u  , -P-  I (®/e,’F.  4>«)  N>/o"  ~5?"/ 

fC  *< 


TL--  |0 


S<^*) 


(30) 


(31) 


(32^ 


where  I(D/E,‘'V*  , Ols  ,•  4^-m  ) is  evaluated  by  using  subroutine  I NT 
and  the  error  functions  are  evaluated  In  EHP. 


Jl.2.4  CZS  subroutine 

The  first  part  of  this  subroutine  performs  the  same  cal- 
culations listed  under  the  SD  subroutine,  equations  (18)  through 
(23)  and  (26).  In  addition,  the  following  are  calculated: 


(33) 


(£zY  ‘(j; 


* ( j;  * £ 


(3*0 


(35) 
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Now  the  program  can  decide  if  the  receiver  is  in  the 
convergence  zone  or  not.  if  the  following  criterion  is  not 
met,  the  tag  'EXIST'  will  be  no,  'N',  and  the  subroutine  will 
return  to  the  main  program.  If 

ft*  * « 1 4 - ft,  i ft  ' * 

(36) 


then  'EXIST'  will  be  'Y'  (yes)  and  the  following  calculations 
are  performed : 


' v L «& J 


(37) 


(38) 


f(*u-*u)'*  <*u  * a;, 

* L J 

T/e' 

A*  . 

* ‘ * /0 


7*l  * /? 


■'>  ($) 


(39) 


(^0) 


(Hi) 


(H2) 


a -in 
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A. 2. 5 CZSR  subroutine 


This  subroutine  uses  values  calculated  in  subroutine  CZS 
as  well  as  data  from  the  main  program.  The  following  values 
are  calculated: 


COS 


C$.) 


'ft.,  ( *4.*  - nj  ^ 


(43) 


N * K1 


? fe.’-e.llei-CI 


TL  * 10 


(46) 


A. 2. 6 CZ  subroutine 

xhe  first  part  of  this  subroutine  performs  the  same  calcu- 
lations listed  under  the  SD  subroutine,  equations  (18)  through 
(23).  Then  the  following  are  calculated: 

C*  “ 

ft*  (^7) 
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Now  the  subroutine  can  decide  whether  or  not  the  receiver 
is  in  the  convergence  zone.  If  it  isn't,  the  tag,  'EXIST*  is 
set  to  'N',  no,  and  control  returns  to  the  main  program.  If 
one  of  the  following  criteria  are  met: 


* * ' * °r  e)  ^ R.'  £ 


(50) 


then  'EXIST'  is  set  to  • Y ' vps  anrf  t-ho 

. 1 » yes,  and  the  following  are  calculated: 


/,  - A * ,4  -/41-  ry  - at]' 

Ul  L (v c-*y  j 

A • l/t'- per  - jCT 

Vie  R 6t  (i- 

K CK-W 

& , jj n»  £ ..-(£) 


TL 


(51) 


(52) 


(55) 


(5*0 


(55) 
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A. 2.8  WOUT  subroutine 

The  following  values  are  converted  from  ft 


to  Jcyds 
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The  value  of  R'  is  calculated  at  ta  and  then  exprecsec  in  kyds: 


i 

3«m 


(62) 


The  incoming  TL  is  now  corrected  from  dB  re  1 ft  to  dB  re  1 yd 
by  the  factor 


(63) 


and  by  a factor  depending  on  temperature: 

Corps  1.25T  K.  10°*  p « /u 


C o*r 


- 7.4gK»cT 
ibH.oS 


4-cm^)cf‘«.4t4V4.  4 So*¥ 
+tm^«.r*.4u.rc  > 5b*  f 


(6ft) 


The  resulting  transmission  loss  is: 

TL&  a T‘U+  tL^  - corr 


(65) 


j 
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A. 3 Directions  for  Running  the  Program 

The  main  program  and  all  the  subroutines  are  collected 
together  under  the  name  WATER.  It  may  be  run  in  either  of  two 
ways: 

1.  The  data  may  be  entered  by  the  user  when  running 
the  program. 

2.  The  data  may  be  set  up  in  a BLOCK  DATA  file.  When 
the  user  is  ready  to  run  the  program,  the  BLOCK  DATA 
file  is  loaded  at  the  same  time  as  WATER.  This 
option  is  useful  if  many  parameters  remain  constant 
for  different  values  of  frequency  and  time. 

An  illustration  of  the  first  method  is  shown  in  Example  1. 
The  ucer  tells  the  computer  to  EXecute  WATER  F4.  All  user  typing 
in  the  examples  has  been  underlined.  The  computer  now  will 
compile  the  program,  if  the  compiled  version  was  not  stored. 

Then  the  compiled  version  is  loaded  into  the  working  area  in  the 
computer  and  the  program  is  started. 

The  first  question  is  whether  the  data  has  been  included; 
since  it  hasn't,  the  user  would  type  N or  NO.  Next  the  user 
must  indicate  if  the  air  temperature  is  greacer  than  50°F. 

If  the  answer  Is  to  be  "yes”,  the  user  would  type  Y or  YES. 

Now  the  program  will  ask  for  several  parameters: 

• wind  speed  (kts)  and  direction  cosines  (dimensionless) 

• aircraft  velocity  vectors  (kts) 

• aircraft  (x,y)  position  and  height  above  water  (ft) 
at  t=0 

• submarine  velocity  (kts) 

• submarine  position  and  depth  (ft)  at  t=0 

The  program  now  asks  If  there  is  a surface  duct.  The  user 
says  no  by  typing  N or  NO.  In  answer  to  the  next  question,  the 
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EXAMPLE  1 


«CX  WATFR.fa 

i oapino 

i oapfr  *.*  p.orf 

A*3W  MAX  SP*  WORMS  FRFF 

rxFniTiov 


IS  DATA  IN  COMHIJIFR?  ^ 

TrvprRATllRF  FRFATfr  than  SO  OFF-  f?  Y, 

FNTFR  THF  APPROPR ! A Tr  PARAM  FTFRS  IN 
THF  PImfnisIOvs  INDIFATFp. 


WIND  SPFFp=  r.h 

X PIRFCTIOV  msjVF  OF  ».i | mi).  . 7Q 7 
Y DIRFFTION  FOSlNF  OF  w|\ip=  . -Q7 

VFl.OFITY  V/FCTOR  X - fM  kFCT  |Oim(KTS)  * Ppp. 

'•‘FLOCITY  VICTOR  Y- PI  RFC T I ON C K TS>  = £j, 

INTIAL  X-rnORHlNATF  OF  AIRCRAFT  ft>  = cu  # 
INITIAL  Y-COORPINJATF  nr  A ? RFKAFT<  FT)  = 1 POfO.* 

HFIPHT  Of  AIRCRAFT  fr.iv  rfa  SURFAGF(  F 1 > « i^POP. 

VFI.OCITY  VFCTOR  X-DIRFCTIOm  SUHCKTS)*  Jj. 

INITIAL  X-FOORPINATF  OF  SUR(FT)s  ^ • 

PFPTH  OF  Slir<(FT>=  AOO. 

IS  THFRF  a SIIRFAFF  nilFT? 

ROUGH  OR  SN001H  SURFAFF  CR  OR  S)?_S 

SIIRFACF  SOUND  SPFFp<  FT/SFO  ) r <;pSP, 

ROTTOM  SOUND  SPFFp< fT/SFF >= 

PFPTH  FOR  ROTTOM  SOUND  SPFFO<FT)=  I 


MIN.  SOUND  SPFFD(  FT/SFF)= 

PFpTH  AT  M|N.  SOUND  SPFFD(FT)r 

INITIAL  T|mf=  - l up  j> » ^* 

FINAL  T|mf-  iron. ~ 

TIMF  jNFRFMFNTSr  '>*> . 

FRFO.  OF  AIRCRAFT  HAO I AT  I Ov < HZ ) « | SO. 

• Here,  "initial''  refers  to  t*0 

t In  this  statement,  "initial"  refers  to  the  starting  time  of  the 
encounter.  A_pn 
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SMOOTH  SURFA*F»  NO  SHKrACF  piifT 
T0=  P.O#  R‘<TO)r  IPOOO.#  R(TO)r 

h*ihooo.ft,  ws=  k.rkts#  nx=  o.7|. 


ippo7., 

N'T*  0.71 


l so. ooo 


*OPF  T TA  101 F k ' 

TO  ChANOF  a hum  f'ANAMr  l fh  , r\jjrH 


1 

0 

3 

4 

s 

A 

7 

R 

9 

10 
1 1 
12 
13 

I 4 

1 s 

fntfk 


VrL . VFC  TORS 
initial  rosniov 
HFJOHT 

RAIMOTFO  frfu. 

VfCTOH 

POSITION 


A | HC.U AFi 
AIRCRAFT 
A I UCHAFT 
AIRCRAFT 
Sl'R  VF1  . 

.SI  IH  INITIAL 
Sl'R  PFPTH 

SIJRfaCF  SOI  imp  SRFrp 
ROTTOm  SOIR'D  SLEEP  ANP  PFPTH 
mJmimiim  joim'i'  SREFD  ANP  pFpin 
VAX  . SOMVp  SpFFfl  AMfl  PfrTh 
T 1 MF  OF  F VrN T S 
RUN 
STOP 
WI*P 
THF 


H 

1 H F 


nop  n/F 

AKRROPHI  A7F 


RM I R 
'M*IHFl<: 


ItA  1 1 o 


1L 


RARAm FTFRS 
AHRKOpKI ATF 


N'^HFk  = J_P 


!wl  UAL  T 1 VF*  - \ I'MH. 
FINAL  T IMF*  mop. 
TlMF  INCRFmfmTS= 


CHAN OF  PAKAMFTFR  s n 


SMOOTH  SHRFACF:  no  siikfacf  PMCT 

T0=  0.0.  R*(T0)r  IPOOO. . R<T0)r  l?007.,  Ft  , SO.OOO 

H«IOOOP.FT.  WS«  H • RM  7S»  Na*  0.71,  Nye  0.71 


MOpF 

T 

TA 

CZ 

-44S. 

- SR9 

cz 

- 44*. 

r.7 

— 4 40  . 

-SR* 

r.r 

- 4 Mil  . 

cz 

440.  * 

49S 

rz 

4 40. 

cz 

44S. 

701 

cz 

44S. 

TO  IF  H • K 

-*01.  P3PP1  4.  PX2?\  4. 

77.  AOS  77.  40 S 
- *94. 230 4? 1 .P30  4P1  . 

7 4. 807  74. HO? 
AH4.P3IW.PI  . P30  4P1  . 

7A.R07  74.H07 
4M9.P3PP1 A.P3PP1 a. 
77.40*  7 7, 40  S 


f)Op  |>/F  phI 

I .OH  I 0.44)77. 

70. 7S? 

I • OR  I 1 . 09 i 7 7. 
70.  PL>  4 

0.93S  1.09  3. 

H 3 . 4 1 | 

0.93S  0.44  3. 

H4.0S3 


I*  RATIO  il 

1 . 70  4F  — 1 3.  I pH . h 
-119.4 
l . 304F- I 3- 1 PH. H 
-119.4 
I . 30 4F-  1 3.  1 PH.  H 
-119.4 
1 . 30  4F-  | 3-  | PH.  H 
-119.4 


CHAN OF  pARAMFTFR  E j a 


FXIT 

• Here,  "initial"  refers  to  the  encounter  starting  time. 
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uner  types  S to  indicate  a smooth  sea  surface  (i.e.,  windspeeci 
■ 0)  or  R to  indicate  a rough  surface  (windspeed  ft  0).  Now 
the  program  will  ask  for  the  appropriate  sound  speed  profile 
quantities : 

• sound  speed  (ft/sec) 

• at  depth  (ft) 

Finally  the  program  asks  for  the  starting  and  final  values 
of  the  time  parameter,  the  step  size  in  tine,  and  the  aircraft 
frequency  of  radiation.  This  concludes  the  data  input  and  the 
calculations  now  proceeds.  In  the  example  illustraced,  no 
results  are  outputted  because  the  convergence  zone  is  very 
narrow  and  the  time  step  selected  was  too  coarse. 

After  going  through  all  values  of  time  input,  the  program 
comes  back  to  the  user  to  ask  if  any  changes  in  parameter  are 
desired.  If  the  user  wants  to  change  the  times  used,  12  ;i,ust 
be  typed,  and  then  new  values  for  the  times  used  are  entered. 

If  no  other  changes  are  desired,  the  program  can  be  rerun  with 
the  new  values  by  typing  13 • Some  typical  output  is  shown. 
Assuming  that  the  user  gets  satisfactory  results  and  no  more 
changes  are  desired,  the  program  may  be  stopped  by  typing  1*1. 

(An  explanation  of  the  output  format  is  given  in  Section  A. 4.) 

There  is  much  less  input  at  the  time  of  running  the  program 

when  BLOCK  DATA  is  used.  The  BLOCK  DATA  file  must  be  created  in 

the  format  of  Fig.  A-l.  Data  are  entered  in  the  file  using 

the  same  conventions  and  units  as  when  the  program  is  run  with- 
out the  file,  with  the  following  exceptions  and  additions: 

1.  source  and  receiver  velocities  are  entered  in  ft/sec. 

2.  Cl  is  the  sound  speed  in  air  (ft/sec) 
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Figure  A-l.  BLOCK  DATA  FORMAT 


C 

C 


C 


C 


C 

C 


BLOCK  DATA 

COMMON  /{"/  VAX,  VAY,  X 1 A,  Y I A,  H , F,  VS,  X J S,  D.'AXV,  VS,  X'.-T,  YWC 
COMMON  /HOI/  T»  H»  KP,  TA,  TD1  r»  TL  » OOP,  nr,  Ph  I , PKSU,  A.MOnF 
COMMON  /CC/  Cl  »C2,CS»  CB,  DR,  DS,  CM  I v,CMAX,  DM  I N , F.TA,  OS,  OB,  KO»  TO,  PI 
1 01 

COMMON /MM /M 

DATA  VAX/371. 3*/ 

DATA  VAY/O.d/ 

DATA  XI  A/0.0/ 

DATA  Y I A/ 1 21)00  • 0/ 

DATA  H/ 10000.0/ 

DATA  F/lSM.d/ 

DATA  VS/1  1 .8  I X/ 

DATA  XI  S/0. 0/ 

DATA  0/400.0/ 

DATA  AXN/O .22/ 

DATA  KS/8.8/ 

DATA  XWC/.707/ 

DATA  YWC/.  707/ 

DATA  Cl  /I  100.0/ 

data  C2/scno.n/ 

DATA  CS/4990.4/ 

DATA  CR/50so.»/ 

DATA  DR/ 15540.0/ 

DATA  DS/4L’0.t>/ 

DATA  CMIN/4K8A.5/ 

DATA  Cy,AX  / 4«Q8  . S/ 

DATA  DM  IN/3720 . 0/ 

DATA  FTA/1.0K-P7/ 

DATA  PI  /3.  1 4J  592*8/ 

DATA  M / j / 

FND 
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3.  C2  is  the  water  sound  speed  in  the  vicinity  of  the 
surface  (ft/sec).  Any  number  close  to  5000  ft/sec 
will  be  adequate . 

4.  ETA  is  a test  number,  against  which  the  computed 
value  of  the  ratio  of  the  squared  received  pressure 
to  squared  source  pressure  at  unit  distance  is 
compared  before  the  TL  is  computed.  If  the  ratio 
is  less  than  ETA,  the  TL  is  set  equal  to  -999. 

5.  M is  an  internally  used  tag  which  tells  the  computer 
to  delete  printing  out  the  full  list  of  data  questions 
when  M = 1 . 

In  Example  2 the  BLOCK  DATA  file  is  stored  as  B2,  therefore 
execution  of  the  program  is  begun  by  typing  EX  WATER,  pi* ,B2.pij . 

In  reply  to  the  question  on  whether  data  is  in  the  computer  or 
rot,  the  user  would  type  Y or  YES.  The  next  question  concerns 
the  air  temperature.  If  the  temperature  is  less  than  50°F  the 
user  would  type  N or  NO. 

The  program  now  skips  all  the  beginning  input  and  asks  if 
there  is  a surface  duct.  If  there  is,  the  user  would  type  Y or 
YES.  The  user  must  be  sure  that  this  answer  agrees  with  the 
sound  speed  profile  included  in  the  BLOCK  DATA  file.  The  type 
of  sea  surface  is  determined  and  then  the  program  skips  over 
the  entry  of  sound  speed  profile  data  directly  to  the  values  of 
time  and  frequency  to  be  used.  In  the  example,  an  error  was  made 
in  typing  in  the  time  increment.  It  was  corrected  by  typing 
(control  A)  for  each  character  to  be  deleted  and  then  supplying 
the  correct  values.  As  soon  as  this  data  is  entered,  the  calcu- 
lations begin. 
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EXAMPLE  2 

• EX  natfr.fa,  no.r/. 

L0A01  no 

LOAOFK  rK  CORF 

*♦3*  WAX  52*  WORDS  6-RFF 

fx  edition 

IS  DATA  TN  COWPOTS-R?  X 

TFMPFRATURF  PRFATfR  THAN  50  DFA  F?  £ 

IS  Thfrf  a surfacf  pijot?  x 
ROUAH  OR  SMOOTH  SIJRFACF  (R  OR  S)  ? J* 

INITIAL  TlMFs  - a A o ( 

FINAL  TIMFs  - *(’p. 

TIME  INCREMENTS* 

FREO.  OF  AIRCRAFT  RAO [ AT I ON < HZ >*  1 50. 


ROUGH  SURFACFi  SURFACE  DUCT  * 


TO*  0.0,  R ' C 1 0 ) s I2P00.,  R < T0 ) e 1 2007 • , F*  150.000 
H*10000.FT,  WS*  8.8KTS,  NX*  0.71,  NY*  0.7! 


MOOF  T 
SO  -AAO. 
SO  -**0. 
SDR  -**0. 
SOR  -**0. 
CZS  -A*0. 
CZS  -**0. 
CZSR  -AA0. 
CZSR  -AA0. 
SO  -A30. 
SO  -*30. 
> -SOR  -*30. 
jv  SOR  -*30. 
" CZS  -*30. 

CZS  -*30. 
CZSR  -530. 
.pZSR  -*30. 
SO  -*00. 
SO  -*00. 
2 SOR  -X00. 
SDK  -*00. 


TA  TOIF  R*  R 

*03.  -*1 5.237*02.237*03. 

79.20!  79.20! 


*03.  -*15.237*02.237*03. 

79.201  79.20! 


57* » -587. 22*830. 22XB3I . 

75. *|0  75. *10 


-57*.  -587.22*830.22*831. 

75. *10  75. *10 


- 5*8  • -559.21*0*0.21*0*0.*’ 
72.020  72.020 


OOP  D/F  PHI  P RA 110  1L 

1.081  3.2*177.  0.000F+00-999.0 

72.425  .999,0 

I • 000F+00  0.000F+00-999. 0 

-999.0 

1.081  A. 8*177.  8.2*3F- | 3- | 20.8 
72. 425  - I I 3.2 

I.000E+00  | .0*4F- 1?- 1 19. 7 

-112.1 

1.081  3.2*177.  0.000E+00-999.0 

*9.094  .999,0 

I • 0O0E+00  0 . 000E*00- 999.0 

-999.0 

1.08!  7«  70 1 77.  8.280E- 1 3- 1 20.8 

*9.094  -J13.2 

t . 000E+00  I .0**F- 12- 1 I 9.7 

-112.1 

1.081  3.2*177.  0.O00F^00- 999.0 

*5.7*2  -999.0 

1.000F+O0  0 • 000 E* 00- 999.0 

-999.0 
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CHANGE  PARAMETER  - 4 

FRFO.  OF  AIRCRAFT  RAOI ATI ON<HZ ) a 8000. 
CHANOF  RARAMFTFR  = | AV**  3 


ROUfiH  SHRFACF.J  SURFACF  DUCT 


T0=  0.0#  R * ( T0 ) = I 2000 * # R<  TO  > * 12007.#  F=  5000.000 

H» 1 0000 .FT#  WS=  8.8KTS#  NX  = 0.71#  NY=  0.71 


MODE  T 
SO  -660. 
SO  -660. 
SOR  -660. 
SDR  -660. 
CZS  -660. 

r.zs  -*f0. 
CZSR  -660. 
CZSR  * 660 • 
SO  -630. 
SO  -630. 
SDR  -630. 
SOR  -630. 
CZS  -630. 
CZS  -630. 
CZSR  -630. 
CZSR  -630. 
SO  - 600 . 
SO  -60O. 
SOR  -60O. 
SOR  -600. 


. TA  TOI F R'  R OOP  O/F  PHI  P RATIO  TL 

-603.  -61 5. ?3 7 60?. 237603.  |.0R|  3.26177.  6. A0OF- 29- 2H I . 9 

. 79.201  79.201  72. APS  -33a. 9 

1 .08SF+00  I • 398F- 28- 278 . S 

-331 . S 

-603.  -61  S. 237602. 237603.  1.081  A. 86177.  I • 798F-  1 5-  I A7.  S 

79.201  79.201  72. APS  -200. A 

I.08SF+00  2. SI AF- I S- I A6. 0 

- 199.0 

- S7  6 • -587. 226H30. 226831 . I.0HI  3.P6|77.  3. A06F- 28-27  A. 7 
7C  '10  75.610  69 . 09  A -327.6 

J . 085E+00  7.  A28F-28-27  I . 3 

- 32  a.  2 

-576.  -587.226830.226831.  1.081  7.  70177.  2. 379F- I 5-  I A6  . 2 

7'.. 610  75.610  69.09a  -|99. 2 

I • 08 5E*00  3. 32AF- 1 5- I aa. 8 

-197.7 

-SAR.  -559.2I606pi.pl  A06R.  |.08l  3.261  77.  1 . 8 1 7F- 27- 267  . A 

72.020  72.020  68.76?  -320. a 

1 . 085F  + 00  3. 95SF-27-26/I.0 

-317.0 


CHAN OF  PARAMFTFR  a I A 
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In  the  example  the  transmission  loss  in  the  surface  duct 
Is  given  as  -999.  This  is  not  the  actual  value,  but  is  a 
convention  used  to  show  that  the  loss  is  so  small  it  is  not 
possible  for  the  computer  to  work  with  such  small  numbers. 

At  completion  of  the  time  values  given,  the  program  asks 
if  there  are  any  parameter  changes  desired.  If  a different 
frequency  is  needed,  the  user  would  type  H and  then  give  the 
new  value.  If  no  more  changes  are  desired,  the  user  would  tell 
the  program  to  rerun  by  typing  13.  After  completion  of  this 
new  output,  the  user  would  indicate  the  end  of  the  run  by 
typing  14  in  answer  t.o  the  change  parameters  question. 
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A. 4 Output  Description  - 1 

A description  of  the  output  produced  by  the  program  is 
given  in  detail  below.  An  example  is  given  on  which  each  item 
has  been  tagged.  These  tags  are  used  in  the  following  explana- 
tion to  indicate  the  item  under  discussion. 


9 9 

(ROUGH  SURFACr)  \SURFATF  nUCT] 

. — £>a ©h , 

[TO-  P.ol  Ir'ctcd*  i poop  .I,  ^Tto)= 

1 1 ~ n*rarJ»i‘— — 'xTi  » f • " ■ mmma»  ■ m — — rsi 


►T*  I ^01)0*.  S'!],  (VS  “pf.  o k TSV  Rx  = (1  .TTl#  Iny  = "'[T.  7' 
~CSP  ’ — <&>  -Y 1 — ~<&J— 

MODF  T Tfl  TDir  R*  R DOP 


» f=  Honn.ono 


a 


<u>\ 

so 

-*00. 

so 

-*00. 

snr< 

l-*00. 

SOR 

l-*oo. 

Ow.  CT**  S5<s 

L-  jL-  SS9  .[?.  I *0*0  .[■?  I f(.)  .! 


d/f  phi 


(19<S 


S>  I £0*0.  P 


P RATIO 

'22ti. 


TL 


» »0H  ||  G.»  77.1(7.  37RF-PHf-p>7n. 


• 3 A/i* 

otlUP . I Q3P-  p 
05? 


EXAMPLE 


Output  Description  - 2 

As  soon  as  all  input  has  been  entered,  the  program  prints 
out  a heading  which  indicates  whether  the  sea  surface  is  smooth 
or  rough  (1),  and  if  there  is  a surface  duct  or  not  (2).  Next 
several  parameters  which  describe  the  data  as  well  as  some 
parameters  which  do  not  depend  on  time  are  printed,  tg  is  given 
in  seconds  (3)»  R*  and  R,  calculated  for  t=tg  are  printed  in 
ft,  (4)  and  (5).  Frequency  is  given  in  Hertz  (6),  aircraft 
height  in  ft  (?)•  Finally  the  wind  parameters  are  printed! 
wind  speed  in  kts  (8)  and  the  x-  and  y-direction  cosines,  (9) 
and  (10),  which  are  dimensionless. 


c 1 xt  cv. 
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The  program  now  types  out  a heading  for  the  time  dependent 
calculations.  Next  a line  of  data  starts  with  a tag  indicating 
the  type  of  propagation  (11).  The  names  are  SD  for  surface  duct, 
CZS  for  convergence  zone  in  the  presence  of  a surface  duct,  and 
CZ  for  convergence  zone  when  there  is  no  surface  duct.  The 
addition  of  the  suffix  R means  the  rough  sea  surface  case;  the 
absence  of  a suffix  refers  to  the  smooth  sea  surface  case.  The 
next  item  in  the  line  is  the  time  used  in  seconds  (12). 

Data  listed  beyond  these  two  items  varies  with  the  line. 

(In  every  case  a minimum  two  lines  of  data  are  output.  If 
there  is  a rou.;;h  sea  surface  an  additional  two  lines  are  output.) 

In  the  first  line  th-?  third  item  is  tA  in  seconds  (13), 
followed  by  tA~tAMAX  in  seconds  (14)  R»  in  ft  (15),  and  R in 
ft  (16).  Next  come  the  doppler  value  (17)  which  is  dimensionless, 
D/E  in  degrees  (18),  and  azimuthal  angle,  4,  in  degrees  (19). 
Finally,  the  pressure  ratio  is  printed  (re  1 ft  and  without  air 
losses)  (20)  and  then  the  transmission  loss  (21)  in  dB  re  1 ft 
(without  air  losses). 

The  third  and  fourth  items  in  the  second  line  are  R'  and  R 
in  kyds,  (22)  and  (23).  Next  the  range,  calculated  for  t - tA 
is  printed  in  kyds  (24)  and  finally  the  transmission  loss  is 
printed  (25).  This  transmission  loss  is  in  dB  re  1 yd  and  is 
corrected  to  include  air  losses. 

The  third  line  only  appears  if  the  sea  surface  is  rough. 

The  third  item  in  the  line  is  i7  (26).  Next  the  pressure  ratio 
(re  1 ft  a^d  without  air  losses)  is  printed  (27)  and  finally 
the  transmission  loss  in  dB  re  1 ft  without  air  losses  (28). 

The  fourth  line  also  appears  only  for  rough  sea  surface  con- 
ditions. The  only  term  after  the  tag  and  time  is  the  transmission 
loss  in  dB  re  1 yd,  corrected  for  air  losses  (29). 
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.TYPF  WATFR.F4 

TRANSMISSION'  LOSS  CALCULATIONS  (100341) 

COMMON  /IN/  VAX.  VAY.XIA, Y1 A.H,  F,  VS.  XI  S.  P.AXN.WS.  XVC.  YWC,  TI  . If.  PT 
COMMON  /OUT/  T.  K.  HP.  TA.TOIF,  TL»  POP,  Dr.  PH  I . PhSO.  A.  M 01  if 
COMMON  /CC/  CI.CP.CS.CH.nR.riS.CMtN.CM  .,DM|N,FTA»r.S,r.p,hO,  10, PI. 
I Gl 

COMMON /MM /m 

COMVON/RS/  SW.CPW.PSI  R,  FPN,  Cl,  PRS,  TLRS,  RN 
C 

CI*II00.0 
CP*S0OO.0 
AXNS0.PP 
FTA* I . F-3P 
PI=3. I4IS9PA5 
SPI«SORT(PI ) 

APSO«-  3P.  1 7 4*3P.  | 7 A 

SWC=I . I sr4*SPI  / < SORT(?.(l)*(‘»i:fl  •0*'*3)  > 
FRSC=P.3Pii'»4.n*SfjRT<P.ri)/(3.0i'9R|  .0) 
r 

TYPF  9R7 

9R7  FORMA TC  IS  DATA  IN  COMPUTER?  •*) 

ACCFPT  3.DATAIM 
TYPF  Q*iP 

9 *»P  FORMATC  • TFMpFRATIjRF  C.RFATFR  THAN  SO  PFR  F?  •«} 

ACCrPT  3. 0TFMp 

IF  < (PATAIN.FO. ’Y ’ ) .OH. (OATAIN.FQ. 'YFS') ) AO  TO  9RS 
MaQ 

C 

TYPF  1 

I FORMATC  • ENTFR  THF  APPR0UR1ATF  PARAMFTFRS  IN"/ 

I*  THF  PI  Mrv.SI  OWS  INDICATFf).  •,  / ) 

CALI.  WINfUWS.XWC.  YWC) 

CALL  AIRV( VAX.  VAY) 

CALL  AIRC(XIA.YIA) 

. • CAl.l.  AIRH(H) 

CALL  SIIRV(VS) 

CALL  SI  IRC  (X  IS) 

CALL  SI'RfKO) 

9RS  TYHF  p 

? (prmatc  is  thfrf  a surface  nurT?*.px,«) 


0 
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3 


*81 


4 


5 

98* 

30 

C 


crcc 

r. 

101 

c 

10S 

c 


980 

74«. 

97V 

C 

740 


ACCEPT  3.I.ANS 
FORMA T< A3) 

RSS*  »S' 

TYPF  981 

FORVAT< '♦ROI'PM  OR  SMOOTH  SURFACF  ( R OH  S)?  •?) 
ACCFPT  3.RSS 

1F<<LAMS.F0. •YF5,).OR.<LOMS.FU. ‘Y*))  AO  TO  4 
IF  <<OATAIN.FO. ,Y,).OR.<nATA|M.FQ. ‘YFS*))  PO  TO  98* 
CALL  FN'VS(CS) 

CALL  FNVCCCH.OP) 

CALL  FN«/n<CMIN»  pviN) 

CO  TO  *> 

IF  ( < pATAf FO.  * Y * ) • 08  *<DATA1N.F0. * YFS  * ) ) PO  TO  98* 
CALL  FNVS<CS> 

call  fnvcccp*  r>n  > 

CALL  FNVn<CMIN»  DM  IN') 

call  fnvacpmax, ns> 

CONTINUF 

CALL  TIM<T I, TF.pl ) 

CALL  A!RF<F) 

CONTI  MI»F 


TOa<<XlA-XIS)*<VS-WAA)-YIA*VAY)/<<<\  :>-VAX)**P)*V/AY*  AY) 

XO«VS*T0*XI S-VAX-TO-XIA 

Y0«VAY*TO*YIA 

RPO*  SO^T  <XO*xo*yo*Y0) 

RO«S0RT<RP0*HH0*r)*n) 

IF  (RSS.MF.  •«< } ■ PO  TO  740 

WSF«|.*88*WS 

AL«C 1 /F 

X I a ( 8 • 0*M*AL* APSO/ < H I *PI  * ( WSF**a>  ))•*().  PS 
CALL. 

U»WSF*30.4K 

S W> SOR T ( S VIC ♦ U*  < I .O-FRFXI  )/P.O) 

XI?«XI*XI 

FRSaFRc.o<spi*<  I •P-FHFXl  ) ♦FXP<  - X I 2)*  C 1 . 0/<  ?.P*A  IP)-  | .0) 

I /X  I ) / < It*  a 3 ) 

SMa0. 

IF  (FRS.PF.O.)  SMrH*30. 4R*S0RT(FRS) 

00M« 1 • F I R 

IF  <ARS<SM).PT« I • E- 18)  OOMa|.0/SM 

OOM  tfm  onvt  a f»OM 
CALL  FRF<0Om, frfim) 

FOOm.r.o 

IF  <OOMp.LT.*0. ) FnnvaFXP( -OOMP) 

FPMap. a. 

IF  <SM.PT. | .F-18)  FpM«(SMaSM*<F0OM^SP!*PPM*< J.O*FRFIM)/ 

1 P.fl)**2)/(P.U*PI  ) 

PSIMa<p.H*F-3)*IJ*PI/|RO.O 

TYPF  101  XI.SW.SM.PSIR 

FORMAT<  * xia  •IPFI0.3.  •#  SIP<W)t  •,  I PFI0.3,  •,  «ia  • , I pr  1 0 . 3, 

1 •»  PMRa 1PFI0.3) 

TYPF  lOS.onM.FRS.FOOM.FPM 

FPR*AT<  ’ OnMa*.  IPFIH.3.  •»  FP  S*  * » IPl“l0.3»  *#  FOOMa  • , | RF I M. 7, 
1 *,  FRNa  •, 1PFI0.3) 

IF  <<LA.MS.FQ.  ’YFS*  ) .OR.  (LANS.  fq.  *Y  •))  PO  TO  74S 
TYPF  980 

fprmat<"’  hoiipm  smrfacfi  mo  smkfapf  nncT’/) 

PO  TO  7S | 

TYPF  979 

FORMAT<  //  ’ ROIIPh  SMRFACFI  SURFACF  0MCT*) 

CO  TO  7 SI 

IF  < <I.AVS.  FO.  * YFS  • ) .OR  . <L  AMS.  FQ.  *Y  * ) ) PO  TO  7S0 


- - - . . 
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TYPF  701 

701  FORMAT!  //  • SMOOTH  St'kFACFi  NO  SliRFACF  OUCT’/) 

RO  TO  75| 

750  TYPE  70? 

70?  FORMATS  / ' SMOOTH  SI'RFACFf  SIRFACF  DICT’/) 

751  TYPF  70 A#  TO*  RPO, kO»  F 

70  A FORMA T< * TOa'FA.I.*.  h • » < TO ) « • F7.0,  • * R(TO>» 'F7.0.  •*  F«',F9.1/> 

TYPr  7A3,M,  WS,XVC.,YkC 

7A3  FORMAT<  '♦Hs'r*,pl  *FT»  WS«  * FA.  | , »XT$»  NX*'F5.P,  ',  NY«'F5.P/> 

TYPF  700 

703  FORMATdH  » 'MOpF  T TA  TDI r K * • . R OOP  P'F» 

I * PHI  P RATIO  TL  * ) 

C 

T«  TI 

700  XS>VS»T->XIS 

XA*  VAX*  T*X I A 
Y4*VAY*T«Y1A 
XSA«XS-XA 

RP«SORT( XSA*X$A*YA*YA) 

R«SORT(RP*RP+n«D) 

FK?*F*F*< I .or- A) 

a«< <0.  i »fkp>  /( i .o*fk?>*<  /jo.o*fh?)  /(  a i oo.o+pk?) > / soon. o 

IF  < RP.LT « 59pp.  C ' RO  TO  900 

THFP»AC0S(C5>/CP) 

RB«(f.A-r.MiN)  /( nn-nM  i v ) 

IF  < .NOT.  C <I.A\S.  ro.  *Y  • > .OR.  (LANS.  EO. ‘YFS’ > ) ) Rn  TO  BOO 
ir(O.OT.nS)  RO  TO  7 AO 
CALL  SP(XSA,  YA.AkP*') 

TYPF  705.MOpr,  T.  TA,  TDI  r,  kp,  k,  HOP,  L*r,  km  | , ppso,  TL 
705  FiiRMATC  IH  , I X«3»  3P* ,0»?r7,o, r*.3, r7, p, fa.0,  | pp | P . 3 , PPF * • I > 

CALL  WOHTCO  ffMp,  TI.,  RP3#  R3.  NPA3,  TI.  3) 

TYPF  953.M0rr,  T,  f:P?.  .•  2,  HP*?,  TL? 

95.1  FORmaTCJH  , IXA3,F5.0,  |?X,?r7.3,Axr7.3,  IAXFA.I) 

I F< RSS.NF , ' R ’ ) RO  TO  7fO 
CALL  SDK(XSA. YA.  AKP".*nprR> 

TYPF  97B, MOOFR, T, AN, PNS, TL  RS 
978  FORMATCIH  ,AA,  cx.O,  33X.  ?<  IPF|  0 . 3),«l|*r  A.  | ) 

CALL  wniiT<OTF«P,  TLPS.nP?.  k3»P8A3,T  3) 

TYPF  954,MOnFk.1.TL3' 

95A  FORMAT! I H . AA, fa.O, 53XFA. | ) 

C 

7*0  CALL  CZStXSA, YA, lMrp, TnrSD, RSn, hM, FXI ST) 

CCCC 

C TYPF  |PO,THFp,TNrsn,  kSd.ro 

100  FORMA  1 ( • (OZS)  TnFR*  *,F7.?,  ',  TmF  SO*  ' , r 7 . ?,  • , RSH*  *f  9.  I , *,  nHc>, 

I F9.  | ) 

IF  (FXIST.Nr.  *Y')  RO  TO  85(1 

TYPF  705,  «nnr,  T,  TA,  Tm  r , |,p.  p,  POP,  f>r , ph  I , pkSW,  TL 
CALL  WOUKQTFvp,  TL.  KPS,  S3,  RPA3,  TL  3 > 

TYPF  953,MOf»F,T.r<P3,i\3.  WA3,  TL? 

IF  ( RSS. vf. ' k • ) on  TO  hso 

CALL  CZSk(X.SA,  YA,  TH»-P,  THFSD,  8S».  KM,  MODcR  ) 

TYPF  97H.M0DfH,T,PN,PkS,  T1.RS 

CAI.I.  W0IIT<  OTFvp,  TL  MS*  RP3,  K3,  KPA.1,  TL  3 ) 

TYPr  VSA.vnprp,  T,  n s 
RO  TO  850 

r 

BOO  CALL  CZ(XSA,Ya,Thfp,rto,rm,fx|  ST) 

IF  (FXI  ST.Nr.  «y  • ) r.O  T<»  M*-0 

TYPF  705,moof,  T,  TA,  TPi  F,  Rp,  H,  DOP,  DF,  pm  I , Pr,SU,  TI 
CAI.I  VOliTtOTfMp;  TI  , Hp3«  83,  kpa.I,  TL?) 

TYPF  95?,M0l>r,  t,RK3,  k.3,  RPA3,  TL 3 

IF  (RSS.NF. 'k')  RO  TO  850 

CALL  CZR(  XSA,  YA,  TMro,  j<TO.  KB , Moprp ) 

TYPF  V 7M , MODFk, T » an, PR  S»  TL  RS 

CALI.  WOIIT(QTFMI»,  TI  MS,  HP  7,  83,  8PA3,  11.3) 
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c 

050 

c 

900 

cccc 

c 

me 

c. 


c 

c 

o 


c. 


c 


r 

cccr 

c 

ion 

r 

r 

r 


r 


typf  vsa.moofh.t.  Tun 

CONTIN1IF 

T*T*OT 

!F<T.LF.TF>  AO  TO  70.* 

TYPr  10?#  R°*  PSP.  KTO 

FORM/T<*  KR**F9.|,*,  RSP*  *.  f>.  I . *«  NTO«*i9.|> 


M«M4  | 

CAI.1  W*TlwCM,»(n> 
ir  .*0.1  j’<  on  TO  30 

rNO 


SIIPROUTJNF  SfWXSA,  YA.ARPmx 

COMMON  /|N/  VAX.  VAY.  XIA.YIA.M.c,  VS.  A I S.  P.  AXN.VS,  XWC.  YWC.  Tl  . IF  . I»T 

COMMON  /OUT/  T»K.  RP»  TA,  TOI  *.  TL.  POP.  n*.  PM| . PI<SO.  A.vopr 

common  /rex  r i . r?.  cs.  r.o.  on.  ns.  pm  in,  r -iax.  p*  in,  *ta.  <^s.  ah,  ko,  to.  pi  , 
I «l 

MOP* • * SO  * 

CON* IPO. /HI 

ta«t*<m/ci  ) ♦<m*/c?) 

Tni*«T-To«<HH-RO)/ce 

HMl»CON*Aros<C-XSA)/RH> 

VA*  SOKT < VAX*VAA*VA I * VAY ) 

CPS  I • < VAX  *X  SA-  VA  Y • Y A ) / ( VA  • H ) 

popr  SORT*  * I w(  I .o.o.p*MA*rp«l  /r?n 

pnr*sonT<  ?.n«<rMAx-rs)/CMAx> 

PF«r.ON*nnr 

SPi»<<r“Ax-r.y|\)/(pM|N-DS)  >••(!.  p/3,  o) 

SSM.9.0 

l*<WS.LT.  in.P)  S5W**.7«. 

I*  (WS.LT.7.0)  SS--A.N 

r l*r/innn.o 

A TOT*  A*<  < 7.  a A*  «*>  /(  ( * I •*  < S.  (1/3.0)  ) 

i *soi *os*ps) «psv*S'jht<  r i 'Pin  / 3000.0 

ARPr-ATOT*RI* 

ARpM*n.r 

ir  (AkP.r.T.<-*<m).p>>  akpm*  io.n**(AKH/in.n) 

PKSO*  < < <H.  0/3.0)*  AX  V*AXN*<  no***3)  ) /(  HP* US)  ) * AM'" 

P«.guo.o 

IF  (HkSO.OF.FTA)  TL*  1(1.  •ALOoin(PkJiO) 


TYP*  ioii.a.atot, akpm 

fORMATt  • <S0>  A*  'FB.  3,  • . AlOT*  *F6.  A,  *.  not  |M  mm.  •,  IPFIC.  A) 
fi*  TURN 

FNO 

sihiromtinf  crstxsA. yo. tmth, tmtsp. hsn. sn, *xi si> 

COMMON  /IN/  VAX.  VAY  ,X  1 A.  Y I A.M.  * , VS.  XI  S.  n.  AXN.  s.  xvr.»  YV.T.  1 I . IF  . m 
COMMON  /OUT/  T.N.NP,  TA.  IMF.  P . OOp,  I'*,  HM  | , PKSV.  ".  MO|tF 

common  /rr  / r. i .cif.cs.rn,  nit,  ns.rv  in,  pmax,  nv i v,  n*!.  rs,  am,  i«o,  to.  pi  . 
1 «l 

MO»F.  TiS’ 

CON* IS0.0/PI 
TA*T«(H/c i )*(kp/re> 

TPI »• T-TO*(KP-HO) /CP 
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560 

C 

*00 


n 


c 

c 

c 


c 


son 

r 

*06 


RMI*CPM*Arns<  <-XSA)/RP> 

VA*SQRT( VAX* VOX* VAY *VAY> 
ftt*S!'»<VAX*XSA-  VAY*YA)/(VA*R) 

OOPBSUkT<  ( I .P-O.fl*  VS*ASA/<r.{>*R)  )/(  I .0-?.0*VA*CPSI  /c?>  ) 
«s*<c*Ax-rs)/PS 

■ 0|  ■ < CMAX- C*  ! » /<  pH  | N-  ,-»S  ) 

OPAC«<  1 .0/01  )♦<  \ ,0/r-P) 
rPACP»«  i .o/f.s)«(  i ,0/R)  > 

RSPa?.n*CMAX*<RrAr*SORT<  ) .o«CM|  v«f;M|N/(C"iAX»CviAX>  >♦ 

I SWHT< I .P-CS*CS/<CHCX«CHAX) )/RS) 

!rS?/ Ji,^r*£S*CS/<f;,Ur‘  j ’ /r*s‘r'rAC?*SQRT<  | .0-CMAX*rMAX/<r.n*r.P)  ) ) 
IP  ((RP.RE.RSO)  .AVP.  (HP.1  f.RH))  RR  to  *00 

I*  <(KP.RP.RH).  AMD.  ( Rp  .L  P.  «.,**D)  ) RO  TP  *00 
PXIST* 

RPTURN 


PXIST* ‘V 

V*  < VS*  VAX ) * ( VS*  VAX  > * VAY*  vAY 
XY*  <X I S-XIA)*(a|S-X|A)*YIA*Y1A 
ST*  T/ABS< T) 

Ti*To*ST*r.ORT<Tn*Tn-cxy-r<R*RB>/v; 
TP«TP*ST*SOhT( TP*  TP-(XY-HSP*KSf))/V) 
TMPSD*ACn.S(CS/CMAX  ) 

ppr»  Thth*  ( T-  T I )*(THFSn-1IHFB)  '•(TP*  Tl  > 

DP*CPV*ppr 

arp-.a*rp 

ARPM«a.O 

IP  (ARP.  RT.  <-  *iOO  .0)  ) ARP'-U  1 0 . 0*  *{  ARP/ | p.p  > 
PKS0*3P.  0«AXN*axn*  ( ~HrR*0*  THPSP*  • 0>  if  VPM  / 
3 <3«0*<  TH^R*  ThcSP)*A»<s<  *»<•  „n-  HSP*  K$|>>  > 
TL*-999.0 

TP  (PRSO.PF.FTA)  TI.*  10.  *ALPPtn<PRSO> 

RPTURM 


PNP 


SIWROUTIN’P  CZ(X5A,  YA,  THm,  RTP,  RH,  Fx  I ST) 


ro^Mnlw  /nnT/VTXBVAI'i  AtYI!/',H'  r'  VJ>» **  H.  AXN. WS.XWC,  YWC . T I . TP.  PI 

Aamw nu  r--T  T»  R.  r«P#  TA#  TPf  p.  n.  » P^'p.  |»P,PH  I ,pR.S(j.  A,>ippr 

CO^MPV  /Ct*  Cl.CP.CS.CO.  PP.nS.R'M  V.rMAX.pvtlV.rTA,  RS.RP.RO,  TO,  p|. 


**Opr**p2  . 

cnv* ipp.p/pi 

Ta*  T*<M/ri  )*<i,p/r?) 

TP I P*  T»  TP* ( HP- Rn j /pp 
|»M|*rPV*ACPS(  <-XSA»/M*> 

VA* SORT (VAX* VAX* VAY*VAT  ) 

C.PS  I ■ ( VA A • X SA  • VM Y * r rt  ) / ( VA  • R ) 

p-pp*  SORT(  ( I .0-?.p«w5*xSA/(rp*K)  )X(  | .M- 

ns*(rs*r'*iM)/n*«t  v 

ftPAC*(  I •1,/PS)4(  | , p/rp) 


y.e*VA*rpsi /rp)> 


RTP* p.o*rs* rpar • svrT ( I . o-c** I v*tm \ w( cs*r  S) ) 
Rn*p.p*rp*(r.fAr*suKT(  i •«i-cmin*c'*i n/( cp*r  mm- 
I Sont( i .p-rs*rsx(rn*rM)>xRS> 

IP  < (HP.  RP.  HTP)  . Af\»p.  ( f,p  r.  up ) ) m TP  *00 

I*  <(RP.r-P.Rn)  .AV0.(RP.LP.  KIP)  ) RP  TP  *110 
rxiST**N‘ 

HP TORN 


PXIST*  *Y * 

V»< VS* VAX)*  (vs- VAX)* VAY* VAY 
XY*(XI S*X|A)*(A| S*XIA)*Y|A*YI A 
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ST«T/APS<T) 

T!«  TO*ST*S0RT<  TO*  70-  (XY-  RP*RB>  /V) 

TP«Tn*ST*SOia(Tn*TO-(XY-KTn*KTO)/V) 

nor«  THrp«( 1 .0-  CT-TI >/l  TP-T1  )) 

np«cn*i*oor 

ARP«-A*RP 

ARPMsfi.n 

fr  <ar^  .RT.(-son.o)  > Ahi'M  = !(t.n**(AKp/in<n) 
PRS0«3?.ll*AXN*AXN«THFP*THFR*AkP!>‘l/( 3.0*ARL( RP*RR- kTO*KTO) ) 
TL«-999.0 

IF  tPRSO.PF.FTA)  TL«  1O.0*Al.OP|  Pl(  PRSO) 

C 

RFTORN 

C 

FNO 

C 

SI!RR«UT!NF  SRKtXSA, YA.ARPm.mOOFR) 

c 

COMMON  /IN/  VAX*  VAY.  X I A*  Y 1 A* M • F , vS.  X I S»  0*  AXN*  VS*  XWC,  YWC . Tl  « TF  * fiT 
COMMON  /out/  T» k.  KP,  TA.  TDI r,  TL  . OOP, OF. PH  I . PKSO. A.MOpr 
COMMON  /CC/  Cl  • CP*  CS*  CP.  HP*  PS.  C«  | N,  CMAX  » OM  I f 7A*  r; s,  PP,  Kf',  TO,  PI  , 
I PI 

COVMON/RS/  SW.CPW.PS1  H»  FPN.C!  *PRS»  ILIvS.BN 

c 

MOOER* * SDR • 

CPW«(XWC*X$A-YVC*YA)/RP 
TNF«nF*3. I A| SYSAS/IRO.O 
CALL  INTCTHF) 

FC«1 .Fjp 

IF  <SU.NF.O.<!)  Ff*  | .0  / < SORTt  2.0  > *Sv’> 

CTAfiCnS(TH»V?.0» 

STA«SlN(T^r/p.m 
ARP | a STA*FC /CTA 
f RP2«rc«C | .r-STA) '(CTn) 

CALL  FRF{  A|{P | , rRF  | j 
CALL  FRF<A«PP, frf2> 

PN«  FPN* ( FRF ) * r KFp  ) 

PR$*<H.0*AXN*AXN«C1 •RN*ARPM)/< RH«nb) 

?LRS«-999.0 

IF  tPRS.PF.FTA)  TLRS*  IP.0»Al  Or.|rt(PKS> 

c 

hr TURN 

c 

FND 

C 

SI'PROUTINF  CZSRCXSA,  YA,  Thf»,  Tmf SO.  hSI*.  hP,-onrK) 

r. 

COvMrtw  /IN/  vax.vay.aIa.yia.h.f,  VS.AI  s.n.Axv.ws.x’-r.  r-«  . 1 1 . i» .,»! 

Common  /CUT/  T*K*M*  1 a*  |H|  f,  71  . oo^,  |,r,  ch  1 1 Pl/Sv,  •»,  MOf.c 

common  /cr/  n.rp,i;s,cM.rH>,ns.rM|'.fMuA,ov|N,»  1a, po,  10, r\ . 

I PI 

COmmon/RS/  SV.CPW,PSIi».*P'’.CI  . PKS.  ILLS.  P*>' 

c 

MOOFRi  ‘C^SR' 

CPV*<XWC*XtA-YVC*YA> /hP 
CALL  iVKTHFr.) 

ClOrCI 

CALL  |N7<TMFSr> 

rc*  1 . fib 

I r (Sw.Nr.n.O)  fc«1 .0/<  sOPT(P.r>«.sy> 

TMS« TMFM«1mfS0 
CtS«CnS< Ims/P.O) 

STS- SIN(IhS/P.C) 

ARP|«STS*FC/CTS 
ARP?.f,-*(  | .O-STS)  ACTS 
CALL  FhF( ARP| , i nr | > 
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c 


c 


c 


c 

c 

c 

r, 

c 

r 

c 

I? 

c 


Cali.  rKF<AKr.{>,  f«f?> 

fw«rpM*crKn  *rR(£>) 

AMNrH.H 

APfcc-A*RP 

tr<Aph'.  Rl.-SAP.lt)  Apple  1O.O»»(APK/|O.0) 

PPSOP.P*  A.\N«AXN«  CLI  n-C|  J*AKpM»f«N/(  7MS*  APS<  RP*kP-hSn«  RM»>  ) 
TLRS«-9A«».n 

!«■  (FkS.r.r.FTA)  71  MS*  U1.0«ALnn|  OtPKS) 

KFTMRV 

r\’n 

SlWpomivr  r.ZR<XSA.YA»THFn.k70,Rn,MnnFk) 

common  / | \ / VAX.  VAY.XI  A.YI  A.m.f,  VS.KIS.n.AXv.PS,  X^c#  YVC.7  I » T‘  . i 7 

CW«ON  /OUT / T.k.fP.  TA,  TP!*.  TV  .PAP.  Pr.  Ph|  , Hi-Mi,  A.moi-F 

CfHMOv  /C C/  Cl  .CP.CS.CP.  P»*.  r>S,CM  |n,C«AX»  iv>  IN,  f TA.nS.  Ah.  ho,  TO,p|  , 

! r.i 

COM«ON/KS/  SW.CPV.PMP.  FFN.C!  .PKS.  TV  KS.PN 

WOPFhe  * CZR* 
r.Pw«txwc«xsA-Y,*r«YA)  /p  ' 

CALL  INTtTHFP) 
rc*i.^iA 

ir  (sv.vr.o.o)  rc> i .n/(SbR7Cp.o)*sw) 

AR Ac 7mFm/?.0 
CT?« CC.SC  AKA) 

STP*  5IN( AHA) 

ARr.|cS7P*«T/C7P 
ARC.pl »r.«<  I .n-S7?)  /C7P 
CALL  rpr<APA| . rj,r  1 > 

CAI  l rpr(AK r.->,  rht?  > 

AN*rpN*(  rpr  I 

ARPMrf1.l1 

ApRc-A«KI> 

I r<ApR. at . - sop.m  AhP-iiP.ii^^eAPK'in.OJ 
PPS*3?.n«AXN*AXv.r  I •PN*AR pv/<  7^Ef-»AMSC  M> • KP-  K TOc  p 7")  ) 

TLRS'-p^h.o 

irCPRS.r.r.FTAl  71  K S*  | 0. 0«  AL  OC  | P<  PK5> 

hriUKN 

FNP 

SUAWUITIVF  |N7(7h) 

canmon/rs/  sw.cpp.pmp, rPN , c i » pk s, ti.ks.av 

IF  <APS(SW).at.  I .*■- |H)  no  TO  | p 
r |cC  ( 7H.PMH.CHW)  ••3-  <PSI»*«f  r»k  >«0>  /.l.p 

RCTiirv 

SCP I * SORT  < "*  • |i|'°?*M 
S0?« SOR7<P.O) 

T7*P.C/?.0 

r7cS.ii/3.0 

pc<*psi M«rpw/t so?«  sw> 

PC  SP* PCS* PCS 
ts«th/<  s'*v«sm 
AKC.c7S»PCS 

ARnp^AKCcApr, 

CALI.  FKC<  Akn,  rpriK) 

CALI  FRFO'CS.  rr.rpi 

CI*<e<S0?»Su)**3>  //>.p)«(  7 S*  77*  < AhT2*  AKC-  l*CSP«PC  S)  ♦ FRF  TP*  < I • 0*  TT 
I •ARC.P)*Awn-FKfp«C  I .♦WT«r'CS?)*PCS*CXF(-AkC?>*(n*  17*AKn?> /SOP! 


A-37 


Bolt  Beranek  and  Newman  Inc 


TM  W307 


c 

c 

c 

c 


200 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 


2000 


2 -exf<-pcs2>*<et*ti*mi:s2)/sopi  > 

RETURN 

END 

SUBROUTINE  FRF ( ARP#  FRFS) 

SARO«ARO/AB5( ARC) 

X*ABS<ARP) 

FR*I .0 

IE  (X.PT.A.2)  DO  TO  200 

F*  I .0*  . 0705230  7rl/»*X'-.0A22«2P  I 23«X*X*  . U0927C5272*<X**3)  ♦ 

I .00015201  A3*  (X**  a)  <•  .0002  7X5X72*1  X«* 5) *.0000 A3UX3K*(  X***) 

ER«|. 0-1 • 0 ' < F*  * I X) 

FRFS=FR*SARr. 

RETURN 

END 

SUBROUTINE  WOUTtUTFMP,  TLO.  RP3,  R3.  RPA3,  TL3) 

COMMON  /IN/  VAX, VAY.XIA.  YIA.H, E,  VS, X ’ S, D» AXN,  WS» XVC,  YVC 
COMMON  /OUT/  T.R.RP.TA.Tnir.TT-.rJOP.Dc. PHI, PfiSO.A, MODE 

RP3=RP/3000.0 

R3«R/30O0.O 

X3*< VS*TA  + X IS)- ( VAX  * TA  + X 1 A) 

Y3* VAY*TA*YTA 

KPA3tS0RT(X3*X3*Y3*Y3) /30O0.O 
TL3e-999. 

IE  f TLU.E0.-999. ) RETURN 
TL9*  IO.0* ALOr-1 0(9.0) 

CO  HR  a ( I • 25c- X ) * F*H 

I E(  (OTFMP.EO.  'Y').OR.<OTFMP.FO.  'YES')  )CORRa(  7.  AE-R)*F*F*H/|  x/,.05 
TL3=TLU*TL9-C0RR 

RETURN 

END 

SUBROUTINE  WATIN< INM.KA) 

COMMON  /IN/  VAX,VAY,X!A,YIA,H,f,  VS,  X I S.  U.  AXN,  US,  XV*C,  Y WC  , T I . 1 ►',!>! 
COMMON  /OUT/  T,  R,  RP»  TA,  TDI  F,  TL,  DOP.  BE,  hh  I » PKJ.U,  A,  MUI'E 
COMMON  /CC/  Cl  ,r?,cs,  CH,  OR,  ns,  CM  IN,  CM  AX,  F'min,  FTa,  rs»  Mi.  ho, 

I PI 

COMMON /MM /M 


Me  | MM 

IP(M.PT. I)  PO  TO  3000 
TYF  E 2000 

FORMAT* I H , ’TO  CHANPF  a RUN  PARAMETER, 


ENTFR  THE' 


I • APPROPRIATE  NIIMPFRs*/ 

I * I AIRCRAET  VEI. . VECTORS'/ 

I * 2 A I RCRAET  INITIAL  POSITION'/ 

I * 3 AIRCRAET  HEIGHT'/ 

I * A AIRCRAET  RADIATED  ERFO • * / 

I • S SUD  VEI..  VECTOR'/ 

I » X SUD  INITIAL  POSITION'/ 

I • 7 SUP  DEPTH'/ 

I • 8 SURFACE  SOUND  SPFFD’/ 

I * 9 BOTTOM  SOUND  SPEE.i  AND  EPTH  * / 

r • in  minimum  sound  speed  «m>  dfpth'/ 
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ptitn 

3000 

3001 

300? 

47 

49' 

49 
1 0 
1 I 
12 
n 

14 

1* 

I A 

19 

20 

18 

904 

r 

30 

r 

C 

?no 

?oi 

POP 


P 03 


200 

201 


I * II  WAX.  .SKHI)  ANIJ  limn*/ 

J * 1?  T 1 1^ nr  EVFNTS'/ 

1 » 13  RIN '/ 

1 • 14  STOP’/ 

1 * is  wimp  pakamftfrs*/ 

I » ENTER  THE  APPHOPRIATF  NUMRER  ■ »,*> 

ACCEPT  2001, *3 
rOR^AK  1 ) 

CO  TO  < 47# 48, 47.  10*  1 I , 1 P.  1 3,  I 4,  | 5, | A, | 9# 20, 3C , 18, ?R4)K3 
TYPP  3001 

FORMATS*  ChAmc.f  PAnAMF  TF k * •#») 

ACCFPT  3002# *4 
F0RWAT< I ) 

fin  TO  (47#  40,  49,  10,  II  , \p,  1 3,  I 4,  1 A,  1 A,  19*2000,  18,904>K4 

CALL  A 1 RV( VAX, VAY ) 

fiO  TO  3P00 

CAIJ.  A1RC(X  IA,  Y IA) 

PO  TO  3000 
CALL  AIRH(H) 
fiO  TO  3000 
CALI.  AIRF(F) 
fiO  TO  3000 
CALL  SMRVCLS) 
fiO  TO  3000 
CALL  SUAC(XIS) 
fiO  Tn  3000 
CALL  SliRrKP) 

PO  TO  3000 

CALL  FNVS(CS) 

fiO  TO  3000 

CALI  FNVC<CP,0FO 

OO  TO  3000 

CALL  FWO(wIh.WIM) 

•?n  to  3oeo 

CA»  I FNVA<C*AX,  OM 
fiO  TO  3000 
CALL  T|M<T1.TF.0T> 
fiO  TO  3000 
CALL  FX1T 

CALL  W|Np(WS,XV*C,YVC> 
fiO  TO  3000 


*3=  13 
RFTIIRN 

F«vn 


SUPKO 1 1 T I v «•  ’■>!  M0(  V.’S,XL’0»  Y VC  ) 

TYPr  200 

FORyATdH  ,",.l\r>  Sp^ED-  ',♦) 

ACCEPT  20 1 , NS 
EORMAT(E) 

TYPE  20? 

fokmat<'*x  ojr'ttiov  cosinf  nc  wivn=  •,«> 

ACCEPT  20 1 , X VC 
TYPE  203 

FORMAK'^Y  niRFCTIO'*  COSINE  OF  MIMAS  ',«) 

ACCEPT  20 1 , Y’.-’C 

RETHHN 

FNp 

SMHROMTI^P  A1RV( VAX, PAY ) 

TYPC  200 

FORMA  T(  I H ,'VFIOCITY  VFClOH  X- 0 1 RECTI  OM(  RTS)  ?X,  f ) 

ACC. FJ> T 20I  » VAXK  T 

FORMAT(F) 

VAX* 1 . XflR*VAXRT 
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202 


200 
20 1 
202 

200 

201 

200 

201 

200 

201 

200 

201 

200 

201 

200 

201 

200 

201 

?02 


typf  202 

format*  •♦</*•(.  nr  nr  vfctor  y-di  i>fction<ktsi*  *,2x,  «> 
accfpt  2oi«vaykt 

VAY* I • XRH«  VAYRT 
RFTl'RN 

fn  n 

SlIMROIITINF  A 1 RC  < X I A,,*  Y I A ) 

TYI'F  200 

FORMAT  ( I M « * 1 N T I A| . X-COAKDINATF  OF  AIRCRAFT*  FT)*  ',?X*  F) 
ACr.rpT  20 1 » X I A 
FORMAT*  F) 

TYPF  202 

FORMA T( •♦INITIAt  Y-COORDINATF  nf  A 1 RC RA : T* f T) * ’, 2X» f J 
ACCFPT  201  * Y 1A 
R»  TURN 
FND 

R UPROOT INF  AIRH(H) 

TYPF  200 

FORMAT* IH  t 'HF  I OHT  f)~  AIRCRAFT  FROM  SFA  SuKFACf* F 1) « • , 2X, d 
ACCFPT  201, H 
FORMAT*  F) 

RFIORV 

END 

SUBROUT  VF  A I R r C F ) 

TYPF  200 

FORMAT ( I H , *FRC0.  OF  AIRCRAFT  RADI  AT  I ON* HZ >* '* 2X,  « ) 

ACCFPT  201,  F 
FORMAT* FJ 
RFTURN 
END 

Sl'RROUTINF  SUBV*V5> 

Tvoc 

format*  IH  , M/Fl.ociTY  WECTOR  •-PIRFCTION  S"H<*fTS>*  *,  '’a,  f ) 
ACCFPT  20  I » VS*  T 
FORMAT* F) 

VSr  | . Xffft*  VSKT 

RFTCRw 

END 

SIJOROUTINF  SI)PC*XIS) 

TYPF  2C0 

FORMAT* |H  , ’INITIAL  X-COORDIVATF  OF  SUP*  rT)= '»  ?X»  * ) 

ACCFPT  20  I , X 1 S 
FORMAT* F) 

RFTORN 

FNO 

SUAROIITINF  S'lRDtp) 

TYPF  200 

format*  IH  ,’r>rH7H  Of  Sl’M*  f T)  * ’ , 2X,  % ) 

ACCFPT  201, D 
FOF-'ATt  f) 

RE TORN 
FND 

SlIRROMT  INF  FNVS*CS> 

TYPF  200 

FORMAT*  IH  , ’SURFAC.F  SOUND  SPFCD*  FT/SFC)«  ’,  -X,  * ) 

ACCFPT  20  I , C S 
format*  f) 

RFTI'RN 

FND 

SiipROPTiNF  FNvr.*cn,Dn) 

TYPF  200 

FORMAT  < ( H , ’BOTTOM  sound  SPFFD*FT/SF0* ’,  2X,  1) 

ACCFPT  20  I, CH 
FORMAT*  F) 

TYPF  202 

FORMAT* ’♦DFPTH  FOR  POTTOM  SOIWD  SPFFD* F T) a ’ , 2X, t ) 
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pan 

poi 

S»(1S> 

pop 

poi 

pop 

?«e 

poi 

POP 

P03 


ACC.PT  20l*t>H 
RFTIIRN 

run 

SMBMAtlTINF  ri\»vn(C>,IN,nMJN) 

TYPF  P00 

rORMATUH  # **1U.  SOItN'p  SPFFp< FT/SFC >* *, PX» J > 

ACCFPT  PCI  I # cm  I M 

FOH'IAT(F) 

TYPF  POP 

FORMATt  ‘♦r'FPTH  AT  «1N.  SOUVD  SPFFf)<  FT)  * *#  P7  * f ) 

ACCFPT 

RFTMR* 

F*p 

SIIARflMTI^F  FMVA<r*AX«  nS) 

TYPF  P00 

FORMAT<  I M > *VAX  • SniNf)  SPFFpf  fT/SFC>=  • , ?X#  * ) 

ACCFPT  ?0l#CMAX 

FOR>1AT<F) 

TYPF  POP 

Ff)RMAT( ‘♦PFP7H  AT  MAX.  srmvp  SP^FfH  FT) * * . ?x»  « ) 

ACCFPT  20 1 « PS 

KFTtiRiv 

EVP 

SliRRPIlT IMF  TIM(TI.TF»PTJ 
TYPF  POM 

msrvaTf  I m .mwitiai  T »mf*  ♦ , p*  . « 1 

ACCEPT  PC  I # T I 

FOKPAT<F> 

TYPF  pfi? 

FOHVAT( •♦final  T IMF*  •»  PX»  f ) 

ACCFPT  2(11  ,TF 
TYPF  PP3 

FC|RMAT(  •♦T1*!F  I JChFMFM.TS*  *#2/#  f » 

ACCFPT  POI.DT 
RFTIJHN 

rvp 


4 


I 
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APPENDIX  B 

DIRECT  PATH  AND  BOTTOM  BOUNCE  SUBROUTINE 
Introduction 

To  calculate  the  transmission  loss  as  a function  of  time 
for  a moving  source  above  the  sea  surface  to  a moving  receiver 
below  the  sea  surface  through  varying  sea  states  (smooth  to  very 
rough),  a computer  program  was  developed.  The  sections  which 
follow  will  be  concerned  with  two  modes  of  propagation  from  the 
source  to  receiver;  direct  and  bottom  bounce  rays.  The  bottom 
bounce  rays  will  be  further  subdivided  i.ito  single  bottom  bounce 
or  double  bottom  bounce  rays. 

Finally,  two  types  of  sound  speed  profiles,  one  with  a 
surface  duct  present  and  ore  without,  will  be  utilized. 

To  perform  the  calculations  for  the  different  modes  of 
propagation  involved,  15  separate  subroutines  are  used  with 
one  main  program  as  the  control.  Below  is  a listing  of  the 
main  and  subroutine  programs  used. 


MAIN  PROGRAM 
WATER. FU 


SUBROUTINES 

BOUN(THEB) 

BSSR(THE0 ) 

SBOUN(THE0) 

BDSR(THE0 ) 

BS(THE0 ) 

DP(THE0 ) 

BD(THE0 ) 

DPR(THE0 ) 

BSS(THE0 ) 

ERFS.F4 

BDS(THE0 ) 

BT3.FU 

BSR(THE0 ) 

BT5.F4 

BDR(THE0 ) 

B-l 
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Figure  E-l  is  a block  diagram  showing  the  flow  of  the  cal- 
culations. This  chart  gives  an  overall  picture  of  the  flow  of 
the  program  calculations  during  the  actual  execution. 

B«1  Brief  Description  and  Purpose  of  Program 
MAIN  PROGRAM  - WATER. F4 

The  main  program  has  several  functions.  One  of  these  is 
accepting  all  the  input  parameters  which  initialize  the  problem. 
The  input  parameters  are  divided  into  four  main  categories : 

(1)  source  parameters,  (2)  receiver  parameters,  (3)  environmental 
parameters,  and  (4)  time  parameters.  Once  all  the  parameters 
are  inputted,  a series  of  calculations  are  performed  to  obtain  a 
set  of  constants  for  each  particular  time.  These  are  stored  in 
a common  block  with  all  the  input  parameters  so  that  they  can  be 
shared  by  all  the  subroutines.  The  third  function  of  the  main 
program  is  to  route  the  flow  through  U main  loops  by  a series  of 
conditional  statements.  Finally,  the  main  program  controls  the 
output  of  values  from  each  individual  subroutine  so  that  the 
values  arc  displayed  in  a pre-assigned  format. 

BOUN(THEe) 

This  subroutine  calculates  the  angle  o,n  the  arriving  ray 
which  reaches  the  receiver  with  no  surface  duct  present.  It  first 
determines  whether  there  is  a single  or  double  bounce  which 
reaches  the  receiver.  An  iteration  routine  is  then  used  to 
finally  calculate  the  acquired  ray  if  there  is  one.  (NOTE:  if 
no  single  bounce  exists,  the  program  defaults  to  double  bounce; 
if  the  double  bounce  does  not  exist,  the  control  is  sent  back  to 
main  program  where  DP.F^  is  called.)  Now  depending  on  whether 
the  solution  is  a single  or  double  bounce,  the.  main  control  is 
either  sent  to  BS(THE0 ) or  BD(THE6 ) . 
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SBOUN(THE8 ) 

The  difference  between  3BOUN(THE0)  and  BOUN(THE8 ) is  that 
a surface  duct  is  present  for  the  calculation  of  the  angle  of 
the  arriving  ray  which  reaches  the  receiver.  Again,  it  also  has 
to  determine  first  whether  there  is  a single  or  double  bounce 
which  reaches  the  receiver  with  the  surface  duct  present.  An 
iteration  routine  is  then  used  to  finally  calculate  the  acquired 
ray,  if  there  is  one.  Depending  on  whether  the  solution  is  a 
single  or  double  bounce,  the  main  control  is  sent  either  to 
BSS(THE8 ) or  BDS(THE0 ) . 

BS(THE3 ) 

The  single  bounce  rays  with  no  surface  duct  present,  are 
handled  within  this  subroutine.  Besides  calculations  to  obtain 
the  TL  for  this  particular  mode  of  propagation,  other  important 
quantities  of  interest  are  also  calculated.  The  final  output 
values  from  this  subroutine  are  sent  to  the  main  program  where 
they  are  printed  out  using  the  main  program  format. 

BD(THE8 ) 

This  is  the  double  bounce  subroutine  which  performs  all  the 
calculations  for  double  bounce  rays  with  no  surface  duct  present. 
Many  of  the  formula  are  modified  versions  of  the  single  bounce 
equations.  The  final  output  values  from  this  subroutine  are 
again  sent  to  the  main  program  where  they  are  printed  out. 

BSS(THE8 ) 

For  bottom  bounce  rays  which  are  in  the  presence  of  a surface 
duct,  the  calculations  are  performed  within  this  subroutine. 

The  final  output  values  are  returned  to  the  main  program  where 
they  are  printed  out. 


M 
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BPS (THE 9) 

This  is  the  subroutine  which  calculates  double  bounce  rays 
propagation  with  a surface  duct  present. 

BSR(THEe) 

All  subroutines  with  the  letter  R contained  in  the  name  are 
used  for  rough  surface  calculations.  This  subroutine  modifies 
the  smooth  surface  results  for  the  single  bounce  mode  of  propaga- 
tion. It  has  the  necessary  equations  to  calculate  TL  for  a rough 
sea  surface.  Once  all  the  output  values  are  calculated,  they 
are  sent  to  the  main  control  program  to  be  printed  out.  No 
surface  duct  is  present . 

BDR(THEe) 

This  is  the  subroutine  which  handles  rough  surface  double 
bounce  propagation  (no  surface  duct  present).  Again  it  returns 
to  main  program  with  a new  PRATIO  and  TL  for  this  mode  of  propa- 
gation in  a rough  sea. 

BSSR(THEO) 

To  take  into  account  a surface  duct  in  a rough  sea, subroutine 
BSSR(THE9)  handles  the  single  bottom  bounce  mode  of  propagation 
for  this  condition.  It  calculates  a new  TKATIO  and  TL  due  to  the 
change  in  sea  surface  conditions.  It  returns  to  main  program 
with  output  values. 

BDSR(THE9 ) 

This  subroutine  calculates  the  double  bounce  ray  for  a rough 
sea  surface  and  surface  duct  present.  It  modifies  the  smooth 
surface  formula  and  calculates  a new  PRATIO  and  TL  due  to  the 
rough  sea  state.  It  returns  to  main  program  with  output  values. 
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DP(THE9) 

For  direct  path  propagation  all  calculations  are  handled 
within  this  subroutine  for  smooth  surface  conditions. 

DPR(THE8 ) 

The  direct  path  propagation  TL  wich  a rough  sea  surface  are 
handled  by  the  formula  which  are  modified  to  take  into  account 
the  rough  sea  surface. 

ERFS.Ffr 

To  obtain  the  statistics  for  the  rough  surface  subroutines, 
the  error  function  is  necessary.  This  subroutine  calculates 
the  error  function  for  a given  input. 

E>r3.F*4  and  BTo.Pfr 

These  two  subroutines  are  used  to  obtain  the  bottom  reflec- 
tion loss  as  a function  of  frequency  and  grazing  angle.  Bottom 
types  3 or  5 are  typical  of  most  of  the  ocean  bottoms  of  the 
world.  This  data  was  obtained  from  the  FACT  model. 

B.2  Instructions  for  Execution 

To  execute  the  calculations,  6 programs  must  be  loaded. 

They  are  WATER. FH,  BO.FH,  SB0.F4,  BT3.FH,  BT5.FU,  and  ERFS.F**. 
Once  loaded,  the  main  program  starts  the  control  which  asks  the 
user  to  input  all  the  parameters  necessary  for  the  calculation. 
The  program  asks  for  the  following  Inputs  listed  below:  (The 

Inputs  must  be  floating  point  numbers.) 
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INPUT 

UNITS 

COMPUTER 
PGM  NAME 

Aircraft  velocity  x-directlon 

knots 

VAX 

Aircraft  velocity  y-direction 

knots 

VAY 

Initial  x-coordinate  of  aircraft  (at 

t«0)  feet  . 

XIA 

Initial  y-coordinate  of  aircraft  (at 

t*0)  feet 

YIA 

Height  of  aircraft  from  sea  surface 

feet 

H 

Frequency  of  aircraft  radiation 

He 

F 

Receiver  velocity  x-direction 

knots 

VS 

Initial  x-coordinate  of  receiver  (at 

t*0)  feet 

XI  .s 

Depth  of  receiver 

feet 

D 

Wind  speed 

knots 

V 

x-direction  cosine  of  wind  direction 

dimensionless 

wsx 

y-direction  cosine  of  wind  direction 

dimensionless 

WSY 

Bottom  type  (3  or  5) 

dimensionless 

BT 

Starting  time 

seconds 

T I 

Final  time 

seconds 

TF 

Time  Increments 

seconds 

DT 

la  there  a surface  duct 

yes/no 

Surface  sound  speed 

ft/sec 

CS 

Bottom  sound  speed 

ft,  sec 

~B 

Depth  for  bottom  sound  speed 

ft 

DB 

Minimum  sound  speed 

f t/sec 

CMI*' 

Depth  of  minimum  sound  speed 

ft 

& 

Maximum  sound  speed 

r./sec 

CMAX 

Depth  at  maximum  sound  speed 

ft 

DS 

The  output  is  formatted  such  that  a heading  is  pr’  ted  first 
with  the  following  names:  MODE,  T,  TA,  TDl'F,  ?.*,  R,  D/E, 

PHI,  PRATTO,  TL. 

MODE  - Mode  of  propagation 

T - Start  time  of  ray 

TA  - Time  of  arrival  of  ray 
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difference  between  TA  and  time  of  arrival  of 
direct  path  ray  from  CPA  point 

Lateral  range  bet  teen  source  and  receiver 

Slant  range  between  source  and  receiver 

Doppler 

Arrival  angle  at  receiver 
Azimuthal  arrival  angle 
Received  mean  square  pressure 
Transmission  loss  betv/een  source  and  receiver 

The  corresponding  value  for  each  heading  is  given  in  a 
column  below  the  headings. 

Once  the  set  of  calculations  have  been  performed  the  program 
prints  out  a message  to  ask  the  user  whether  he  would  like  to  run 
the  program  again  with  any  of  the  input  parameters  changed.  The 
user  may  opt  to  run  with  a new  set  of  parameters  or  stop  the 
execution  totally  by  typing  the  pre-assigned  number. 

The  following  pages  are  a listing  of  a sample  run  showing 
how  the  input  parameters  are  given  and  the  actual  formatted  output. 


TDIP 

R» 

R 

DOP 

D/E 

PHI 

PRAT 10  - 

TL 
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- n: ; iinnrr.r  i.G0.r4.s?0.r4OT3.F4.  Ts.r4.EFrs.r4 

rirrrivn  EPrs.r4 

LUlIJiHG 

loader  in  rope 

13*3»  nr  : 47?  wor-i'S  ffee  ikput 

E:rCUTIOM  PARAMETERS 


enter  tmc  RrrporriRTE  w pneteps  *n 

THE  DIMENSIONS  INDICATED. 

UCLOCITV  "ECTOP  N-DIFECT ION N TS«  * £M.O 
MELOCITY  "CCTOP  Y-niFECTIONH  TS1  = 0.0 

IHTIflL  : -COORDINATE  or  AIRCRAFT i FT  l*  0.0* 
INITIAL  Y-COOPBIMnTE  Of  AIRCRAFT  (FT  i * 1000.0* 

height  or  aifcfaft  rr on  sc  a soprACEim*  so.o 

FPEO.  or  AIFCPhTT  r hIiIm!  IOHiHZ I = 1000.0 

1 ■CLOCITV  "CCTOF  : !— I* I F EC  T I 'tl  SUB  U TS  ' “ 10.0 

INITIAL  .-coordinate  or  SUB i rT  I = 0.0* 

depth  or  suB'FTi=  300.0 
MIND  SPEE3N  TSi-  04.  G 


X-DIPCCT ION  COSINE*  1.0 

Y-DJPCCTION  COS  I ME*  0.0 

BOTTOM  TYPE  CITHCF  3.0  0 T 3.0*  3.0 

IMITIRL  TIME*  1100.0 

FIMRL  timc*’  0.0 

TIME  INCREMENTS*  1.0 

IS  THERE  R SURFACE  DUCT  YES 

SURFACE  SOUND  SPEED ■ FT  GEO-  4??G.O 

BOTTOM  SOW  IB  SPEED «TT  GEC»=  5 100.0 


DEPTH  TOP  BOTTOM  SOUND  GREED (FT • = 


1 3000 . 0 


Mill.  SOW  IS  SPEED  i FT  SEC)* 


4 


C*v  i • 


0 


DEPTH  RT  Mill.  SOUND  SPEED ■ FT • - .3300.0 

MR::.  SOUND  SPEED  ITT  SCO*  4??  7.0 

DcrTH  rt  i m.:.  sound  speed itti  = 1300.0 
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• Here,  "Initial'’  refers  to  1=0 

i 

t In  this  statement,  "initial" 
refers  to  the  starting  time  j 
of  the  encounter. 
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1PVIHC  A SURFACE  DUCT  lilTH  DOUBLE  BOUNCE  SOLUTION 


1170.  33992? . 3S9929.  0.929 
1 1 72.  S3? 929. 32992'?.  0. 930 


0.324E+01  0. 300E+01 
THEO*  0. 131E+02 
NODE  T TA  TIi  IF 
SDS  1100.  11 GO.  1179. 

0. 3344E+C»0 
0. 3333E+00 
0. 1134E-01 
0.0292'E-C'l 
0. 1 S77E-00 
0.  G003C«-00 
0. 02 1 7E+00 
0. 3074E-01 

0. 324 1 E+0 1 0 . 3000E+0 1 
0. 1997E+00 
2 BSP  1100.  II  GO. 

DP  1100.  117:3. 

0. 3544E+00 
0. 3G3GF >00 
0.  11S4L-01 
0.C292E-01 
0. 1S77E-00 
0. t OOG'E+OO 
0.521 7E+ 00 

DPP  1100.  !1~S.  1173.3:3? 

TO  CHANGE  A FUN  PARAMETER * 
flPPROPP I ATE  NUMBER: 

1 AIRCRAFT  1 :EL.  "tCTOF'S 

2 AIRCRAFT  INITIAL  FOSITION 
AIRCRAFT  HEIGHT 
RIRCPAFT  PRBIflTEB  FPEP. 

SUB  "EL.  "ECTOR 
SUB  INITIAL  FOSITION 
SUB  DEPTH 

SURFACE  SOUND  SPEED 
BOTTOM  SOUND  SPEED  AND  DEPTH 

10  MINIMUM  SOUND  SPEED  AND  DEPTH 

11  HIND  SPEED  AND  DIRECTION  COSINES 

12  BOTTOM  TYPE 

13  MR::.  SOUND  SFEED  AND  DEPTH 

14  TINE  OF  D 'ENTS 

15  PUN 
10  STOP 

ENTER  THE  APPROPRIATE  NUMBER  = 10 


OUTPUT  PARAMETERS 

R"  P DOP  D E PHI  P PATIO 

339929. 3S9929. 0.929  13.11  0.  4.220E-14- 


EXTRANEOUS : ONLY  FOR 
CHECKING  MECHANICS  OF 
PROGRAM 


TL 

•147. 


13.  1 1 

0 . 03 


u. 

0. 


1.212E-15- 
3. 330E-1 9- 


149. 

r.*3. 


1 •' '4  -''J'J  -"4  II  '4 

v . ■ u . • U • . 

ENTER  THE 


0.03 


0.  2. G9GE-15-14! 


o 

4 

5 

0 

< 

3 

Q 


EMIT 


• k'.’OS 

U l-r  1217?  11]  uuc  27  I SC  It  2 nrr  TTY  13  fit 
Cunnuct  t i pie  0:12  CPU’S  354 


10:24  AM  Thu  1S-Dcc-75 
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B . 3 Example  of  Program  Flow 

For  the  user  to  understand  the  program  flow,  an  actual  case 
will  be  simulated.  This  way  the  user  may  follow  through  to  see 
which  subroutines  and  formula  are  used  for  one  specific  case. 

As  shown  in  the  section  on  execution  of  the  program,  all  the 
inputs  are  manually  inserted  first.  Once  this  is  done  the  main 
program  calculates  a series  of  constants  set  up  by  the  main 
program  internally.  These  are  constants  which  do  not  change  with 
time.  A.  listing  of  these  constants  are  given  below: 


Constant 

GS  * ( CS-CMIN ) /DMIN 

GB  ■=  ( CB-CMIN ) / ( DB-DMIN ) 

Cl  * 1100 

AXN  - C1/C2 

C2  - 5000 

t « ( XIA-XIS ) ( VS-VAX ) - ( YI A ) ( VA  Y ) 

° [ (VS-VAX) 2+  VAY2] 


Computer  Prognm  Name 

GS  (surface  gradient) 

GB  (bottom  gradient) 

Cl  (sound  speed  in  air,  in  ft/sec) 

AXN  (ratio  of  air  sound  speed 
and  water  sound  speed) 

C2  (average  sound  speed  in  water 
in  ft/sec) 

TO  time  of  closest  approach 


R(tQ)  * C(Xs(t0)-Xa(t0))2+  Ya2(tQ)  + D2]1*  R0  slant  range  at 

closest  approach 


A 


r .iffK.p1  „ ] i 

*■  * 4 HlOO  4 F2^  J 


volumetric  atten  (DB/ft) 


PkHs  " P/100° 


B-ll 


Bolt  Beranek  and  Newman  Inc. 


TM  W307 


R’(t)  - [{x3(t)  - Xa(t))J  + 2 ( t ) ] ^ RP  Lateral  Range 

R(t)  * C{Xs(t)  - Xa(t))2  + Ya2(t)  + D2]^1  R slant  range 

If  a surface  duct  present  then: 

G1  * ( CMAX-CMIN ) / ( DMIN-DS ) G1 

GS  * (CMAX-CS)/DS  GS 


With  all  of  the  necessary  parameters  Inputted  and  all  the 
constants  Initialized,  the  main  program  sets  up  a common  storage 
area  so  that  all  this  data  may  be  shared  with  all  other  subroutines. 

Now  with  all  the  necessary  parameters  initialized  and  set 
up  so  that  they  can  be  shared  by  all  subroutines,  the  main  pro- 
gram now  tests  for  whether  there  is  a surface  duct  or  not. 

For  our  case  there  is  no  surface  duct  so  the  subroutine  BC'UN(THEft) 
Is  called.  This  subroutine  calculates  the  arrival  angle  of  the.  • 
ray  which  reaches  the  receiver.  At  a time  instant  t,  the  sub- 
routine first  tests  whether  there  is  a single  or  double  bounce 
and  this  is  done  by  testing  whetner  the  distance  to  the  conver- 
gence zone  Is  less  than  or  greater  than  the  lateral  distance 
between  the  receiver  and  the  virtual  source. 


Lateral  receiver/virtual  source  range  at  t 


b<b)=  Ve.)  ]'A 
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If  B(0)  > then  it  is  a double  bounce,  if  net  it  is  a 
single  bounce.  Now  an  iteration  routine  is  used  to  finally 
calculate  the  acquired  ray  angle.  Once  this  is  obtained*  and  for 
purposes  of  simulation  a single  bounce  solution  will  be  used  so 
that  subroutine  BS(THE9)  is  called.  The  following  calculations 
are  performed  within  BS(THE0). 


Boumce  elisramce:  B(eu')  = 


-sim®. 


0 


&.[  id'-c-Ff'*-* 


Travel.  Ti*e  a Long  8B  Ray  : ffB.)  * JL  f coah~'f  . — ^5—  ) - cosk"Y ! — . 

* 7 9s  L X.mW  ICOSG.JJ 


ARRIVAL  ANGLE:  * THEO 

ffl*)  * BteJ  tan  % +2cJ  J 

* 


•If  no  solution  for  either  single  or  double  bounce,  the  control 
is  sent  back  to  main  program  where  DP.F4  is  called. 
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AZIMUTHAL  ARRIVAL  ANGLE:  Cos  ) r { *ct  ( - Xs  lt  + t (e.))} 

-A)  - X,(t 4 tf*))}*  + yl  U - ±.)  ] y* 

AUXILIARY  ANGLE  FOR  DOPPLER: 

COSsp  = ****  + AT^  [-^Ct-h/e,)] 

1^1  [{*,(*+*♦.>)  )}*  + c/Vt-Ji)  +C2D  -D)2T/2 

C|  **  C|  ‘ 8 


DOPPLER  SHIFT: 

f.  _ 1 

ft  i m 

* p+^-ZIf^cos^J1 

- *.<*-*.)}*+  +avi>){]>i 

• 

TIME  OF  ARRIVAL  MINUS  TIME  OF  ARRIVAL  OF  DIRECT  PATH  FROM  CPA  POINT: 

- son  -Jl 

c»  cl 
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At  this  point,  the  angle  of  the  grazing  ray  to  the 
bottom  is  calculated  and  depending  on  which  bottom  type  (3  or 
5)  is  inputted,  subroutine  BT3.F4  or  BT5.F4  is  called  with  the 
grazing  ray  angle  and  returns  with  a bottom  loss*.  The  subroutine 
BT3.F4  or  BT5.F4  has  a series  of  linear  equations  which  are  good 
for  specific  frequency  and  grazing  angle  ranges.  The  subroutine 
interpolates  to  find  the  best  fit.  The  bottom  loss  values 
were  programmed  using  LRAPP  data. 


It  is  now  possible  by  combining  several  of  the  parameters 
calculated  already  to  obtain  the  mean  square  pressure  at  the 
receiver  and  the  TL. 

Jj®*  _ 8(nl  sim©Bcos8b  - + KBLfe4)J//0 

p*  - |0 

B(e»)  | asfe.) 

I de. 


TL=  io  Lo 3 J 


8 s in©ocos0, 


| aBfe0)| 


-p  f - ASf®#)  - RBLCefc) 


A correction  for  air  attentuation  is  added  to  the  Tl.* 


The  control  is  now.  sent  back  to  the  main  program  where  the 
values  of  MODE,  T,  TA,  TDIF , R»,  R,  DOP,  D/E,  PHI,  PRATIO,  TL 
are  printed  in  a preassigned  format. 


The  main  program  now  goes  into  a series  of  conditional 
statements  to  find  out  which  mode  was  Just  calculated  for  the 
smooth  sea  surface  case,  (which  In  this  case  the  mode  ■ BS)  and 
routes  the  program  to  call  the  subroutine  for  the  rough  surface 
case  which  for  BS  is  BSR(THE0 ) . 
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* See  Eq.  64,  Appendix  A. 
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Once  within  this  subroutine  a series  of  rough  surface 
formula  to  set  up  the  statistics  for  the  rough  surface  case 
are  calculated  as  given  below: 

Azimuthal  angle  between  plane  of  acoustic  path  ahd  wind  direction 
at  arrival  time  ta 

AUXILIARY  VARIABLE: 


Having  a value  of  zeta  (5),  the  subroutine  now  calls  a 
separate  subroucine  ERFS(  ) to  obtain  a value  of  the  error 
function  for  that  particular  value  of  £.  Once  the  value  is 
obtained  it  is  returned  to  the  calling  subroutine. 

Mean  square  sea  slope  (fully  arisen  sea) 

e1  = /..s*/ov^r  V Ci-er/a)] 

rms  slope  in  wind  direction 


tf„  = ovfe 
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Most  probable  slope  estimate 

V = U(szi)JL 

I ooo  sec  ' (go 

Mean  square  sea  surface  curvature: 


Auxiliary  variable: 

im  = hff 


Again  the  subroutine  erfs(  ) is  called  and  values  of  the 
error  function  are  obtained  for 


fin 


tan  9 


Once  these  values  are  obtained,  the  equation  for  the  average 
number  of  refracting  paths  (fully  arisen  sea)  can  be  calculated: 
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For  the  rough  surface  case  only  the  mode,  the  mean  square 
pressure  and  the  TL  are  modified  from  the  smooth  surface  case; 
all  other  output  Is  exactly  the  same.  The  new  modified  mean 
square  pressure  and  TL  are  given  in  the  equations  below: 


P* 

1m  = 8(n*c.s«. 

p*-  Brej  i aewi 

l-itf.  1 


{if1  * + 

VT  <rw  L 


_ (tan^+fljcos  fa,)* 


(taro  6. + c|>  cos  4>w)*]  -4-  £w  cos1  (tan  e.  + t|>  ces4>w)  e 2 

fzrT 


} 


- (fc Sf«0  + RBL(Gl) )/( o 

X N io 


With  the  two  new  values  the  control  Is  referred  back  to  the 
main  program  where  the  output  values  are  again  printed  of  MODE, 
T,  TA,  'TDIF,  R\  R,  DOP,  D/E, PHI,  PRATIO,  and  TL. 

The  main  control  program  now  automatically  goes  to  the 
subroutine  DP.F^  to  calculate  for  direct  path  propagation. 

Within  DP.F4  the  following  equations  are  solved. 

Depression  angle 


si  n 6ft)  = ±— 

[faw  - X,W}'+  tjlit)  +£>-JVl 
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Azimuthal  arrival  angle 

costyt)  = -Xcfr)} 

Lateral  range 


Slant  range 


Time  of  closest  approach 

to  _ ( Xu- X„ ) (w,  - <r)  - 

Time  difference  between  arrival  time  and  arrival  time  of  direct 
path  from  CPA  point 

At  * -t  -To  -f 

c* 

Time  of  arrival 

t * t + il_  + 

* c,  C* 

Auxiliary  variable: 

cos  W+J  = f Xs  W.  - X.rt)J  - %(t) 


<K  R<t) 
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Doppler: 

r 

22s>wsMrtO  j 

1-  zflr,  r /,(t) 

C*  J l 

c*  1 Rft)  J J 

Squared  Pressure : 

\ 

%. 

^mD 

Transmission  Loss: 

TL  * 10  -AR(t) 

The  control  is  now  sent  to  the  main  program  where  the  values 
outputted  are  printed  MODE,  T,  TA,  TDIr , R',  R,  DOP,  D/E,  PHI, 
PRATIO,  and  TL.  Now  the  final  main  loop  is  entered  which  is  the 
calling  of  the  subroutine  DPR.F^  which  does  the  calculations  for 
the  rough  surface  case  of  direct  path  propagation. 

The  initial  rough  surface  formula  for  the  statistics  of  the 
rough  sea  surface  are  the  same;  thst  is 

5 * , N 

are  the  same  exact  equations  only 

cos  $w(+ft) 

is  different  as  shown  below 

cos  4wf*A)  = 

I /*»«■>  - 4 
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Therefore  the  modified  mean  square  pressure  and  TL  for  th ; 
direct  path  rough  surface  case  are: 

p*  2 r 

-£r  = ) [=H 1 + gr * c°ste  fa* + frA"«  +&“*+«>*] 

-(tine  + $cos4>„f  — 

+ cosze  (-tine  -Hpcos^)  £ 

OiT  J 


TL  = 10  L«3,.J'^rj 

A correction  for  air  attenuation  is  then  added  to  the  TL.* 

The  final  values  are  returned  to  the  main  program  where  they 
are  printed  out  in  the  specified  format,  the  values  printed  at 
MODE,  T,  TA,  TDIF , R’,  R,  DOP,  D/E  PHI,  PRATIO,  TL. 

The  main  control  program  has  now  completed  all  of  the 
necessary  data  reduction;  it  then  asks  the  user  if  he  would  like 
to  change  anything  and  rerun  the  program;  if  not,  the  user  can 
exit  and  end  the  calculations. 


• See  Eq.  64  of  Appendix  A 
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B.4  Output  Description 

The  output  of  the  programs  differs  slightly  depending  upon 
(1)  whether  a surface  duct  is  present  or  not  and'  (2)  if  the 
program  is  outputting  single  or  double  bounce  results.  Examples 
of  the  output  for  each  of  the  four  possible  cases  are  given  and 
described. 
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CASE  t,  OUTPUT  (A):  NO  SURFACE  DUCT,  DOUBLE  BOUNCE  OUTPUT 


TRYING  ft  DOUBLE  BOUNCE  SOLUTION 


}^>.5£5E+01  1 L 

THE0= 

= 0.537E 

MDBE 

T 

BD 

-900. 

lOE+Oo! 


ff(2> 


CL596£E+  0 1 


0.  1000E+01  ® 


Tfl  TD IF 

k' 

R 

DDP 

D^E  [PHI 

>" 

RflT ID  1 TL  | 

83£.  -846. 

3£4 478. 

3£4479. 

1 . 076 

5.37(176.  | 

8. 

353E- 14-130.81 

! 

108. 159 

|100.  044 

-121.3J 

lpi.0  0 00E+  0 Ojg)1 

j 0.226  OE-~01  <g)| 

- £864E-££  © 


_0. 163 IE- 05  © 

drsoooE+bQ”®- 

*0.'  5£43E+  0 Jo .’  00  0 OE+'O  0 (&  | ® ® ($ 

(S>BHR  (2) -9 00 . fc.832  j# -J46 .J3£4 4 78 . 3£44 7971 . 076fa&“5737|f?6T 

®*[l  .08 . 1 59ll  0 0 . 044*3 

OP  -900.  -8£6.  -840. 3£4478. 3£4479. 1 . 081 

. 108. 159  99.319 

U_2.596£E+01@ 

0.  1 000E+01  (f 
0.  OOOOE+OO  (i 
0.  ££60F-01 
tL£864E-££  & 
e^iesiE-os  ©. 

£75000E+00,@ 

DPR  -900.  -8£6.  -840. 3£4478.3£4479. 1 . 081 

108.159  99.319 


± 


1.  1 16E-13hl£9.5j 

— 120. 

0.07176.  £. 63 IE- 18- 175. 8 
-166.4 


0.07176.  4. 4£4E-16-153.  5 
-144. 1 


For  each  particular  mode  of  propagation,  two  main  lines  of  output 
are  printed.  Numbers  1-14  are  the  main  output  values.  They  are: 


(1) 

Mode  of  propagation 

(2) 

Running  time  parameter  t 

(3) 

Time  of  arrival  (secs)  tA 

(4) 

Difference  in  time  between  arrival  time 
time  of  arrival  of  direct  path  from  CPA 

of  mode  t. 
(secs)  A 

and 

(5) 

Lateral  range  between  source  and  receiver  (ft)  at 
signal  leaves  source  (t) 

time 
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(6)  Slant  range  between  source  and  receiver  (ft)  at  time 
signal  leaves  source  (t) 

(7)  Doppler  shift  at  arrival  time  tA 

(8)  D/E  arrival  angle  at  receiver  (deg)  at  tA 

(9)  Azimuthal  arrival  angle  (deg)  at  tA 

(10)  Mean  square  pressure  (includes  bottom  loss,  water 
volumetric  attenuation,  no  air  losses) at  tA 

(11)  Transmission  loss  re  1 ft  (dB)  at  tA 

(12)  Lateral  range  (kiloyards),  at  time  signal  leaves 
source  (t, 

(13)  Lateral  range  at  time  of  arrival  (kiloyards)  at  tA 

(14)  Transmission  loss  re  1 yd  (dB),  includes  air  attenuation, 
at  tA 

Also  outputted  are  some  supporting  data  used  as  checks,  they 
are;  (15-27): 

(15)  Grazing  angle  at  bottom  (deg) 

(16)  Bottom  loss  from  LRAPP  data  for  this  grazing  angle 

(17)  D/E  angle  at  receiver 

(18)  £= 

(19)  ERF (^ ) 

(20)  (mean  square  sea  slope  (fully  arisen  sea)) 

(21)  Most  probable  slope  estimate  (radians) 

(22)  r f = 2-3XIO**  £ jL£frT  +£}f  _L 

v jtx*3  5 ‘ 

im  — k/r|P 
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(21!)  First  part  of  N (Average  number  of  refracting  paths) 


(25)  Grazing  angle  at  bottom  (deg)  (same  as  (15)) 

(26)  Bottom  loss  from  LRAPP  (same  as  (16)) 

(27)  Same  as  (2*1) 
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CASE  1,  OUTPUT  B:  NO  SURFACE  DUCT  — SINGLE  BOUNCE  OUTPUT 


THERE  IS  fl  SINGLE  EDUNCE 


0. 582E+00  0. OOOE+OO 
THE0=  0.128E+01 

MDUE  T TR  TDIF  R"  R DDP  D^E  PHI 

BS  -600.  -555.  -569. £17057. £17058. 1. 076  1.23174. 

72. 35£  66.967 

0. 5962E+01 
0. 1 000E+01 
0. 0000E+00 
0. 2260E-01 
-.  £864E-£2.  . 

0. 1631E-05 
0.  5000E+00 

0. 5821E+000. OOOOE+OO 


BSR 

-600. 

-555. 

-569.217057.217058. 1. 076 
'72.352  66.967 

1.28174. 

DP  -600. 
0. 5962E+01 

-547. 

-561. £17037.  217058. 1. 081 
72.352  66.102 

0.11174. 

0. 1000E+01 
0.  OOOOE+OO 

• 

0.2£60E-01 
-. 2864E-££ 
0. 1631E-05 
0. 5000E+0P 

• 

DPR 

-60  J. 

-547. 

-561.217057.217058.  1. 081 
72.352  66.102 

0.11174. 

P RRTIO  TL 
£. 258E-14-136. 5 
-127.0 


5.  983E-14-132. £ . 

— 1 2 c!  • 8 

1.34 0E-17-16S. 7 
-159.3 


1. 005E-15-150. 0 
-140.5 


For  Case  1,  Output  B,  the  format  Is  exactly  as  In  Output  A 
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CASE  1,  NO  SURFACE  DUCT 

INPUT  USED  IN  EXAMPLE 

CUTER  THE  APPROPRIATE  PARAMETERS  IN 
THE  DIMENSIONS  INDICATED. 

VELOCITY  VECTOR  X-DIRECTION <KTS> = ££0.0 

VELOCITY  VECTOR  Y-DIRECTION  <KTS> = 0.0 

INTIAL  X-COORDINATE  OP  AIRCRAFT < FT>  = 0.0 

INITIAL  Y-COORDINATE  OF  AIRCRAFT <FT> = £4000.0 

HEIGHT  OF  AIRCRAFT  FROM  SEA  SQRFACE  <FT) = 10000.0 

FREQ.  OF  AIRCRAFT  RADI ATION <HZ> = 150.0 

VELOCITY  VECTOR  X-DIRECTION  SUB<KTS>  = 7.0 

INITIAL  X-COORDINATE  OF  SUE <FT> = 0.0 

DEPTH  OF  SUB <FT)  = 400.0 

IS  THE  AIR  TEMP.  LESS  THAN  50F?  NO 

WIND  SPEED <KTS>  = 8.3 

X-DIRECTION  COSINE-  .707 

Y-DIRECTION  COSINE®  .707 

BOTTOM  TYPE  EITHER  3.0  OR  5.0=  3.0 

INITIAL  TIME®  -1000. 0 

FINAL  TIME=  1000.0 

TIME  INCREMENTS®  100.0 

IS  THERE  A SURFACE  DUCT  NO 

SURFACE  SOUND  SPEED <FT/SEC,  = 505£. 0 

BOTTOM  SOUND  SPEED <FT ✓SEC) = 5053.0 

DEPTH  FOR  EOTTOM  SOUND  SPEED<FT>=  15660.0 

MIN.  SOUND  SPEED  <FT ✓SEO  = 4875.  0 

DEPTH  AT  MIN.  SOUND  SPEED<FT>=  3440.0 
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CASE  2 OUTPUT  (A):  SURFACE  DUCT  PRESENT,  DOUBLE  BOUNCE  OUTPUT 


TRYING  A SURFACE  DUCT  WITH  DOUBLE  BOUNCE  SOLUTION 


( 0. 114E+02I 0, 116E+ 
(THE0t_P._144E+0£  j-'& 
NODE  T TA  TDIF 

BDS  640.  688. 


p/  r hop 

676.230421.230421. 0.931 
76.807  82.519 


D/E  PHI 
14.41  3. 


P RATIO 
4.11  0E-13- 
-116.2 


TL 

123.9 


0. 5962E+01  (?*J 


[pTlOOOE+Ol  ® 

0^000 0E+ 00  (ft 
0.2260E-01 
-. 2864E-22 
07l'63TE-05  (g 
0.5000EV00<S 
0.5418E-01  @ ^ 

0.  U45E+_02j0.  1 157E+0QI 
O^llSGE+OA^r 


? $ 


a 

t 


t 


I1-1,  x IStSpT.’-'P'**'  ! , r ■ - 1 * - • * 

(§}DSR  (26.40.  |3)688.  »676.:230421 .1230421  .|0.  931  0 14. 41 

L-  jgj)7”6 . 8 07F 8~2 .51 9 

634.^30421 . 230421 . 0. 935 


.T- 


& 


DP 


640.  695. 


0.10  3. 


76.807  83.411 


3.595E-i:-j 
QT  £i  16.8 
1 . 053E-17 
-162. 1 


-124. 4 | 

1 

-169.8 


0. 5962E+Q1  (/$> 


.0^1.0  0JE_±.01.^2>, 
Oi.OOOOE+00 
0.2260E-01 
-.2864E-22 


0. 163 IE- 05 


0. 5000E+00 


DPR 


640.  695. 


684.230421.230421. 0.935 
76.807  33.411 


0.10  3. 


0. 000E+00-999. 0 
-991.3 


For  each  particular  mode  of  propagation,  two  main  lines  of  output 
are  printed.  Numbers  1-14  are  the  main  output  values.  They  are: 

(1)  Mode  of  propagation 

(2)  Running  time  parameter  t 

(3)  Time  of  arrival  (secs)  tA 

(4)  Difference  In  time  between  arrival  time  of  mode  t.  and 
time  of  arrival  of  direct  path  from  CPA  (secs) 

(5)  Lateral  range  between  source  and  receiver  (ft)  at  time 
signal  leaves  source  (t) 
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(6)  Slant  range  between  source  and  receiver  (ft)  ct  time 
signal  leaves  source  (t) 

(7)  Doppler  shift  at  arrival  time  tA 

(8)  D/E  arrival  angle  at  receiver  (deg)  at  tA 

(9)  Azimuthal  arrival  angle  (deg)  at  tA 

(10)  Mean  square  pressure  (Includes  bottom  loss,  water 
volumetric  attenuation,  no  air  losses) at  tA 

(11)  Transmission  loss  re  1 ft  (dB)  at  tA 

(12)  Lateral  range  (kiloyards),  at  time  signal  leaves 
source  (t ) 

(13)  Lateral  range  at  time  of  arrival  (kiloyards)  at  tA 

(1*0  Transmission  loss  re  1 yd  (dB),  includes  air  attenuation, 
at  tA 

Also  outputted  are  some  supporting  data  used  as  checks,  they 
are;  (15-28): 

(15)  Grazing  angle  at  bottom  (deg) 

(16)  Bottom  loss  from  LRAPP  for  this  grazing  angle 

(17)  D/E  angle  at  receiver 

(18)  § )* 

(19)  ERFU) 

(20)  cTx  (mean  square  sea  slope  fully  arisen  sea)) 

(21)  Most  probable  slope  estlma  (radians) 

(22)  Jjf  * A.3x/o*  J £ A [/?{ l-erf(j)}  +£} [j-  - 1}  ] 
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(24)  First  part  of  N ("average  number  of  refracting  paths) 

W mO  cost*) } cos1 4, EffJ -f  (trto^,+  (pc6s^w)j  + qrw  cos'€jt*n9‘  +'Pces4#)* 

««f  [-  +<pe»s4vi)x/2°'Jl  J y 

(26)  Grazing  angle  to  bottom 

(27)  Bottom  loss  "rom  LRAPP 

(28)  Path  length 
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CASE  2,  OUTPUT  (B):  SURFACE  DUCT  WITH  SINGLE  BOUNCE  OUTPUT 


SURFACE  DUCT  WITH  SINGLE  BOUNCE  SOLUTION 

r — ■■  - Mi)*-, ©"V 

• I 0.21 73E+ Oq-.  151 3E+08  | _ _ 

13)0.034  ' 7<a>  0.  019  |(2J>  *0.014] 

GHU  0 0 . 0 00 ©499  0 . 4 0 0&4998. 5 0 Gj©4836 . 5 0 OQ,  5 05  0 . 00  0 | 
0.  422E+0  0 IQ . 0 0 0E+  0 0 i*-g» 


NODE 

BSS 


0.S32E+01  I 
T Tfi  TDIF 


R DOP 


620.  664. 


653.. 22324 0.223240.  0.  929 
74.413  79.695 


D/E  PHI 
8.82  3. 


P RATIO  TL 
1 . 713E-14-137,  7 
-130.0 


0i59_62E+01©  j 
|_0._1  OOOE+Ol  (g)j 
0TO000E+00<5S)i 
0. 226QE-01  Kb  1 
-.2864E-22<S)j 
™571 63  IE- 05 1 


0. 50QQE+0Q 
0. 1865E-01  _ ... 

0.42l6E+QClo.  OOOOE+OO^i 
0.2201E+06@|  Qfr 

fl)3SSR  foS20.i(B 

DP  620.  674. 


<2K/ 


664.10)  653 T £2324 0^22324 0 Jo. 9291(g)  3 .82j®  3 J” 1 .358E-14,'- 133/71 


0. 5962E+0 1 0 1 
O.JJOOOE+Ol  0| 
0.  0000E+00  ©J 
6. 226  0E- 01  © ' 
-. 2864E-22  & 
6. 1631E-05 i 


|®74.413[  79. 695Q>[ 
662.223240.223240. 0.935 
74.413  80.345 


<t$  -131.01 
0.10  3.  1.  19‘,E-17-169.2 

-161.6 


0.5000E+00 


DPR  620.  674. 


662.223240. 223240. 0.935 
74.413  80.845 


0.10  3. 


0.  000E+00-999. 0 
-991.3 


For  each  particular  mode  of  propagation,  two  main  lines  of  output 
are  printed.  Numbers  1-1 '•  are  the  main  output  values.  They  are: 

(1)  Mode  of  propagation 

(2)  Running  time  parameter  t 

(3)  Time  of  arrival  (secs)  tA 

(4)  Difference  in  time  between  arrival  time  of  mode  tA  and 
time  of  arrival  of  direct  path  from  CPA  (secs) 

(5)  Lateral  range  between  source  and  receiver  (ft)  at  time 
signal  leaves  source  (t) 
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(6)  Slant  range  between  source  and  receiver  (ft)  ai  time 
signal  leaves  source (t) 

(7)  Doppler  shift  at  arrival  time  t'A 

(8)  D/E  arrival  angle  at  receiver  (deg)  at  tA 

(9)  Azimuthal  arrival  angle  (deg)  at  tA 

(10)  Mean  square  pressure  (includes  bottom  loss,  water 
volumetric  attenuation,  no  air  losses) at  tA 

(11)  Transmission  loss  re  1 ft  (dB)  at  tA 

(12)  Lateral  range  (kiloyards),  at  time  signal  leaves 
source  (t ) 

(13)  Lateral  range  at  time  of  arrival  (kiloyards)  at  tA 

(14)  Transmission  loss  re  1 yd  (dB),  includes  air  attenuation, 
at  *A 

Also  outputted  are  some  supporting  data  used  as  checks,  they 
are;-  (15-38)j_ 

(15)  B(0q)  “ bounce  distance 

(16)  *B/*e0 

(17)  '01 

(18)  GS 

(19)  GB 

(20)  Cl  - air  sound  speed 

(21)  CS  - surface  sound  speed 

(22)  CMAX  - max.  sound  speed 

(23)  CMIN  - min.  sound  speed 

(24)  CB  - bottom  sound  speed 

(25)  Grazing  angle  to  bottom 

(26)  Bottom  loss  from  LRAPP 
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(27)  D/E  angle  at  receiver 

(28)  l 

(29)  erf  (5) 

(30)  Mean  square  sea  slope  (fully  arisen  sea) 

(31)  Most  probable  slope  estimate  (radians) 

(32)  ? 

(33) 

(34)  First  part  of  N (average  number  of  refracting  paths) 

(35)  + **  ecs1* +*CM*  7 

V5r 

(36)  Grazing  angle  to  bottom 

(37)  Bottom  loss  from  FACT 

(38)  Path  length  (SO.)) 
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CASE  2,  WITH  SURFACE  DUCT 

IN°UT  USED  IN  EXAMPLE 

..  LORD  WATER • REL » BO . REL  > SBO.  REL  > B T 3 . REL » BT5 . REL  > ERFS . REL 
LOADING 

* 

LOADER  UK  CORE 
13+3K  MAX  84  UORDS  FREE 

EXIT 

.ST 

ENTER  THE  APPROPRIATE  PARAMETERS  IN 
THE  DIMENSIONS  INDICATED. 

VELOCITY  VECTOR  X-DIRECTION <KTS) = £20.0 

VELOCITY  VECTOR  Y-DIRECTION  <KTS) = 0.0 

INTIAL  X-COORDINATE  OF  AIRCRAFT < FT)  = 0.0 

INITIAL  Y-COORDINATE  OF  AIRCRAFT <FT) * 12000. 0 

HEIGHT  OF  AIRCRAFT  FROM  SEA  SQRFACE  <FT>  * 10C00.  0 

• 0 

FREQ.  OF  AIRCRAFT  RADI AT  ION <HZ) = 150.0 

VELOCITY  VECTOR  X-DIRECTION  SUB<KTS>  = 7.0 

INITIAL  X-COORDINATE  OF  SUB <FT> = 0.0 

DEPTH  OF  SUE <FT) * 400.0 

IS  THE  AIR  TEMP.  LESS  THAN  5 OF?  YES 

WIND  SPEED <KTS)=  8.8 

X-DIRECTION  COSINE*  .707 

Y-DIRECTION  COSINE*  .707 

BOTTOM  TYPE  EITHER  3.0  OR  5.0=  3.0 

INITIAL ’TIME*  600.0 

FINAL  TIME*  700.0 
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Case  2 continued. 

TIME  INCREMENTS®  20,0 
IS  THERE  A SURFACE  DUCT  YES 
SURFACE  SDUND  SPEED (FT/SEC) = 4990.4 

BOTTDM  SDUND  SPEED <FT /SEC) = 5050.0 

DEPTH  FDR  BDTTDM  SDUND  SPEED<FT>®  15540.0 
MIN.  SDUND  SPEED  (FT'SEO  ■ 4836. 5 

DEPTH  AT  MIN.  SDUND  SPEED <FT) = 3720.0 

MAX.  SDUND  SPEED <FT /SEC) = 4998.5 

DEPTH  AT  MAX.  SDUND  SPEED <FT>*  420.0 
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B.5  Program  Listing 
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TYPE 


1 


6666 

6665 

6667 

6664 

2 

3 

4 

5 

30 
91 1 


WATER* \*n.F4 

COMMON  / I N/  VAX . VA Y * X I ft  * Y I A * H * F * VS  * X I S ? B * AXN  * BT  * U * WSX  * WSY 

COMMON  /OUT/  T*  R*  RP*  Tft*  TDIF*  TL* BOP  * BE  * PH I * PPSG * A > MODE 

COMMON  /CC/  C 1 * C£>  CS  > CB  * DB  * BS  > CM I N* CMAX  * DM I N*  ETA  * 6S * GB  * RO>  TO . P I 

1 G1 

M*0.  0 

TYPE  1 

FORMAT  <1H  , ''ENTER  THE  APPROPRIATE  PARAMETERS  IN'/ 
i'  THE  DIMENSIONS  INDICATED. ' * /) 

CALL  AIRV  <VAX» VAY) 

CALL  A IRC CXIft* VIA) 

CALL  AIRH <H> 

CALL  AIRF  <F) 

CALL  SUBV  CVS) 

CALL  SUBC  <XIS) 

CALL  SUED  CB) 

TYPE  6666 

FORMAT <1H  >'IS  THE  AIR  TEMP.  LESS  THAN  50F?' *£X> $> 

ACCEPT  6665* LANS 
FORMAT  <A3> 

IF  CLANS. EQ. 'YES')  GO  TO  6667 

AL=  C.  000000074^ <F**£) ♦H) / <50. 0*3. £81) 

GO  TO  3664 

AL=  <• 0 0 1 £5*F^H) / 1 0 0 0 . 0 
CONTINUE 

CALL  WIND  <U*  WSX*  WS‘  ) 

CALL  EOT  <ET) 

CALL  TIM <TJ * TF*  BT) 

TYPE  £ 

FORMAT <1H  , 'IS  THERE  A SURFACE  DUCT'*£X*$) 

ACCEPT  y*LMNS 
FORMAT  <A3) 

IF  CLANS. EQ. 'YES')  GO  TO  4 
CALL  ENVSCCS) 

CALL  ENVC  CCB*  BB) 

CALL  ENVB CCM I N * BM I N> 

GO  TO  5 
CALL  ENVS CCS) 

CALL  ENVC  CCB*  BB) 

CALL  ENVB CCM IN* BM IN) 

CALL  ENVft CCM AX*  BS> 

CONTINUE 
CONTINUE 
T=T  I 

CONTINUE 

GS= CCS-CM I N) /BM I N 
GB=  CCB-CM IN)/ CBB-BM I N> 

Cl*l 100.0 
AXN=. ££ 

C£=5000. 0 
PI *3. 14159265 
TX1=  CXIA-X1S) 

TX2= CVS-VAX) 

TX3=<YIA*VAY) 

TX4*CTX1«TX£)-TX3 
TX5=TX£**2 
TX6=TX5+ CVAY^S) 

TD=TX4/TX6 
RX1*  CVS^TO) +XIS 
RX£*CVAX*TO)+XIA 
RX3=  CRX1-RX2)  ♦♦£ 

RX4=  < <VAY>  TO)  +Y I A > ♦♦£ 

RO=SORT CRX3+RX4+ <&♦♦£> ) 
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333 


334 

31 

703 

705 


706 


4444 


4446 


THE 0=0. 0 

AX1=<F/1000.  0>~£ 

RX£=  <0. 1#AX1>/<1. 0+AXl> 

RX3=<40.  0*RX1>/<4100.  0+ftXl) 

R=  <AX£+AX3)  -'3000.  0 
RY1=<VS-*T)+XIS 
RY£=  <VRX*T>  -t-Xlfl 
RY3=  <RY  1 -PY£)  ♦♦£ 

RY4= < <VRY*T> +YI  A> 

RP=SOPT <RY3+PY4> 

R=SQRT (RY3+RY4+ (»♦♦£) > 

IF  <LRNS.  EQ.  "YES"!)  GO  TO  333 
CALL  BOUNCTHEO) 

GO  TO  334 

G 1 = CCMAX-CM IN)/ <DM I H-DS> 

GS= <CMRX-CS> /ns 
CALL  SBOUN  <THE0> 

TYPE  31 » THEO 

FORMAT <1H  ? 'THEG=' ? El  0.  3) 

TYPE  703 

FORMAT  <1H  ? "MODE  T TA  TDIF  R"  R DOP  D/E' 

1 ' PHI  P RATIO  TL'> 

TYPE  705? MODE?  T?  TA?  TDIF?  RP?  P? DOF  ? HE?  PHI ? PRSQ?  TL 

F0RMAT<1H  ? 1XA4? 3F6.  G?£F7.  0? F5. 3? F7. £? F4. 0? 1PE10.3? OPFb. 1> 

RPK=RP/3000. 0 

RY1= (VS^TA) +XTS 

RY£*  <VAX»TA) +XI A 

RY3=  <Rt  i-k’i  c>  ♦♦£ 

RY4=  < <VA  r>TA) +YI A) ♦♦£ 

RPA=SQRT <RY3+RY4> 

RPT=RPA/3000. 0 

TL1=TL+<10.  DIALOG  1 0 <9.  0>>-AL 

TYPE  706?  RPK? RPT  ? TL1 

FORMAT  < 1 H ? £3X ? 2F7 . 3 ? £ OX ? 0PF6 . I > 

I F <MODE . EG . ' BS ' ) GO  TO  4444 
IF  <MODE.  EC).  ' ED  ')  GO  TO  4446 
IF<MODE.EG. 'ESS')  GO  TO  4447 
IF  (MODE.  EC!. ' BDS'>  GO  TO  444S 
GO  TO  4445 
CALL  BSP (THE  0) 

TYPE  705? MODE?  T? TA?  TDIF?  RP?  R?  DOF/  DE? PHI ? PFSG?  TL 

RPK-RP/3000.0 

RY1=<VS*TA)+XIS 

RY£= <VAX*TA>  +XIA 

RY3=  <RY 1 -RY£> ♦ ♦£ 

RY4= < <VAY*TA) +YIA> ♦♦£ 

RPA-SQRT  CRY3+RY4) 

RPT =RPR/3  000. 0 

TL  1 =TL+  < 1 0 DIALOG  1 0 <9.  0>  ) -AL 

TYPE  7 06? RPK? RPT? TL1 

GO  TO  4445 

CALL  BDR<THE0) 

TYPE  7 05 ? MODE ? T ? T ' ? TD I F ? RP ? R?  DOP?  DE ? PH I ? PRSO ? TL 

RPK*RP/3000. 0 

RY1=<VS»TA)+XIS 

RY£=  <VAX*TA> +XI A 

RY3=  <RY 1 -PY£> ♦♦£ 

RY4=  < <VAY*T  A > +Y 1 A)  ♦♦£ 

RPR*SC!RT<RY3-tRY4> 

RPT=PPA/3000. 0 
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TL1=TL+  <1 0. DIALOG 1 0 <9.  0>>-AL 
TYPE  7 06  * RPK * RPT * TL 1 
GO  TO  4445 

4447  CALL  BSSR<THE0> 

TYPE  705* MOPE > T * TA*  TPIF  * RP*  R * POP > PE*  PHI > PRSQ , TL 

RPK=PP/3000.  0 

RY1=  <VS*TA> +XIS 

RY£*  <VAX*TA> +XIA 

RY3- <RY1 -RY£)  ♦ ♦£ 

RY4= < <VAY*TA> +YI A) ♦♦£ 

RPA*=SQRT  <RY3+RY4> 

RPT=RPA/3000.  0 
TL1=TL+<10.  O^ALOGl 0 (19.  0>)-AL 
TYPE  706* RPK * RPT » TL 1 
GO  TO  4445 

4448  CALL  BPSRCTHEO) 

TYPE  705>  MODE*  T * TA* TPIF* RP* R* POP*  PE*  PHI  * PRSQ*  TL 

RPK*RP/3000. 0 

RY 1 * <VS*  T A) +X 1 3 

RY£= <VAX*TA>  +XIA 

RY3= <RY1-PY£> ♦♦£ 

RY4*=  < <VAV*TA)  +YI  A)  ♦ ♦£ 

RPA=SQRT <RY3+RY4) 

RPT*RPA/3000. 0 
TL1=TL+  <10. O+ALOGl 0 <9. 0>)-AL 
TYPE  7 06 . RPK  * RPT  * TL 1 
4445  CALL  PPCTHEQ) 

TYPE  705*  rmriF . T.  TA.  TH I F.  RP,  R»  POP , PE  * PH  I * PRSQ  - TL 

RPK®RP/300G.  0 

RY1«<VS*TA>+XIS 

RY£*  <VAX*TA) +XIA 

RY3=  <RY1 -RY£) ♦♦£ 

RY4=  < <VAY»TA> +YIA> ♦♦£ 

RPA*SQRT <RY3+RY4) 

RPT«RPA/'3000.0 

TL1=TL+  <10. O^ALOGl 0 <3. 0>)-AL 

TYPE  7 06  * RPK  * RPT  * TL 1 

CALL  PPRCTHEO) 

TYPE  7 05  * MODE  * T * T A * TP I F * RP  * R * POP  * PE  * PH I * PRSQ * TL 

RPK=RP/3000. 0 

RY1=  <VS*TA> +XIS 

RY£=<VA»TA>+XIR 

RY3=  <RY 1 -RY£) ♦♦£ 

RY4=  < <VAY*TA»  + YIA> ♦♦£ 

RPA=SQRT <RY3+RY4) 

RPT =RP A/3  0 0 0 . 0 

TL1=TL+  <1 0. 0»ALOG1 0 <3.  0) > -AL 

TYPE  7 06 * RPK * RPT  > TL 1 

T*T+PT 

IF  <T. LT. TF>  GO  TO  911 
M=M+1 

IF <M. GT.  1>  GO  TO  30no 
TYPE  £000 

£000  FORMAT  <1H  *'TO  CHANGE  A RUN  PARAMETER*  ENTER  THE'/ 

1 ' APPROPRIATE  MUMPER:'/ 

1 ' 1 AIRCRAFT  VEL.  VECTORS' / 

1 ' £ AIRCRAFT  INITIAL  POSITION'/ 

1 ' 3 AIRCRAFT  HEIGHT'/ 

i 

K:- 
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2001 

3000 

3001 

3002 

47 

48 

49 
10 
11 
12 

13 

14 

15 

16 

50 
1? 

19 

20 
IS 

200 

201 


202 


1 ' 4 AIRCRAFT  RADIATED  FREQ.'/ 

{ ' 5 SUE  VEL.  VECTOR'^1 

1 ' 6 SUE  INITIAL  POSITION'/- 

1 ' 7 SUE  DEPTH'/- 

1 ' 8 SURFACE  SOUND  SPEED'/ 

1 ' 9 BOTTOM  SOUND  SPEED  AND  DEPTH'/ 

1 ' 10  MINIMUM  SOUND  SPEED  AND  DEPTH'/ 

1 ' 11  WIND  SPEED  AND  DIRECTION  COSINES'/ 

1 ' 12  BOTTOM  TYPE''' 

1 ' 13  MAX.  SOUND  SPEED  AND  DEPTH'/ 

1 ' 14  TIME  OF  EVENTS'/ 

1 ' 15  RUN'/ 

1 ' 16  STOP'/ 

1 ' ENTER  THE  APPROPRIATE  NUMBER  = M) 

ACCEPT  2001 >K3 
FORMAT <I> 

60  TO  <47s48s49s 10s 11? 12» 13* 14s 15s 16s 50s 17s 19s20s 30s 13>K3 
TYPE  3001 

FORMAT <1H  s 'ENTER  THE  HUMBER  OF  THE  PARAMETER'/ 

1 ' THAT  YOU. WISH  TO  CHANGE  = '»$> 

ACCEPT  3002s K4 
FORMAT <I> 

GO  TO  <47s43s49« lOslls 12s 13s 14s 15s 16»50s 17s 19s20s30s 18)K4 
CALL  A I R V (VAX  s VAY) 

GO  TO  3000 

CALL  A IRC  <XI  As  YIA> 

GO  TO  3000 
CALL  AIRH<H> 

GO  TO  3000 
CALL  AIRF  <F> 

60  TO  3000 
CALL  SUBV  <VS) 

GO  TO  3000 
CALL  SUEC  <XIS> 

GO  TO  3000 
CALL  SUED  <D> 

GO  TO  3000 
CALL  ENVS  <CS> 

GO  TO  3000 
CALL  ENVC<LL.l1B) 


GO  TO  3000 

CALL  ENVIKCMIN-DMIN) 

GO  TO  3000 

CALL  W I ND <U s WSX s WSY> 

GO  TO  3000 
CALL  EOT  <BT) 

GO  TO  3000 

CALL  ENVA<CMAXsDS> 

GO  TO  3000 

CALL  TIMCTIs  TFsDT) 

GO  TO  3000 
CALL  EXIT 
END 

SUBROUTINE  AIRV CVAXs VAY) 

TYPE  200 

FORMAT <1H  - VELOCITY  VECTOR  X-DIRECTION <KTS> =' - 2X< S> 
ACCEPT  201s  VAXKT 
FORMAT  <F) 

VAX® 1 . 689*VAXKT 


TYPE  202 

FORMAT <1H  s 'VELOCITY  VECTOR  Y-DIPECTION <KTS> =' s 2Xs *> 
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ACCEPT  201  » VAYKT 

VAY»1.683»VAYKT 

RETURN 

END 

SUBROUTINE  A IRC  <XIA>  YIA> 

TYPE  £00 

200  FORMATdH  .'INTIAL  X-COORDINATE  OF  AIRCRAFT  < FT^ =' > 2X. i> 
ACCEPT  201 > XI A 

201  FORMAT <F)‘ 

TY°E  202 

202  FORMATOH  .'INITIAL  Y-COORDINATE  OF  AIRCRAFT <FT> =' > 2X» $> 
ACCEPT  201 > VIA 

RETURN 

END 

SUBROUTINE  AIRH  <H> 

TYPE  200 

200  FORMATdH  .'HEIGHT  OF  AIRCRAFT  FROM  SEA  SORFACE  <FT)  =/  * c’X.  l) 
ACCEPT  201.H 

201  FORMAT  <F> 

RETURN 

END 

SUBROUTINE  AIRF  <F> 

TYPE  200 

200  FORMAT  OH  . ' FREQ.  OF  AIRCRAFT  RADIATION  CH2  > ='  . £'At$) 

ACCEPT  2 0 1 . F 

201  FORMAT <F) 

RETURN 

FNTl 

SUBROUTINE  SUBV<VS> 

type  n 00 

. 200  FORMA T < 1 H . "VELOCITY  VECTOR  X-BIRECTION  SUB <KTS> = ' . 2X» *> 

ACCEPT  201 » VSKT 
201  FORMAT  <F> 

VS*1.688»VSKT 

RETURN 

END 

SUBROUTINE  SUBC<XIS> 

TYPE  200 

200-  FORMAT dH  . "INITIAL  X-COORDINATE  OF  SUB <FT> = £X. *) 

ACCEPT  2(1 . XIS 

201  FORMAT  <F> 

RETURN 

END 

SUBROUTINE  SUBD<D) 

TYPE  200 

200  FORMAT < 1 H .'DEPTH  OF  SUB <FT)  = ' » £X»  *> 

ACCEPT  201. D 

201  FORMAT <F> 

RETURN 

END 

SUBROUTINE  ENVS(CS) 

TYPE  200 

200  FORMAT  <1H  . •'SURFACE  SOUND  SPEED <FT/'SEC>  =' > 2X. 

ACCEPT  201.CS 

201  FORMAT <F> 

RETURN 

END 

SUBROUTINE  ENVC <CB. DE> 

TYPE  200  : 

200  FORMAT  < 1 H ,' BOTTOM  SOUND  SPEED CFT^SEC) *' . £X. 1> 
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201 

202 


200 

201 

£02 


£00 

201 

£02 


£00 

£01 

202 

203 


200 

£01 


200 

£01 

£02 

. 203 


ACCEPT  £01 ? CB 
FORMAT <F> 

TYPE  £02 

FORMAT  <1H  ? 'DEPTH  FOR  BOTTOM  SOUND  SPEED <FT> = '? cX? t) 

ACCEPT  £01 » DI- 

RETURN 

END 

SUBROUTINE  ENVD <CMIN? DMIN> 

TYPE  £00  _ , / 

FORMAT  <1H  ? 'MIN.  SOUND  SPEED  CFT^SEO  = ?£X?$> 

ACCEPT  £01  ? CM IN 


FORMAT  <F> 
TYPE  £02 
FORMAT  <1 H 


? 'DEPTH  AT  MIN.  SOUND  SPEED <FT>  = ' ? £X?  $> 


ACCEPT  £01?  DM I N 


RETURN 

END 

SUBROUTINE 
TYPE  £00 
FORMAT <1H  ? 
ACCEPT  £01 » 
FORMA  T <F> 
TYPE  £0£ 
FORMAT <1H  ? 
ACCEPT  £01? 
RETURN 
END 

SUBROUTINE 
TYPE  £00 
FORMAT <1H  ? 


ENVA (CMAX?  DS> 

'MAX.  SOUND  SPEED <FT/SEC>='?£X?*> 
CMAX 


'DEPTH  AT  MAX.  SOUND  SPEED <FT> =' ? 
DS 

TIM  <TI ? TF?  DT> 

'INITIAL  TIME®' • £X?  $> 


£X?S> 


ACCEPT  £01 ? TI 
FORMAT  <F> 

TYPE  £02 

FORMAT <1H  ? 'FINAL  TIME='?£X?$> 
ACCEPT  £01 ? TF 
TYPE  £03 

FORMATCIH  ? 'TIME  INCREMENTS®'  ? £X? 

ACCEPT  £01?  DT 

RETURN 

END 

SUBROUT I ME  EOT  <BT> 


TYPE  £00 

FORMAT <1H  ?' BOTTOM  TYPE  EITHER  3.0  DK 
ACCEPT  £01 ? ET 


5.  0='  ? £X?  $:> 


FORMAT  <F> 
RETURN 


END 

SUBROUT I NE  W I MD  <U? MSX ? WSY  > 

TYPE  £00 

FORMAT  <1H  » 'WIND  SPEED 0(TS> = ' »£X? *> 
ACCEPT  £ 0 1 ? U 1 
FORMAT  <F> 

U«U1» 1.688 
TYPE  £02 

FORMAT  <1 H ? ' X-DI RECT I ON  CD  SINE®'?  £X ? 
ACCEPT  £01? MSX 
TYPE  £03 

FORMAT  <1H  ? 'Y-DIRECTIDH  CDSINE='?£X?t> 

ACCEPT  £01? MSY 

RETURN 

END 
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SUBROUTINE  DP(THE0> 

COMMON  ✓IN-'  VAX ? VAY ? X I A?  Y I A?  H?  F ? VS ? X I S ? B ? AXN  ?BT?U*  USX  ? WSY 

COMMON  /OUT/  T ? R ? RP»  TA?  TDIF?  TL ? POP ? BE  ? PH  I ? F'PSO . A ? MODE 

COMMON  /CC/  C 1 ? C2  ? CS  ? CB  ? DB  ? OS?  CM I N ? CMAX * DM I N?  ETA ? GS ? GB ? RO ? TO ? P I 

1 G1 

X1=<VS*T>+XIS 
X2=  <VAX*T>  +XIR 
X3=  <X1-X2> ♦♦£ 

X4=  <VflY*T>  +YIA 
X5«X4**2 

SIN=D/SGRT  (X3+X5+ <D**2> > 

DE1*ASIN(SIN> 

BE=  <DE 1*180.  0>/PI 
MODE* 'BP' 

C01=<VAX*T>+XIA 

C0£=X1 

C03=  (COl-Xl)  /'SORT  <X3+X5> 

C04=AC0S  <C03) 

PHI=(C04*1©0.  0>/PI 
TDIF=T-TO+  <R-RO>  /C2 
TA=T+  CH/Cl  > + <R/C£> 

ZX=  <VS*T> +XI S 

2X1  = <VflX*T> +XI A 

2X2=  <VAY*T> +YI A 

2X3=  (VAX* <2X-2X1 > > - <VAY*ZX£> 

ZX4=SQRT(VHX**2*VAY**£> 

ZX5=ZX3/  <ZX4*R> 

DOl = <2. 0*ZX4*ZX5> /C2 
BOZ-CQRTC.  O-DOl.' 

503=1.0/1102 
B05=  (ZX-ZX1J /R 
DOG*  <2.  0*VS*BD5'>  /C£ 

DO?=SORT  <1 . O-BOt..) 

DOP= <B03*B07 ) 

QZ1=-  <A*RJ/1  0.  0 
QZ2= <10. 0**CJZ1> 

PR1=4.  0* <AXN**£> ♦ <D**2> *022 
PRS0=PR1/<R**4) 

ETA=1 . 0/(10.  0**£0.> 

IF(PRSG.LT.ETA>  GO  TO  999 
GO  TO  993 
TL=-999. 0 
GO  TO  100 

TL=  (1  0.  0*ALOG  1 0 (F'RSGO  > 

RETURN 

END 

SUBROUTINE  DPR (THE 0> 

COMMON  / 1 H/  VAX ? VflY ? X I A. Y I A?  H?  F ? VS ? X IS?  B* HXN . BT * U ? USX , WSY 
COMMON  /OUT / T?  R?  RP?  TA?  TDIF?  TL?  BOP?  BE?  PHI . PPSO? A? MODE 
COMMON  /CC  ' C 1 ? C2 ? CS ? CE ? BE ? DS? CM I N * CMflX ? DM I N ? E T A ? 6 i?6B*  RQ«  TO?  F I 
1 G1 

MODE='BPR' 

COS 1 = <VS*T> +X I S 
COS£=(VAX*T>+XIA 
C0S3= (VAY* T ) ♦ V I A 

C0S4*  ( (WSX*  (C  0 : 1 -CQS£>  > - <WSY*C0S3>  :>  /PP 

2ET=C1/F 

2ET1  = <3.  0*H*ZET ♦ (32. 174**2))/,<<PI»*2>#(U**4» > 

ZET£=S0PT(ZET1>  , 

2ET3=SGRT  <ZET£> 
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TYPE  1111,  ZET3 
1111  FORMAT  <1H  , El  0. 4) 

ZET4-IU30.48 
. CALL  ERFS (ZET3, ERF) 

TYPE  111 If  ERF 

SIG*  <1 1500.  0*ZET4>  / <981. 5** 3) 

S1G1=SQRTCPI/£.  o> 

SIG£*  (SIG^SIGl) ♦ <1 . 0-ERF) 

TYPE  1111.3162 
SIGW=SORT  <SI6£>  x'3C!RT<£.  0) 

CY-  <£.  86*ZET4)  y'l  000.  0 
CY1-<CY*P1>/'180.  0 
TYPE  1111, CV 1 

ATI**  <£3000.  04  (SORT  <£.  0)  ) *4.  0)  ✓ <881 . 5*  <ZET4~3>  *3.  Q> 
AT£*  (SORT  <PI  > ) ♦ < 1 . 0-ERF) 

AT3=2ET3~£ 

AT4=EXP<-AT3)^ZET3 
AT5=<1. 0/ <£. ©♦ATS)  > -1 . 0 
ATN=AT1 * <AT£+  > AT4*AT5> ) 

TYPE  1111, ATM 
HX*H*30. 48 
EM=HX4SQRT  <ATM> 

TYPE  1111, EM 

EM1=1 . O'EM 

CALL  ERFS  (EMI,  ERF.) 

BN1*  < <1 . 0+ERF) ♦SORT  <PI)  ) ✓ <£. ©♦EM) 

BN£*<1. IV  <£!>!♦♦£)  ) 

BH3* < (EXP  <-BN£ ) +EM1 ) ♦ ♦£> ♦£. 0 
BN4*  < <EM««£)  ♦EH3>  / <4.  ((♦PI  > 

TYPE  1111, BM4 

SXN*D'R 

SXN 1 *AS I H (SXN) 

SXN£*<180.  0*SXM1.)/PI 
SXN3=S I NO  <SXM£> / COSO  <SXN£) 

SXM4*SXN3''  < (SORT  <£.  p>  .>  ♦SI GW) 

CALL  ERFS  <SXf  !4, ERF> 

SXN5=£RF 

SXN6*  <1 . 0-SIND  <SXM£>  > ^CtlSD  <SXM£> 

SXN7*SXH6/'  < (SORT  <£.  0»  > ♦SISW) 

CALL  ERFS <SXN7, ERF) 

BN5*SXN5+ERF 
BM5*EMS*PM4 
F01=SXN3+ <CY1#C0S4) 

F0E*F01  / < (SORT  <£ . 0 > > *31  GW) 

CALL  ERFS ‘"FOE.  ERF.) 

FO  3*  ( 1 . 0+E':’F  ) -£.  0 

F04=  < <S 1 6W**£)  + > FOl ♦♦£) ) ♦ (COSH <SXM£> ) ♦ ♦£ 
F05*(F04^F03' 

F06*  < <F01 ) ♦ ♦£»✓  <£.  0^'.i  I6W,,£)  > 

F07*EXP (-F06 ) 

F08* ( ( (COSD <SXM£)  ) ♦♦£> ♦F07»i IGW)  'S©PT<£.  ©♦PI) 

F09= (F05+F08) ♦BN6 
FOl 0=- <A*R) s\  0. 0 
FOll* (10. ©♦♦FOl 0) 

F0l£=F09^F0l 1 

pqi  * <£.  ©♦a:  :m  » -'R 
PRS©*<F01^E) ♦foie 
ETH= 1 . 0- <1 0. ©♦♦£© > 

IF (PRSO.LT. ETA*  60  TO  999 
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GO  TO  998 
999  TL=-999. 0 

GO  TD  100 

998  TL=1 0.  O^RLOGl 0 (PRSQ) 

100  RETURN 

END 

SUBROUTINE  DSP  (THE O') 

COMMON  / 1 N/  VRX * VRY  • X I ft  * Y I A * H*  F * VS * X I S * D « RXM  * E T * U * WSX  * WSY 
COMMON  /OUT''  T »R»  RP»  TR*  TD I F * TL  * DOR * BE  * PH I * PRSQ  * A * MODE 
COMMON  /CC/  C 1 . C£  * CS * CD * DE * DS  * CM 1 N * CMRK* DM IN* ETA* GS > GE  * RO  * TO  * P I * 
1 G1 

Bl  = ( (CMIN/CS) ♦ ♦£) ♦ (COSD (THEO) ) ♦ ♦£ 

Ec=  <,SORT  <1 . 0-El > .>  -S I ND  (THE 0> 

B3=  <,  (C8/CS)  ♦♦£)  ♦ (COSD  (THEO)  > ♦ ♦£ 

B4=S0PT  <1 . 0-B3) 

B5*  (SORT  <1 . 0-D1) ) -E4 
B6-(B£>GS)+(B5/GB) 

B7= (£. 0*CS*D6) /CO 3D (THEO) 

T 1 =CS/ (CM I N«CDSD (THE  0> ) 

T£*AL06  (Tl+SORT  ( (T 1 ♦♦£) -1 . 0.)  > 

T3*l . O/COSD (THEO) 

T4*AL0G (T3+SQRT < <T3**£> -1 . 0)  ) 

T5= ( (T£  T4> ♦£.  0)/6S 
T6*CS/  (CB^COSD  (-THE  0>  ) 

T7*RL0G (T6+S0RT ( (T6^£)  -1 . 0)  > 

T8=  ( (T7-T£) ♦ £.  0)/6E 

T9=T5-T8 

MODES'-  BSR' 

Pl=£7* (SI ND (THEO) ,'C OSD (THE 0)  ) 

P£=  ( (CMIN/CS)  ♦♦£>  ♦SIND  (THEO 
P3«  ( (CMIN/CS)  ♦«•£)  ♦ (COSD  '.THEO) ) ♦♦£ 

P4«SOPT ( 1 . 0-P3* 

P5* ( (P£/P4) - 1 . 0) /GS 
P6* (CB/CS) ♦♦£ 

P7*P6^ (COSD (THEO) ) ♦♦£ 

P3=SQRT (1 . 0-P7) 

P8=P6/P8 

PI  0= ( ( (CMIN/CS)  ♦♦£) /P4) -P9 

PI  1 = ( (S I ND (THE  0) > ♦P 1 0) /SB 

Pl£= ( (P5+P1 1)  ♦£.  0»CS) +P 1 

COS1* (VS*<T+T9)  ) +XIS 

COS£=  < VR»  (T-  (H/C 1 ) > ) +X I R 

C0S3* (VRY* (T- (H/C 1 > ) ) +Y I R 

C0S4*  (WSX*  (CO-S1  -C03£>  » - <WSY*C0S3)  ) /F7 

ZET*C1/F 

2ET1=(8. 0*H*ZET*(3£. 174**£>)/ ( (PI**£> ♦ <U**4) ) 

ZET£=SORT (2ET1) 

ZET3*S0RT (2ET£) 

TYPE  Till?  ZET3 
1111  FORMRT (1H  *E10.4> 

ZET4*U*30. 48 
CALL  ERFS (ZET3* EF:F) 

TYPE  1111*  ERF 

SIG=  d 1500.  0*ZET4> / (981 . 5**3) 

SIG1*SQRT  (PI  -'£.  0) 

SIG£=  (SIG*SIG1)  ♦ •'  1 . 0— ERF) 

TYPE  1111  * SIG£ 

SI  GW* SORT (SIG5> /SORT (£. 0> 

CY= (£. S6*ZET  4 > / 1 0 0 0 . 0 
CY1*(CY*PI)/180.  0 
TYPE  IUIjCYI 
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ATI*  <23000. 0* <SQRT <2. 0>>^4. 0) ✓ <981 . <ZET4^3> *3. 0> 

RT2*  <SQRT <PI)  > ♦ <1 . 0-ERF) 

AT3*ZET3~£ 

AT4*EXP <-fiT3)^2ET3 
AT5=<1. 0/ <2.  0*RT3>)-1. 0 
RTH=AT 1 ♦ <AT2+  <AT4^AT5> ) 

TYPE  1111 »RTN 
HX* <H^30. 4S> 

EM=HX^S0RT  <RTN  > 

TYPE  1 1 1 1 * EM 
EMI*  <1 . O^EM) 

CALL  ERFS  <EM1 » ERF> 

EMI*  < <1 . O+ERF.)  ♦SORT  <P I > > ✓ <2 . O^EM > 

BN2*  <1 . O-'  <EM^*2>  ) 

BN3* < <EXP  <-EM2> +BN1 ) ♦♦£> ♦£. 0 
BM4= < <EM^£> ♦BN3> / <4. 0*PI> 

TYPE  1 1 1 1 » BM4 

SXN3=SIND  <THE0> ✓COSD <THE0) 

SXM4*SXM3/  < <SQRT  <2.  0)  > ♦SIGW) 

CALL  ERFS  <SXN4«  ERF> 

SXN5=ERF 

SXN6=  <1 . 0-SIMD  <THE0>  > /'COSD<THEO> 

SXN7=SXNG/  < <SORT  <2 . 0>  > ♦$ I GW.) 

COLL  ERFS  <SXM7>  ERF> 

BN5*$XN5+ERF 
BN6=BN5*BN4 
F01-SXN3+ <CY1 ♦CQS4> 

F02=F0 1 ✓ < <SORT <2 . 0) > ♦$! GW) 

CALL  ERFS <FD2> ERF) 

FGX1.  CiERH.'C.  C 

F04*  < <SIGW**2>  + <F □!♦♦£> ) ♦ <CDSD <THE0) > ♦♦£ 

F05*  <FD4»FD3) 

FOG*  < <F01 ) ♦♦2)  ✓ <2.  0^  <SI6W~£>  ) 

F07*EXP <-F06) 

FDS*  < < <COSD  <THE0)  ) ♦♦£)  ♦FO?*SIGW)  /'SORT  <2.  O^PI) 

F09* <F05+F08) ♦ BNG 
TB=  <CB«COsD  <■*  HE  0)  ) ✓CS 
TB1=AC0S  <TB) 

TB2* (TB1M80. 0)-PI 
IF <BT. EO. 3. 0>  GO  TO  9993 
IF<ET ■ EO. 5. 0)  GO  TO  9999 
TYPE  999* 

9996  FORMAT <1H  « 'DIB  MOT  HAVE  A VALUE  FOR  BT*V> 

9993  CALL  BT3  <TB2»  RBL) 

60  TO  9997 

9999  CALL  £T5<TB2jRBL) 

9997  TYPE  9995 » TE2»  RBL 

9995  FORMAT <1H  >2E10.4) 

SI*  <2.  0*CS>  ••'COSO  < THEO) 

S2*  <1 . 0/GS>  + <1 . O/'GB) 

S3*ASIN<  <CMIM^COSD  <THE0> ) 'CS> 

S4=ASIN  < <CB*COII*  -TTHEO') ) /CS> 

S5=  < <PI '2.  0)  - < <THE0»PI>  ^ ISO.  O.O^'GS 
S6* < <-S£*S3)  + <S4/GE>  +S5) ♦SI 
FOl  0*-  <F:BL+  <A^S6)  >y\ 0.  0 
FOl 1*<1 0. ©♦♦FOl 0) 

F012=F09*F01 1 

P01*<8.  ©♦  < AXM»»2)  ♦COSH  <THE0)  > < tSIND  <THEO)  > ♦B7^ABS  <P12)  > 

PRS0=P01»F012 

TL=1 0. ©♦ALOGl 0 <PRSO> 

RETURN 
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SUBRQUT I HE  EDR  (THEO) 

COMMON  / I N/  VAX  * VAY  * X I A • Y I ft  » H*  F » VS  , X I S > D«  AXM * ET . U * USX  * WSY 
COMMON  /OUT/'  T * R • RP*  TA * T D I F » TL * HOP » DE  * PH I • PPS 0 • ft • NODE 
COMMON  /CC/  C 1 * C £ * CS  * Cl  * DB  * DS * CM 1 N > CMAX , DM 1 N * E TA . 6 3 * GB  * F:D*  TO • P I * 
1 G1 

Bl=  < <CMIN /CS> ♦ ♦£>  ♦ (COSD  (THEO)  ) ♦♦£ 

B£* (SORT <1 . 0-B1 ) > -S I HD (THE 0 ) 

B3=  < (CE/CS)  ♦ ♦£)  ♦ (COSD  (THEO)  > ♦ ♦£ 

84*S0RT(1. 0-B3) 

B5*  (SORT  <1 . 0-El  .»  ) -B4 
B6= (B£/GS) + (B5/GB> 

E7= (£. 0*CS*E6) /COSD (THE  0) 

T1 *CS/  <CMIM*COSD  (THE  0)  ) 

T£=ALOG  <T  1 +SQRT  < (T 1 ♦♦£)  - 1 . 0> ) 

T3*l . O/COSD (THEO) 

T4=AL06 (T3+SQRT ( (T3**£> -1 . 0)  > 

T5= ( (T£-T4> ♦£.  0) /GS 
T6=CS/  (CB*CDSD  (THE  0>  .> 

T7=AL0G (T6+SQRT ( CT***£) - 1 • 0) > 

T8=  ( (T7-T£) ♦£,  0)/GB 

T9=T5-T3 

MODE* 'EDR' 

P1=B7* (SI  ND  (THE  0> /COSD  (THEO)) 

P£= ( (CMIN/CS) ♦♦£> ♦SIND (THEO) 

P3= ( (CMIN/CS) ♦♦£) ♦ (COSD (THE  0) ) ♦♦£ 

P4=SQRT (1 . 0-P3>  ' 

P5= ( (P2/P4) -1.0) /GS 
P6* (CB/CS) ♦♦£ 

^Fov  (Cu jl* \ t ! ICC/  ) * * £ 

P3=SQRT (1 . 0-P7> 

P9*PG/P8 

PI 0=(( (CMIN/CS) ♦♦£> /P4)-P9 
PI  1*  ( (SIND (THE  0 ) ) ♦P 1 0) /6B 
P12*  ( (P5+F'l  1 ) ♦£.  0*CS)+P1 
COS 1 * (VS*  (T  + (£ . 0*  T9) ) ) +XI S 
COS£*  (VAX*  (T-  ( H/C 1 ) ) ) + X I ft 
CQS3* (VAY* (T- (H ' C D > ) +Y I A 

C0S4« ( (USX* (C0S1-C0S2) > - (WSY*C0S3) ) / (£. 0*B7> 

ZET=C1/F 

2ET1*  (8. 0*H*ZET* (32. 1 74 ♦♦£)>/( (PI ♦♦£) • <U**4  > ' 

ZET£*S0RT(ZFT1) 

2ET3=SQRT <ZET£) 

TYPE  1111*  ZET3 
1111  FORMAT (1H  • E 1 C . 4 ) 

ZET4=U*30. 43 
CALL  EPFS (ZET3* EPF) 

TYPE  1111* ERF 

SIG* (11500.  0*ZET4>/>981.5**3) 

SIGl*SQRT(PI/£.  0> 

SI  G£=  (S IG*  S I G1 ) ♦ ( 1 . 0-ERF) 

TYPE  1111*  SIG2 
SIGU=SORT (SIG£) /SORT (£. 0) 

CY*(£.8G*ZET4>/1400. 0 
CY1  = (CY*P I ) 180.  0 
TYPE  1111*  C V 1 

ATI*  (£3000.  0*  ( SORT  (£.  0>)*4.  0)  / *981 . 5* • ZET4**? • *3. 

AT£= (SORT (FI ) ) ♦ < 1 . O-EFF ) 

AT3*ZET3**£ 

AT4=EXP(-AT3> /ZET3 
AT5*(1. 0/(£. 0*AT3))-1 . 0 
ATN*ATl*(AT£*(ftT4*AT5'> 
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9996 
9998 

9997 
aacis 


TYPE  11 11» ATM 
HM-H»30.  48 
EM*HX«SORT<ATN> 

TYPE  111  If  EM 
EM1*(1. O'EM) 

CALL  EPFS <EMlf EPF) 

BN1* ( (1 . 0+ERF) *:ORT  (PI))x(£.  t>EM) 

BN2=(1. 0x< EM* ♦£  ' > 

BN3*  < (EXP  <- END  ♦ EMI >♦♦£>♦£. 0 
BN 4= < <EM**£) *BN3 • x <4 . 0*P 1 > 

TYPE  1 1 1 1 f EM4 

$XN3=SINB  <THE0>  xCOSD  (THEO) 

SXN4=SXN3x  < (SC'PT  <2.  0>  > *SIGW> 

CALL  ERFSCSXM4.EPF) 

SXN5*EPF 

SXN6*<1. 0-SIND(THE0>  > xCOSD <THE0> 

SXN7=SXN6x  < (SORT  <2.  Cu  ) ♦S16W> 

CALL  EPFS  <SXN7t  ERF > 

BN5* (SXN5+ERF) 

BN6«BN5*BN4 
FQ1*SXN3+  <CY1*C 0S4 • 

F02=F01  x < (SORT  (c  0 - • ♦SIGH) 

CALL  EPFS  (FDc'f  ERF  > 

F03*<1. 0+EPF  >x£.  (I 

F04*  < (SIGW**2>  ♦ > FOl  ♦♦£)  > ♦ (COST* ' THE  0>  .>  ♦♦£ 

FQ5*  <F04*F03> 

FQ6=  < (FOl) ♦♦£)  ✓ v2.  (•♦  <SI6W**2>  ) 

F07-EXP (-FOG) 

F03*  < < (COSO  •'  THE  C > ■ ♦♦£ ' ♦FD?»S  I GW)  xS  OF  T <2.  0*P ! ) 

F09* (F05+F03) ♦ ENt 
TB=  <CE*mSTi  (THEO  < > xCS 
TB1*AC0S(TB> 

TB2*  (TDl • ICO.  0)  'PI 
IF  <ET.  EO.  3.  O'  GO  TO  999? 

IF  <ET . EP. 5. 0 > GO  TO  ???? 

TYPE  9998 

FORMAT  (1H  f'OID  MOT  HAVE  A VALUE  FOP  £:T'  ) 

CALL  BTG  ' TESf  PEL  • . 

GO  TO  9997 
CALL  ET5 ' TE2«  PEL) 

TYPE  9995 > TE2» PEL 
FORMAT (1 H .£E10.4) 

SI*  <2.  (•♦CS.'xCOIB  '.THEO) 

S2*  <1 . 0 GD ♦ <1. O GE > 

S3=a:  i n •;  <cm  i n* r a : o • t he  o > * -■  c : > 

S4*Ai  I M ' C E*C 0 SO  • THE  0)  > . C i 1 

S5*  < (PI  x£.  <)>-<  ■ .'THE ('♦PI  • x 15  0.  O'  * xGS 

SG*  ( (-3  £♦£  3 ’ ♦ < 1 4 ■ Gt )♦!•?»  *3  1 

FO 1 0=-  •:  Kci . 0*6  EL  • ♦ •£.  •’•♦A*  1- G • '•  x 1 0.  0 

FOll* <10. 0**F01 O' 

F01£=F09*F01 1 

poi»<8. o* • m/n**2 ‘♦co:® <"the o • > x .4. o*s  1 me ••  t he 0 ' ♦ 1 7*ae : *pie*  • 

PPS0*P01«F01£ 

TL*1  0.  (••ALOGlO-PPiC' 

RETURN 

EMP 

suepoutime  e::p<theo> 

common  x in  ■/a:-:.vay.xia.via.h.f.v:.::i:.p.a::n.et.u.m:  .m:v 
COMMON  'OUT  x T.p.pp.  TA.  TP  IF . TL- POP  • I'E  . PHI . FRIO.  h.  r.OI'F 
COMMON  xCC  Cl  f Cc'»  CS*  CEf  I'E*  PI*  CMIN«  CMA  • PM1 N. ETA. GI . GT • FQ. TO- F I . 
1 Gl 

Bl*  <2.  (•♦CD  xCOIP • THEO' 

E2»  ( <C MA»C  Z > ♦♦£.*  ♦ • C 0 1 D • THE  0 ' > ♦♦£ 

B3-SOFT  a.  O-ED 
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B4»<SIND<THE0>-S3>.'GS 

B5«  < <CMIH/'C$>  ♦♦£)  ♦ tCOSD  <THEO>  > ♦♦£ 

BG«SDPT<1.0-B5> 

B?*a:6-F2)^G1 

B8*  < <CB/'Cr>**£>  ♦ <CDSIKTHE0> ) ♦♦£ 

B9«SGRT <1 . 0-B8> 

B1  0«  <.F&-E9<  'Gt 
Bll»<B4+£r*Bl(i'*Bl 

TH«flLDG  U1 . O'CO:  Ik  THE  O'  > ♦ iC’PT  < < < 1 . O^CDTIt  (THE  O'  ' ♦♦£>-!.  fr.  . 
TH1  *C:'^  <CMfl»CD  I Ii  ■:  THE  0 > > » 

TH£=ftLDG<  THl+'.OF  T < <TH1 ♦♦£> -1 . 0> > 

TH3* < <TH-TH£)  ♦£.  0>  G3 
TH4«CS^  <CMIN*C0n*  < THE  0*  ' 

TH5=fiLDG  |'TH4  + iOPT  < . TH4**£.-l . 0*  > 

TH6*  < <TH5-TH£ • ♦£. O'.'Gl 
TH?*CS^  <CF*CO:  rKTHEO-’  > 

TH8=HLDG  ‘,THr ♦'!  OF'T  • <.TH?  ••£ ' — 1 . O') 

TH9=  < (TH5-TH3  • ♦£.  to  /£% 

TH10=TH3*TH**TH9 

MDIiE*,BrSF' 

PI *E 1 1 ♦ < i I MB  1 THE  to  CO 'Ii  - THE  O'  ) 

P£*  < <CMPX^CS>  ♦♦£)♦  (SIHD  '.THE O'  > 

P3«  < <CMftX/C3  '♦♦£>♦ <C  0 3H • THE  O')  >♦♦£ 

P4*SGPT<1. 0-F?'  ' 

P5=  <1.0- <P£/P4 ’ ' Gi 
P6=  <CMIN.yCS"*  ♦♦£ 

P7»SC*RT  <1.0-  <P*»  < < C □ I D 1 THE 0>  >♦*£>> ' 

P8*PGP? 

P9=  <CMhX.'C:  ' *♦£ 
r: o-r9'ri 

PI  1 * < <P8-P  1 0 • ♦ l- 1 HD  •,  THE  o :•  :•  ✓ G 1 
• Pi£*\.:i.'C  :,<♦♦£ 

pi 3®  soft  <1 . o-  • pj £♦  < <co "■  r • the o*  *♦♦£''  > 

P14=PlcyPl  •: 

p 1 6=  < <p8-p i a ' ♦ i i nr<  < the  o • * gf 

PI  5* < <P5*P1 1 +P ! 6 ' ♦£. 0*0 : > +p  1 

con*  tvs* < t*thi  to  • ♦xir 

COS£=  (VAX*  • T-  - H - c 1 ' • ) ♦ s;ih 
CD  3 3 = tVR  Y*  < T - * H Cl')  > *Y  I h 

CDS4» t <w*» <co.  i -co j£)  •-  •i.i:y*-‘o.-*o  * fi  i 
zet«i;i/f 

ZE T 1 = < 8.  0*H*ZE T ♦ <?£ , 1 74**£ > >/(  <.PI  ♦♦£  * ♦ <i  l**4 .»  i 
ZET£=S0PT<IET1 ' 

ZE  T 3- SOP:  TCET£' 

TYPE  1111.  SET-  ? 

1111  FQPMRTtlH  . FI  0.4' 

2ET  4 =i.i*  ? ft.  49 
CRLL  EPFS • SETS. EPF > 

TYPE  1 111. EPF 

SIG*  <1 1500. 0*2ET4>  <931 . 5 • ♦ ;-o 

SIG1*:CPT‘PI  £.  to 

SI6£« <SIG*I IG1 •♦< 1 . 0-EPF  * 

TYPE  1 1 1 1 « ; IG£ 

S I GW-3CP T < * I G£»  •' : OP T ■ E . O' 

CY* t£. 3G*ZET4* - l 000. 0 
CY1*  <CY*PI  > 130.  ft 
TYPE  1111.CY1 

RT1«  <£3000.  0# • :opt < £ . 0 • •♦4.0)  <931 . 5* • ZET4**3 • *3. ft • 
RT£«<iOPT<PI  ' )*<  1. 0-EPF.* 

RT  3=ZET '*:•♦£ 

FlT4*EXP<-RT3>/ZET3 


# 
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RTS*  <1 . O''  <£.  0^RT3')’>-1 . 0 
RTM*ATn  (AT.;+  (RT4^RT3)  > 

TYPE  1111. ATM 
HK*H»30. 4S 
EM«HX*SOPT (RTM) 

TYPE  1111, EM 
EMI* <1. 0 'EM) 

CALL  ERF 3 < EMI, EPF) 

EN1*<<1.  (•♦EPF'1  ♦SORT  <PI  >>✓(2.  (•♦EM.) 
EN2*  <1 . 0 ' vEM**c  >:> 

BM3*  < (EXP  (-BH2)  +BN1  ) ♦ «£>  ♦,?.  (• 

£M4=  < (EM**S  • ♦£M3>  <4.  fi»PI  > 

TYPE  1111, EM4 

SXH3=  3 1 MD (THEO)  'CD _D • THE 0 > 
SXN4=rXM3-/  ( ( SORT  <£.  i.o  > ♦?:I6',I> 

CALL  EPFS  (SXM4. EPF, 

SXN5*EPF 


3344 


SXN8* (1 . 0-SIMP  • THEO)  > 'COSP  (THE  0> 

SXN7*SXM6-  ( (SORT  <£.  0 , > ♦SIGM> 

CALL  EPFS (SXN7, EPF) 

RM5*SXK5+EPF 

BM6=EN5*£M4 

F01*$XH3+(CY1^C0S4> 

FD2*F01/' < (SORT (£.  0>  >*SIGW> 

CALL  EPFS (FQ2, EPF) 

FD3*  <1 . 0+EPF)  /■£.  fi 

FD4*  < <i  IGW##2)  ♦ (FQl  ♦♦£•>  > « (COSP  (THE (•>  > ♦ ♦£ 

F05*  <FD4^F0  •:> 

FU6*  < CFO  1 > ♦ ♦£'-'  (•♦  • * I G'aI**£)  ) 

F07*EXP<-F09> 

Frw.r^Cf]-ruTu:ft,  ' * -£  > ♦ C’S7*  S J ■ 'SORT 'C.  . 

FG9*  <F0?.+F38>  ♦ EM8 
TYPE  3344, Fu? 

FORMAT C1M  , El  0.4, 


TE*  <CE*C OiP, THE 0 * ) 'Ct 
TE J *ACQS <TE>  • 

TE2*  (TE1 ♦ 1 80. 0 * 'P I 
IF  <ET.  EO.  3.  (•>  of?  TO  999$: 
IF  (ET.  EC1.  5.  O'*  GO  TO 
TYPE  999* 


??9i  FORMAT  ' 1H  .'HIP  MOV  HAVE  A VALUE  FOc'  ET  ✓ - 
393?  CALL  ET*- TEE* PEL) 

GO  TO  9r97 

399'?  CALL  ET?.  (TEE*  PEL  • 

3 39  7 TYPE  999? , T i E * P F L 

? 3?f'  eOF  MAT  ( 1 H * EE  1 0.  4 » 

Sl*(£.  0*Ci'  'CO'P'THEO', 

32*8003'  < 0 MA  ’*,;  C * P 1 THE  O ' * 'C  3 > 

S3* « •:  ».TH£0*PI » 180.  0 > -*£>  - GS 
i4*AC0S  • 'CMIM*CO:p,  THE 0*  > .'03 ) 

S5*<S4-i£»  61 

36*AC 0 3 < *. C I*C  0*  P • THE  0 > > C “ > 

S?*'S4-S8>  6E 
S8*(S3+S?.^I7*  *S1 
TYPE  8543*: 8 

854  3 FORMA  T ( 1 H , E 1 0 . 4 » 

FOl  (•*-  (PEL+‘A*S8  • ) 0.  M 

FOll  = • 1 0.  (•♦♦FOl  0> 

F01£*F09»F01 1 

PO 1 * <8 . •:•♦  • a:  :m  ♦ ♦£  ( ♦•:  0 : ri  • THE  0 * > / ( <3 1 Min  T HE  ft  > ) ♦ E 1 1 ♦ H E**  • P 1 M . . 

PRSOPOnF  01E 

TL*1  0.  (•♦Ac CGI  0 »FT'  SP‘ 


RETURN 

EMP 


B-5o 


Bolt  Beranek  and  Newman  Inc 


TM  W307 


SUBROUTINE  EDSR CTHEO) 

COMMON  / I N/  VAX > VRY » X I R?  V I ft ? H . F ? VS  ? X I S ? D > RXN » BT » U » WSX ? WSY 
COMMON  -'OUT  ''  T ? P?  F;P ? Tfi?  TDI F . TL ? POP ? DE  ? PH  I . PRSO * R ? MODE 
COMMON  /CC/  Cl » C£?  CS>  CE?  DB?  BS»  CMIN?  CMRXp  DMIN#  ETR?  GS»  GB?  PO«  TO.  PI 
1 61 

Bl=  C£.  0*CS)/CDSD  CTHEO) 

B2=  < <CMRX''CS>  ♦ ♦£)  ♦ CCOSD  (THEO)  > ♦♦£ 

B3=S0RT <1 . 0-B£> 

B4=  <SIND  CTHEO)  -E3)  /GS 

B5-  C <CM I N/'CS>  ♦♦£)  ♦ <COSD  <THE  0>  > ♦♦£ 

B6=SQRT  <1 . 0-1:5) 

B7= <B6-B3) /G 1 

B3=  < CCB/CS)  ♦♦£)  ♦ <COSD  (THE  0)  ) ♦♦£ 

B9»S0RT<1.O-BS)  • 

Bl 0= CBo-B9>/GB 
B1 Is  <B4+B7+B1 0) ♦Bl 

TH=RLOG  C C 1 . 0''" COS'D  CTHEO)  > +SQRT  < < C 1 . 0/CDSD  (THE  0.)  .)  ♦ ♦£.)  - 1 . 0)  ) 
TH1*CS/ CCMfl»CDSD CTHEO)  ) 

TH£=RLOG  CTH1+SQRT C CTH1 ♦♦£) -1 . 0) > 

TH3= < CTH-TH3) ♦£.  CO  /GS 
TH4=CS/  CCMIN*COSD  CTHEO)  > 

TH5=RL0G  CTH4+SQRT C CTH4**£) -1 . CO  > 

TH6= < CTH5-TH£) ♦ £. 0) /G1 
TH7=CSv  CCB*CDSD  CTHE  CO  .> 

TH3=RLDG  CTH7+SQRT C CTH?**£) - 1 . 0 > > 

7H9= C CTH5-TH8)  ♦ £.  CO  /'GB 
TH1 0=TH3+TH6+TH9 
MODE* 'ED SR' 

PI  =R1 1 • fs  ind  cthf  ro  /COSD  cthEO)  ;■ 

PC--  C CCMRX/CS)  ♦♦£)  ♦SIND  CTHEO) 

P3= < CCMRX/CS) ♦♦£) ♦ (COSD CTHEO) > ♦ ♦£ 

P4=SQRT<1. 0-P3) 

P5= Cl . 0- CP£/P4) .) /GS 
P6= CCMIN/CS) ♦♦£ 

P?=SQRT  a . 0-  CP**  < CCOSD  (.THEO)  > ♦♦£)  ) ) 

P3=P6/P7 

P9*  CCMRX/CS) ♦ • 2 

P10*P9/P4 

PI  1=  c (.P3-P1  0)  ♦SIND  (THEO)  ) - 61 
Pl£*  CCB/CS) ♦♦£ 

P13*S0RT  <1 . 0-  CP1£*  (.  CCOSD  CTHEO) ) ♦ ♦£.)  > ) 

P14*Pl£/Pl 3 

P 1 6 = < CP8-P 1 4 > ♦ ; I ND CTHEO > > /GB 
D15=  C CP*>+P1 1+P16)  ♦£.  0*CS)  +F1 
COS Is  CVS*  • T+  C£.  0 ♦ T rl  1 0 > > >+XIS 
COS£=  CVRX*  CT—  (.rl /C 1 .)  > ) +X I H 
CO  S •:=  ( VRY  • CT  - <H/£  1 ;• .)  + y I R 

COS 4=  < CMSX*  CCDS1-CDS2)  > - cMi'Y*CDS3»  > / (£.  i*i ♦ B:  1 1 j 

ZET=C1/F 

ZET1=  (3.  0»H*ZET ♦ C?.£.  174**£>  ) / ( (PI**£)  ♦ i :_l ♦ ♦ 4. :>  > 

ZET£=SORT  (.ZET 1 > 

ZET3=SC*RTCZET£.) 

TYPE  1111?  ZET3 
1111  FORMRT  C1H  ? E 1 0 . 4 > 

ZET4=U*30.43 
CRI.L  EPFS fZET3»  EPF > 

1 YPE  1111?  EPF 

SI6=  Cl  1500.  0*ZET4>  / (.981 . 5**3> 

SI 61 “SORT (PI/£. Co 
S ? 6£=  < S I G*  S I G 1 > ♦ (.  1 , Ci— EPF > 

TYPE  1 1 1 1 ? S I G£ 

SI6W=SC'PT (SI6£)  /SORT C£.  0) 
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CY*  <2.  86^ZET4>  "1  000.  0 
CY1*  CCY^PI)  180.  0 

TYPE  1 1 1 1 > CY1  ^ 

ATI* <£3000.  <SG!PT <£.  0))M.  0>/<9G1.5«'C2ET4^3>^3.  O' 
AT£=  CSQRT CPI) ) ♦ <1 . O-EPF) 

AT3*ZET3~£ 

AT4=EXP  <-AT3)-/ZET3 
RT5=  <1 . 0/  <£.  (•♦AT  3)  ) -1 . 0 
ATN=RT  1 ♦ CAT£+  < RT4^AT5.> » 

TYPE  1 1 1 1 * RTN 
HX*H^30. 48 
EM=HX*SQRT <RTN> 

TYPE  IIUjEM 
EMI* Cl. O/ EM) 

CALL  ERFS CEM1> ERF) 

EN1*  < Cl . O+ERF) ♦SORT  CPI) ) ✓ C£. O^EM) 

BN£»<1.0'CEM»*£>> 

£N3= < <EXP C-Bh...  EMI) ♦♦£>♦£.  0 
EN4=C  <£[■!♦♦£)  ♦EN3>/C4.  0*PD 
TYPE  1 1 1 1 > EM4 
SXN3=SIHD  CTHEO)  .'COSE  CTHEO) 

SXN4*S.XN3/'  C (.SORT  C£.  0)  ) ♦SIGW) 

CALL  ERFS  <SXM4j  ERF) 

SXH5-ERF 

SXN6*  Cl . O-SI MD  CTHEO ' > /'COSO  CTHEO) 

SXN7*SXM6-'  < (SORT  C£.  0>  ) ♦SIGW) 

CALL  ERFS  CSXH7 > ERF) 

BM5= CSXM5+ERF) 

BN6*EH3^£N4 

FL)l  = SXii:-r  CCY  i »CC34> 

FO£*FO  1 ✓ < CSQRT  C£ . CO  ) ♦ S I GW) 

CALL  ERFS  CFO£>  ERF) 

F03*  Cl . O+ERF)  •'■£.  0 

FD4*  < <S I GW*»£)  + CFO 1 ♦♦£> ) ♦ CCOSD CTHE 0) > ♦♦£ 

F05*  <F04^F03> 

F06=  <■  CFOl > ♦♦£)  / c£.  ('♦  cSISW»»£>  > 

F07=EXP C-F06) 

FD8*  C C CCOSD  CTHEO)  > ♦♦£)  ♦F07»SI6W)  ■''SORT  C£.  0*PI  > 

F09* CF05+F08) ♦ BK6 
TYPE  3331 > F09 
3331  FORMAT C1H  - El  0.4) 

TE= CCE^COSD CTHEO) ) 'CS 
TE1 =ACOS  CTE> 

TE£*CTE1»130.  CD 'PI 
I F CET . EC* . 3 . 0 > GO  1 0 999% 

IF  CET.EQ.  5.  0)  GO  TO  9999 
TYPE  9896 

8986  FORMAT  < lH  •'DID  MOT  HAVE  A VALUE  FOR  IT'-') 

9998  CALL  ET3 CTE£« PEL) 

GO  TO  9997 

9999  CALL  ET5 < TE£» PEL) 

9997  TYPE  9993 * T E£ • P EL 

9993  FORMAT  Cl H . £E 1 0. 4 > 

SI*  <£.  0»C i ) 'COSII  CTHEO) 

S£=ACOS  < <CMA»COS  P 1 THE  0>  > ' C S> 

S3*  < C CTHE 0*P I ) / 1 8 0 . 0.)  - S d>  / GS 
S4«ACO S •:  • C M I H*C  0 SB  1 THE  0 > > 'C S> 

S3*  CS4- S£  ) /G1 

S6*AC0 S < • CEK  OSD  < THE  0 * ) /CS> 

S7=  CS4-S6'  GE 
S3*  CS3+S5+S7) ♦£! 
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3332 


TYPE  3332*38 
FORMAT  <1H  * El  0.4) 
FOl  0®-  <.  <2.  (•♦PEL)  + <2. 
FOl 1*<10. O^rFOlO) 


0*A*S8>)^10. 0 


F012®F03»F01 1 

PQ1=  <3.  0*  <.AXM»*2.)  ♦COSD  (THEO)  > s 


<4.  O^SIND <THEp> ♦£! 1 ♦AP; 


PRSQ=PC!l*F0l2 

TL®  1 0 . 0*ALD61(KPRSC» 

RETURN 

END 


;•  ip  15>  > 
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.TYPE  B0.F4 

SUBROUT I ME  BOUN  (THE0> 

COMMOM  / 1 M " VAX . VAY  . X I A > Y I A . H . F . VS  * X 1 S » 0 > AXN . E T 
COMMON  /OUT/  T»  R»RP»  TA«  TDIF. TL» DOR. BE « PHI « PPSG. A.MOOE 
COMMON  /CC'  Cl  * C2.CS.CB.  DB.BS.CMIN. CMRX.  Ml  IN. ETA. GS«6B»  F'0«  TO*  PI » 
1 61 

X1»<VAX^  <T- <H/C1)) ) +X1 A 
X2=<VS*T)  +XIS 
X3=  <VAY»  <T-  <H  'C  l ) >)■*•¥  I A 
X4*  < <X1 -X£>  ♦ ♦£>  -*-X3^2 
X5=SGRT <X4 ) 

Rl*  <CMIN/'CB)  ♦♦£ 

R2=SQRT<1. 0-R1) 

R3=SQRT  <1.0-  <CS/'CB.>  ♦♦£> 

R4=  <R£-R3> /'GS 
R5=R£-'GB 

R6* <R4+R5) ♦£. 0»CB 
IF <X5. GE. RG)  GD  TD  100 
TYPE  1 

I FORMATCIH  > ' THERE  IS  A SINGLE  BOUNCE'/’) 

THE8=ACr)S<CS/Ct> 

THE* <180.  0»THE8i/'3. 14158265 
ZN*. 01 

II  DO  10  I = li.9000.1 
Z2=R_0AT<I) 

2*THt+  < <Z2- 1 . O')  / 1 00.  0)  ♦ . 0 1 
IF <Z. GT. 90. 0)  GO  TO  400 
TH*CS/ (CMIN^CUSD O > 

1 HI  =HL(JG  ( IH+  ;l.(K  i ( i H**2.-*  -i  . if*  t 
TH2=1. O/COSIKZ) 

TH3*AL0G  <TH2+ SORT  < <TH£**£)  -1 . O')  > 

TH4*CS/  <CB»COSB  <Z)  "> 

TH5*f  rJG  <TH4+S0RT  < <TH4»»£)  -1 . 0)  ) 

TH6*  (v  '♦  <TH  1 -TH3)  > 'GS 
TH7=  <2.  0*  <TH5- TH1) > /GB 
TH8*TH6-TH7 

XI  = <VAX* <T- <H/C 1 > > ) +X I A 
X2* <VS>  <T+TH8 > > +X I S 
X3*  <VA'r>  <T  - (H  ''C 1 > > > +Y 1 A 
X4=  < <X1  -X2.)  ♦♦£)  +XS**2 
X5=SC'RT  <X4 ) 

181=  < <CM!M /CS>  ♦ ♦£.'  * < 'COSO  <2> > ♦♦£) 

B2«S0RT»1.  n-f,l) 

B 3=  <£:£-. S I N 0 1 Z> /'G  S 

E4=  < <CI:/’CS)  ♦♦£.)  ♦ < <C OS  0 '•  Z>  > ♦♦£) 

B5*SPRT <1 . 0-E“ • 

Fb*  (BE— 1:5)  ' GF 

FT*  < <B3*Fb 1 »2 . 0*CS)  / (.COSO  < 2)  > 

F3=<E:?-X5>/E:7 

I:9*AI:S(E:8> 

IF (B9. LE. ZN)  GO  TO  41 
10  CONTINUE 

2N=ZH+. 01 
G*.  1 

IF  <ZN.LE. 0)  GO  TO  11 
400  TYPE  40 

40  FORMAT (1H  .'HO  SOLUTION  TO  SINGLE  BOUNCE'-) 

GO  TO  1U0 

41  THE0=Z 

CALL  FSdHEO) 

GO  TO  42 

B-54 


Bolt  Beranek  and  Newman  Inc 


TM  W307 


100  TYPE  43 

43  FORMAT <1H  * " TRYING  A DOUBLE  BOUNCE  SOLUTION' /) 

THE9=AC0S<CS./CB> 

THE1*  <180. 04THE9>/3. 14159365 
2N 1 = . 01 

13  !«□  99  1 = 1 > 9000*  1 

23=FL0AT  <I> 

Zl«THEl+<<23-l.  O.'/IOO.  0)4.  01 
IF  <Z1 . GT. 90. 0)  60  TO  988 
THB=CS-/  <CM1M*CDSD  <Z1)  > 

THB1 =ALOG  <THB+SGRT < <THE4*£> -1.0) ) 

THB£=i. 0/C03D<Z1) 

THB3=RL06  <THE:£+SCRT  < <THE244£>  -1.0)) 

THB4=CS/  <CB*COSD <Z1) ) 

THB5=AL0G <THE4+SORT  < <THfe4*»£) -1 . 0)  ) 

THE6=  <8.  04  '.THE  1 - THR3> ) -'GS 

THP?=  <£. 04  <THE5-THB1))/6B 

THB8*TH£6-THB7 

XB1 = <VAX4 <T- <H/C 1 > ) > +X I A 

XB£= <VS>  <T4  <£ . 04 THE8) ) ) +X I S 

XB3=  <VAY4  <T-  <H''C  1)  > ) + YIA 

XB4= < <XB1-XE£)44£) +XB344£ 

XB5=SQRT <XE4> 

BB 1 = < <CMIN/CS)  44£>  4 <COSD  <Z1)  ) ♦♦£) 

BB£=SGRT<1. 0-BE1) 

BE3=  <BB£-S I MB <2 1 ) ) / GS 

BB4*  < <CP./CS)  ♦♦£)  ♦ < <COSB  <21 ) > ♦♦£> 

BB5=SGRT <1. 0-EE4) 

BP6- ''PBe'-EE®1'1  /G* 

BB7» < <BB3+BF6 >44. 04CS> /* <COSIi <2 1 > > 

BE8=  <EE7-XE5)  -"EE7 
BE9=AES  <BE8> 

IF  <BB9.  LE.  2ND  GO  TO  97 
99  CONTINUE 

2N1-2N1+.01 
01*.  1 

IF <ZH1 . LE. 01)  GO  TO  1£ 

988  TYPE  98 

98  FORMAT  <1H  * ' NO  SOLUTION  TO  POUELE  BOUNCE*'/') 

GO  TO  4£ 

97  THE 0=21 

CALL  BIKTHEO) 

45  RETURN 

END 

SUBROUTINE  BS <THE0) 

COMMON  /'IN/  VAX < VAY»  XI  A-  YIA*  H*  F • VS*  XIS*  D * A XN * B T 
COMMON  'OUT./  T*R*RR*  TAjTPIFjTL*  PUP  • Be  * F*H  I * F F SO  * A « MOPE 
COMMON  /CC/  Cl  * C£* CS*CB, BB* BS* CM  IN. CMAX* BM1N*  E TA«6S * GB«  *0, TD*  F I * 
1 61 

Bl  = < <CHIN/CS) ♦♦£> 4 < CO SB <THE 0> > ♦♦£ 

B£= <SORT <1 . 0-El > ' -SIND r THEO • 

B3=  < <CB/CS) 44£»  4 (CO IB  - THE 0) > ♦♦£ 

£4= SORT '1.  0-B3 > 

B5= ». SORT  < 1.0-El)  > -E4 
B6=  <B£ /GS ) 4 • E5GE > 

B7* <£. 04CS4E6* /COST"  THEO) 

T 1 =CS / <CM 1 N4  c 0 I T"  THE  0 • > 

T£=AlOG  <T1*I0RT  < '.T  l44£>  -1 . 0>  > 

T3*l.  0 /CO II": THEO) 
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T4aRL0G CT3+SQRT  < 0>  > 

T5* C <T£-T4> ♦ £. 0> 'GS 
T6*CS',CCB*COLBcTHEO>  > 

T7=RL0G  CTG+SOPT  < 0)  > 

T8* C <T7-T£> ♦£. 0>'6B 

T9=T5-TS 

DE*THE0 

Tfl=T+T9 

P1=B7*  CSIND  CTHEOJ  /COSB  CTHE 0)  > 

P£=  C CCMIfVCS)  ♦♦£)  ♦SI  HD  '.THE  0> 

P3=  < <CM 1M"CS> •♦£  > ♦ cCDSD 1 THE  0> > ♦ ♦£ 

P4=S0RT  <1 . 0-P3> 

P5=  < <P£/P4> -1 . Q>  /'GS 
P6= <CB/CS) ♦♦£ 

P7*P6* CCDSD <THE 0>  > ♦♦£ 

P8=SQRT  <1 . 0-P7) 

P9*=P6/'P3 

PI  0=  < C COM I H'CS>  ♦♦£>  ■/P4.»  -P9 
PI 1 = C CS I HD CTHE 0 > ) ♦P 1 0 > 'GB 
Pl£= C lPS+PI 1> ♦£. 0*CS> +P1 
X 1 = <V  AX*  < T - CH'C  1 > > > +XIA 
X£= CVS* <T+T9> > +X I S 
X3*  CX1-X£> ♦♦£ 

X4= CVAY*  C CT- CH/C 1 > ' » > + V I fi 
X5* CX1-X£> /SORT  CX3+ <X4**£> > 

X6=RC0S  CX5> 

phi  = ciso.  ■ R i 

CD1  = CVS*  CT*T9> > +XIS 

C02» CVRX*  CT-  cH/C  1 > > > *>'1H 

C03*VRX*  CCD  1 -•:□£> 

C04=- < CVAY*  CT-  <H/C1 ) ) > + YIA* 

CD5*C03+ CVAY*C04 > 

CD6=SQRT  <VRX**£*VAY**£> 

CO?=RBS  CCOG.) 

C08«  CCD  1 -•:□£>♦*£ 

CD9=CDS+  CCD4**£ ) + C c£.  0*PB>  -D.'  ♦*£ 

C01  0=C07*SG'RT  <C09> 

C01 1=C03'CD1  0 

DO  1 = CC  06**£  '• / 1 C £♦♦£  > 

DQ£= C C£. 0*C 07> *C G 1 1 Li 
DD3= 1 . 0+D01-DQ2 
DD4*SQ6'TCDD3> 

DDR  1 * 1. 0 • I'u4 
DD5=CVS/C£.)  ♦♦£ 

DDG=  C C.SClRT  cCQ3>  > ♦£.  0*VS>  - C£ 

DD7=1 . u + P'13-  • DC-:  • SORT  ' CD?  • > 

DlF'=  c'DOP 1 > ♦ *C'R1  CDD7> 

MODE='BS 

T D I F = T + T 9 - T □ - R □ / C £ ' - • H ' C 1 > 

TB=  CCB*C D : D CTHE  0 .>  > -CS 
TBl*RCDScTB> 

TB£=CTB1*180.  0-  'PI 
IFCBT.EC>.  3. 0>  GO  TD  =>998 
IF  CBT.ECU  5. 0>  GO  TD  9999 
TYPE  9996 

9^96  FORMAT  < 1H  .'DIP  HOT  HAVE  A VALUE  FDR  FT  • 

9998  CALL  BT3 CT  B£ . RBL > 

GD  TD  9997 
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9999  CALL  BT5  <TB£»  RBL> 

999?  TYPE  9995 » TB2*  PPL 

9?95  FORMAT  <1H  >£E10.3> 

SI* <2. 0*CS) xCOSD CTHEO) 

S2*  < 1 . OxgS)  + C 1 . O' GB  > 

S3*  AS  I N < (CM I N*COS  D (THE  0> ) xCS) 

S4=ASIN  t (CB*COSD  (.THE to .)  /CS) 

S5*  < (PIx£. 0) - < (THE0*PI>  xlgfl. 0>>xGS 
S6=  < C-S£*S3>  + (S4x6B>  +S5)  ♦SI 
Q21*-((R*SG>+RBL>xl0.  0 
Q22*  <! 0. 0**Q21> 

PQ1=3. 0* CAXN**2> ♦SIMD  CTHEO) ♦COSB CTHEO ♦QZS 
PRSQ=PQ 1 x <£?*A£S tPl£>  > 

TL*  <10. Q*AL061 0 (PRSQ) ) 

RETURM 

END 

SUBROUTINE  BD(THEO) 

COMMON  x 1 NX  VAX j VAY» XI  A*  YIA?  H*  F • VS»  XI S*  IN AXN • BT 

COMMON  xD JTx  T*  Rj  RP>  TA«  TDIFj  TL  j HOP , I»E * FH I • PRSO > A > MODE 

COMMON  xCCx  C 1 » C£>  CS>  CBj  DB» DS  * CM I N * CMAX . DM I N» E T A . GS » G£» RO , TO • P I 

1 G1 

Bl*  < <CM I NxCS> ♦♦£) ♦ CCOSD  CTHEO.) > ♦ ♦£ 

B2=  <SQRT  < 1 . O-Bl > .)  -S I ND  CTHE 0.) 

B3* ( (CBxCS) ♦♦£> ♦ CCOSD CTHEO) ) ♦♦£ 

B4=SQF:T  Cl . 0-B3.) 

B5*  (SORT  Cl . 0-B1.)  > -B4 
B6«CB£xGS>  + CB5xGB.) 

B?*  C2.  0*CS*B6>  xCOSD CTHEO 
Tl*CSx  (CM1N*C05D  CTHEO) ) 

T2-AL0G  CT1+SQRT C CT1 ♦♦£> -1 . 0) ) 

T3*l.  OxCOSD CTHEO 

T4*AL0G  (T3+SGRT  C <T3**£>  -1 . 0.)  > 

T5*  < CT2-T4) ♦£. 0>  xGS 
T6*CSx  (CB+CQSD  CTHEO) ) 

T?*ALOG  CT6+SQRT  C CT6**2>  -1. 0)  ) 

T8*  C CT7-T2.)  ♦£.  .DxGB 

T9=T5-T8 

DE*THE0 

TA*T+ C£. 0*T9> 

TDIF=T+  <2.  0*T9> -TO-  CHxCl)-CRQxC£) 

F 1«B7*  (SIND  CTHEO  -'COSD  c THE  0>  > 

p£*  c ccm i n xcs ' ♦♦? ) ♦sinn c the o > 

P3*  C CCMINXCS) ♦♦£> ♦ CCOSD CTHEO) > ♦♦£ 

. P4*S9RT<1.0-R3> 

P5*  C CP2/P4 ) - 1 . 0) xGS 
P6=  cCBxCS) ♦♦£ 

P7*P6*  CCOSD  CTHEO  ) ♦♦£ 

P5*S0RT  Cl . 0-P7) 

P9*PGxP8 

PI  0=  C C CCMINXCS)  ♦♦£.)  xP4 ) -P9 
PI  1*  C (SIND  CTHEO) ♦PI  0> ' GB 
P1£*C  CP5+P1 1) ♦£. 0*CS> +P1 
X 1 * (VAX*  C T- • H C 1 ) > • +X I A 
X£=  CVS*  CT+  ■ 2.  (•♦T9.) ) ) +XIS 
X3* CX1-X2) ♦♦£ 

X4= CVAY* C <T- cH'Cl ) ) ) > +YIA 
X5= CX1-X2) xSCRT CX3+ CX4**£) ) 

X6*AC0S  CX5> 

PHI* (180. 0*X6>xPI 
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5>?5»€. 

9999 

9997 

4^35 


CDl  = <VS*<T+<£.  0»T9»>)+XIi: 

C02=  <VA»  <T- -H>C  1 ) > .)  +X  I A 
CD3*VftX*<CDl-CD2> 

C04=-  ( <VftV* <T- <H- Cl> > ) +YIR) 

CD5-CD3+ <VfiY*CD4> 

C06=SC‘RT  < V A ♦ ♦ 2 + V H Y ♦ ♦ 2 ) 

CD7*RBS  •'  •:□*£> 

C03=  CCD  1 -CD£>  ♦♦£ 

C09=CD?+ <C04»*2> ♦ < <4. 0*DB) -B> ♦♦£ 

COl 0=C07*S0RT  'X09> 

C011*CD5'C-010 
POl*  <C.O<S~£>  <C2**2> 

P02=  <.  < 2.  0*CD7>  ♦CDl  1>X£ 

D03*l. 0+D01-D02 
D04=SQRT tD03> 

DOP1  *■  1 . 0 DO4 
P05=  <VS'C3> ♦♦£ 

1)06= <.  <CQRT  <C08> > ♦£. 0*VS> /C£ 

H07=t . 0+H05-  ‘.IiOt.  'SQRT  <C09>  > 

DOP=  <DOPl)  ♦SORT  • • H07> 

TB=  <CB*COSli  <THE0>  > - C£ 

tbi=acoz<tb> 

TB2=  CTBl^lCO. 0»/PI 
IF <BT. EQ. 3. 0>  SO  TO  9993 
IF  <BT. EQ. 5. 0>  SO  TO  9999 
TYPE  9996 

FORMAT  <1H  > ' Bill  HOT  HAVE  A VALUE  FDR  BT'x> 
CALL  BT3 <TB£>  RBL) 

so  to  7 

CALL  BT5(.TB2»RBL.) 

TYPE  9995 *1 £2? PEL 
format  an  »£Eio.3> 

S 1 = <2 . 0*C O COSD  < THE 0) 

S2=  <1 . 0-'GS»  ♦ a . 0-'SB> 

S3-ASIM  < COM  I H* CO  S B <THE  O'* ) -"CO 

S4=ASIM < (CE*COSD «.TH! Co  > CO 

S5=  C < P I /£ . 0>  ~ ( *;TH£Ci«F'I  > "ISO.  CO  > --Si 

S6=  < C-J2*S3>  ♦ <S4-"GFo  + i'5>  ♦ $:  1 

ClZZ  1 = - < '.2.  Ci*h*SP.‘  + •: 2 . u«F  £'L  > > / 1 0 . 0 

0ZZ2=  1 0 . u ♦ »QZZ1  ,* 

M0fE=  t D ' 

PQ1  =S.  ()♦  C A XN  «♦:*'*  ♦SIMP  i THE  Co  ♦COSD  '.THE  0.:*  ♦ OZ  Z2 
PRSO*P Cl  1 \ t 7 * a B i *"P  1 2 '■  .*  *4 . O.j 

TL=  '1  Ci.  0*AL0S  1 0 CC'P  VC'  > '• 

RETURN 

EMD 
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.TYPE  SB0.F4 

SUBROUTINE  SBOUN <THE0> 

COMMON  ✓IN-'  VAX*  VRV*  XIR*  YIR*  H*  F*  VS*  XIS*  D*  RXN*  BT 
COMMON  ✓OUT'  T*R*RP*  TA*  TDIF*  TL*  HOP * DELPHI  * PRSQ > R * MODE 
COMMON  .'CC'  Cl  * C£*  CS*  CB*  BE*  DS  * CM I N * CMRX  * DM  I N * ETR * 6 5 * GB * RO>  TO  * P I * 
1 Gi 

Xl*<VRX+a-Of'Cl)>>+XIR 
X£=<VS+T)  +XIS 
X3=  <VRY+  <T-  <H'C1> ) ) +YIA 
X4= < <X1-X£) ♦ ♦£) +X3++£ 

X5*SQRT  <X4> 

R1=SQRT a . 0- < <CS-"CB) ♦ ♦£) > 

R£«SQRT  <1 . 0-  < CCMRX'CB)  ♦♦£)  > 

R3=  <R1-R£> 'GS 

R4=SGRT <1 . 0- < <CMIN'CE) ♦♦£)  > 

R5«<R4-R£>/61 

RC=R4'GB 

R7= (R3+R5+R6) ♦£. 0+CB 
IF<X5.GE.R7>  GO  TO  100 
TYPE  1 

t FORMRT <1H  * "SURFACE  DUCT  WITH  SINGLE  BOUNCE  SOLUTION'.') 

THE8=RC0S <CS'CB) 

THE=<180. 0+THE8)  'PI 
2N=. 01 

11  DO  10  1 = 1*9000!.  1 

2£*FLORT  <I> 

2=THE+  < <2£- 1 . 0)  ✓ 1 0 0 . 0)  + .01 
IFC2.GT.90.  0>  GD  TO  400 

TH=ALOG  < < 1 . O'COSD  <2)  ) +SORT  < < <1 . O-'COSD  <2>  > ♦♦£.»  - 1 . 0)  ) 

— ..  i . n r.  \ 

• '.Wl»  * 

TH£*RLOG  <TH1  +SQRT  <.  <THl++£> -1 . 0)  ) 

TH3=  < tTH-TH£)  ♦£.  U>  .'GS 
TH4=CS' CCMIN+COSD <2) ) 

TH5=RL0G <TH4+SOPT < <TH4++£> -1.0)) 

TH6=  < <TH5-TH£)  ♦£.  0)  -"61 
TH7=CS' <CB+COSD <Z> > 

TH8=RLQG <1 H7* SORT  < <TH7++£> -1 . 0)  ) 

TH9= v <TH5-TH8) ♦£. 0)  'GB 
TH1 0=TH3+TH6+TH9 
X 1 = <VAX+ < T- CH  'C 1 ) ) ) +X1 R 
X£=  CS+  .'T+TH1  O'*  ) +XIS 
X3=  CVAY+  a-  <H."C  1 ) ) ) +YI A 
X4=  < <X1 -X£) ♦♦£) +X3++£ 

X5=S0RT (X4» 

El  = <£. 0+CS) 'CDS D (2) 

B£= < (CMRX'CS) ♦♦£) ♦ CCOSD  <Z>  > ♦♦£ 

B3*S0RT <1 . 0-B£) 

E4=  <SIND  <Z) -B3> "GS 

B5=  < <CMIN'CS >♦♦£)♦ CCOSD <Z>  > ♦♦£ 

B6=SQRT  <.  1 . 0-B5) 

E7=<PG-E3)."G1 

B8=  < ’XB'CS ) ♦♦£) ♦ <COSD <2) ) ♦♦£ 

E9=  SORT  < 1 . 0-E8) 

El  0=  <B6-B9) 'GE 
El  1=  <B4+E7+B1 0>  +B1 
B 1 £=  <B  1 1 -X5  > Ell 
E13=RES  <E1£> 

IF  <E13.LE.ZN)  GO  TO  41 
10  CONTINUE 

ZN=ZN+. 01 
Q=.i 
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IF  (ZN.  LE.  O)  GO  TO  II 
400  TYPE  40 

■10  FORMAT < 1 H p'ND  SOLUTION  OF  SURFACE  DUCT  WITH  SINGLE  BOUNCE'/) 

GO  TO  100 
41  THEO-Z 

CALL  ESS (THE 0) 

GO  TO  4£ 

100  TYPE  43 

43  FORMAT <1H  .» 'TRYING  A SURFACE  DUCT  WITH  DOUBLE  BOUNCE  SOLUTION'  ") 

THE9=ACOS  (CS/'CB) 

THE1«<180. 0*THE9)/PI 
ZN1*. 01 

12  DO  99  1*1.9000.1 

Z3=FL0AT(I> 

Z1-THE1+ < (Z3-1 . 0>/100.  0>+.  01 
IF <21. GT. 90. 0)  GO  TO  988 

THB=AL06  < < 1 . O/'COSD  <Z1 ) > +SGRT  < < (1 . O/COSD  <Z 1 ) ) ♦♦£) - 1 . 0> ) 

THB1=CS/  <CMRX*COSD  <Z1)  > 

THB£=ALOb (THB1+SGRT ( (THBl *♦£) -1 . 0) ) 

THB3*  < <THB-THB£)  ♦£.  0) /GS 
THB4*CS/  (CMIN*CDSD  <Z1  > > 

THE5=AL0G (THB4+SGRT < <1 HB4**9) - 1 . 0) ) 

THE6*  < (THE5-THB2)  ♦£.  0)  /G1 
THB7-CS/  <CB*COSD  <Z1)  ) 

THE8=AL0G (THB7+SQRT < (THB7**£) -1 . 0) > 

THE9*  < (THB5-THB8) ♦£. 0) - GB 
THB1 0=THB3+THB6+THB9 
XB  1 * (VAX* (T- (H/C 1 ) ) ) +X I A 
XES-  CVS*  ( I + (£.  U*  I HB 1 0 > > ) +KlS 
XB3* (VAY* (T- (H/C 1 ) ) ) +Y1 A 
XB4* ( (XB 1 -XB£) ♦♦£) +XB 3**£ 

XB5-SGRT (XB4) 

. BE 1 = (£ . 0* CS ) /COSD  <Z 1 ) 

BB£= < (CMAX/CS) ♦♦£) ♦ (COSD  <Z1 > ) ♦♦£ 

BB3*SGRT (1 . 0-BB£) 

BB4* <S I ND (Z 1 ) • BB3)  /GS 

BB5*  ( (CMIN/CS.*  ♦♦£)  ♦ (COSD  (Zl)  ) ♦♦£ 

BB6-SQRT  <1 . 0-BB5) 

BE7* (BB6-BB3) /G1 

BBS- < (CB/CS) ♦ ♦£> ♦ (COSD (Z 1 ) ) ♦♦£ 

BB9*SQRT (1 . 0-BB8) 

BB10=(BB€-EP9)/GB 
BB1 1 = (BB4+BB7+BB 1 0) *BB 1 ♦£ . 0 
BE 1 £* (BE  1 1 -XB5) / BB 1 1 
BB13-ABS  (Bl'l£) 

IF (BB13.LE.ZN1 » GO  10  97 
99  CONTINUE 

ZN1-ZN1*. 01 
Gl-.l 

IF (ZN1.LE.Q1)  GO  TO  1£ 

988  TYPE  98 

93  FORMAT (1H  ."NO  SOLUTION  OF  SURFACE  DUCT  WITH  DOUBLE  BOUNCE""* 

GO  TO  4£ 

97  THE0-Z1 

CAlL  EDS (THE 0> 

48  RETURN 

END 
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SUBROUTINE  ESS '.THE O') 

COMMON  /IN-'  VOX » VfiY  » X I A»  Y I A . H * F » V S ♦ X I S » B • ftXN . BT 
COMMON  'OUT/  T * R»  RP * Tft • TB I F • TL ♦ I'OP . BE • PHI » PPSO« A« NODE 
COMMON  /CC/  Cl > C£»  CS»CB» BE* BS > CM I N • CMftX • BN I N • £ Tft • L ,3 « GB « PO « TO « P I > 
1 G1 

El=  <£.  0*CS>/COSB(THE0> 

p£= < <CMftX'  CS> ♦♦£>  ♦ • CO SB <THE0> >♦♦£ 

B3=SGPT  <1 . 0-B£> 

B4=  (SIND  < THE  O'*  -B3>  /GS 

B5=  < <CMIN.'CS>  ♦♦£>  ♦ 'COSH  (THE CO  > ♦♦£ 

B6*SQRT <1 . 0-B5> 

B7*<B6-B3>/61 

B8=  < <CB/CS>  ♦♦£>  ♦ (COSB  (THE  Co  > ♦ ♦£ 

B9*SGRT (1. 0-B8> 

B10*(B6-B9> /GB 
B11«<B4+B7*B10>*B1 

TH=ftLOG ( (1 . O'COSIl <THE 0>  > + SORT <( <1 . O'COSB (THE0>  >♦♦£•- 1 . 0>  > 

thi=cs/  (CMft»co:n  <theu>  > 

TH£=ftLOG  (TH1  + SORT  ( < TH1  ♦♦£>  - 1 . CO  > 

TH3= < <TH-TH£> ♦£.  Co  /6S 
TH4=CS ■"  <CMIN*CDSB  (THE  CO  ) 

TH5=ftLuG (TH4+S0PT < <TH4**£> -1 . U>  > 

TH6=  < <TH5-TH£>  ♦£,  G>  /G1 
TH7«CS/ <CB*CDSB (THE  0) > 

TH8*ftL06 (TH7+S0RT ( ( . H7**£>  - 1 . 0> > 

TH9--  < (TH5-TH8)  ♦£.  0>  - GB 
TH1 0=TH3+TH6+TH9 
BERTHED 
Tfi=T+THl 0 

Pl-Bl  1*  (SI MB  (THE CO  'CO SB  (THE CO  > 

P£*= < (CMftX/CS) ♦♦£> ♦SIND (THE  CO 
P3* < <CMftX/CS> ♦♦£-' ♦ (CO SB (THE  CO  > ♦♦£ 

P4*SGRT<1. 0-P3) 

P5= (1 . 0- <P£/P4> > 'GS 
P6*=  <CMIN/CS>  ♦♦£ 

P7=S0RT ( 1 . 0- (PE* ( (COSH (THE  0 > ) ♦♦£> > > 

P8=P6/P7 

P9=  (CMAX/CS.)  ♦♦£ 

PI 0-P9/P4 

PI  1=  ( (P8-P1  Co  ♦SIND  (THE CO  ) '61 
Pl£*=  (CB/CSO  ♦♦£ 

P 1 3*S6RT <1.0- (P 1 £♦ < (COSB (THE O' > ♦♦£> > > 

P14*P1£/P13 

P 1 6=  < (PS-P  14>  ♦S I MB  - THE  Co  > /GB 
P 1 5=  < (P5+P  1 1 +P 1 6 ) ♦£ . ©♦C  S > +P 1 
TYPE  3334»B11*P15 
3334  FORMftT (1H  >£E10.4> 

TYPE  3S:c;» 61  * GS  * GB 
S3 35  FORMftT (1H  »3F10.3> 

TYPE  3336 » C 1 > CS  * CMftX » CM I N » CB 
3 336  FORMA T ( 1 H . 5F 1 0 . 3> 

x i = (Vft:  :♦  (T-  (H  • c i > > > ♦;*:  I ft 
X£= (VS* (T+TH1 0>  > +X1S 
X3= <X1-X£> *♦£ 

X4= (Vft Y ♦ (T- <H  C 1 ) ) > *Y I ft 
X5* (XI -X£> /SORT (X3* <X4**£) > 

X6=AC0S (X5) 

PHI  = <130.  0*X6)/PI 
C01*<VS*(T*TH10)>+XIS 
CO£= (VAX* ( T- (H/C 1 > ) > +X I ft 
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9?96 

999? 

9995 


C03*VftX*<C0l-C02> 
C04*-<VftY**T-<H-Cl>>>*Ylft 
C05*C03* ‘ VftY*C04  > 

C06*j0RT  «VfiX**c'+VHY**c;> 

C07*fiF3‘C06> 

C03*<C0l-C02>**2 

C09*C03*  < C 04 ♦ • 2 ’ ♦ • 2.  o*DF  > -D’ **2 
C010*CD?*St>PT<CD9> 

C0il*C05'C0l 0 
SOI*  <COt  s *,C2**2 » 

002*02.  d*C0D*C011)/C2 

003*1.  o>nai-nD2 

004* SORT  <003 J 

OOP  1*1. 0/-D04 

005*  (VS'C2>**2 

006*  < < l OP  r <C08>  > *2.  'C2 

007*1 . 0*005-  'DOG' SORT  'C09»  > 

OOP*  aiOF  1 > ♦SC'PT  aiD7  ‘ 

MODE*xFSS' 

TDIF*T*T  *1 O-TO- <P0'C2>- (H/Cl > 

TO*  <CF*C  D'  D • THE  Co  :•  'CS 
TFl*ftC03 ( TF > 

TF2*  < T B 1 ♦ l 0 . 0 • 'p  I 
IF <£T. EO. 3. 0>  GO  TO  999$ 

IF  i'FT.  EO.  5.  O'  GO  TO 
TYPF  999 6 

FOPMfiTflH  . DIO  MOT  HAVE  ft  VftLUE  POP  feT',-) 
CftLL  FT  3 >.  TF2.  PFl  i 
6Q  TO  999 ? 

CftLL  FTC -TIC.rrL. 

TYPE  99y5.TF2.PFL 
FOPMftT  ‘ 1 H » 2'E  1 0.  3) 

si*  *2.  o*c.?'  x osi«  ■ THE  0> 

S2*fiC0S  < •.C«ftX*COSO  < 1 Ht  0 ■ > sCS > 

S3*  < <<THEO*PI  > '•  ISO.  0>--ic,,./'GS 

S4*ftcos  c <cmih*cdsb  ttheo > > ,-cc> 

C5*  (?4~  (3 1 

S£*ftCOS  <,  >XF*COSD ' THEO>  ,»  ✓CS> 

S7*  <S4-SC.>/GF 
S8*  (S3*  i':.*37»  *S'i 
QZ1*-  < <«♦  tPBL.)  r 1 0.  0 
QZ2*  <1 0.  ('♦♦Or  1 > 

P01*s.  0*  (ftxr  *♦2)  ♦ ' IMP f THE 0> ♦COSO <THEO>  *0^9 
PRSO*FO  1 ^ < fc  1 1 ♦ftj:  i ..  p i f.  > > 

TL*  <1  0.  ('♦HLOGl  t'i . F P:'C> ’ 

RETURN 

END 


SUBROUTINE  FO  3 < THE  0 > 

COMMON  'IN  ' Vft  X . Vft  Y , J : I ft . V I ft , H , F . VS  , X I S , D . ft/ri . y ^ 

CnMHnr!  ' - - T/ r ] ' ? ’ ™ ' T 11 1 F ' TL  » rDf'  * OE  • PHI . EPIC',  ft.  MODE 

COMMON  HC/  Cl . C2.CS.CF.  OF.  u S . CniM.CMftX,  DMIN.  tTft.GS  «GF«  PQ.  TO. PI . 

FI*  <2. 0*CS> 'COSO < THE O’ 

F2* < < C MftX'C S * ♦ ♦Zv ♦ r r OSD i THE  0 > • ♦♦2 
F3*SC'PT  • 1 . 0-F2 ' 

®4*  I f ID  • THE 0 ) -F ’ '<91 

F5*  ‘ • C M I N ' • **2  • ♦ • COCO « THE  0 » • 

F6*S0PT <1. 0-F5 • 


F7*  <F*-F3 • '61 

F8*  < <CFX  1 > **2  • ♦ 'COSI"  THtO*  • ♦♦Z 
F9*3QPT  < 1 . U-fciEu 
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El  0*<E6-B9>/GB 
B11*<B4+B7«-B10>*B1 

TH*HL06 < <1 . 0/C030<TH£Ci>  >«-SGRT<<U.  O/COSDcTHEO;  >♦♦£>-!.  0>  > 
TH1  *03/  < CNfiX^CO  SD  THE  0 > > 

TH2*flL06  (.TH1+  SORT  < <TH1 ♦♦£>-! . 0>  > 

TH3*  < aH-THE)  0*  /GS 
TH4*CS/  <XM  I 0 a HE  0 > > 

TH5*ftLOG<TH4+SGRT < <.TH4»*2>-1 . 0> ) 

TH6=<<THf>-TH2>  ♦£.  OV'Gl 
TH7*CS/<CB*CQSD<THE0>.> 
rH3*HL06  fTHT+SORT  <.  < rH7^2>-l.  0>  > 

TH9*  <.  (TH5-  THE  ’ ♦£.  0>  'GB 
TH1 0= TH3+ TH6+TH9 
D£*THE0 

Tfl=T+  <2. (•♦  TH 1 0 > 

P 1 *B1 1 ♦ <S IHD  <THE  0 » /CO  3D  < THE  0 > > 

Pc*  < <CNflX/CS>  ♦♦£ > ♦SINE  '.THE 0> 

P3*  < <Cf!rlX/CS>  ♦♦£  > ♦ CCOSD  < THE 0>  > ♦ ♦£ 

P4*SQRT  <1 . 0-P3> 

P5=  <1 . 0-  CP2/P4 ) .)  /GS 
P6*  <CMIM'CS>  ♦•£ 

P7=30R’’  <1. 0-<PG»  <<COSD<THEO >>++£>>  > 

P8=PG'P7 

P9= <CMHX/CS> ♦♦£ 

pi  f.* P9/P4 

PI  1 * < CPS-P 1 0.)  ♦ i I ML  •-  THE  0 > .'61 
PI  2* (CB/CS.) ♦ ♦£ 

P13*3QRT  <1.0- <P1 2* < (CDSD <THE0> >♦♦£> >> 

P14*P12/P13 

p«  c - / ('po^c i rj£|  itu: iv>  > .i;v 

P15*  < <P3*P  1 1+P16 • ♦£ . 0»CS/'  +P  1 
X 1 * < < T - <H  'C  1 ) > > * X I H 
X2* <VS*  <T+ (2. 0*TH10>> >+XIS 
X3* <X1-X2> ♦♦£ 

X4*  <VrlY*  <T-  <H'C1 > > >*YIH 
X5*  <X1 -X2>  /SORT  <X3* « X4**2>  > 

X6*H'S0S  Z> 

PHI*  <.130.  0*X6> /PI 
COl*  <VS*  (T*  <2.  0* TH 1 0>  > > ♦ XIS 
C02*  /v‘P.->  •! T—  <H'C  1 ) > > ♦XlH 
C03*VPX*  <C01-CQ2> 

CQ4*- CVR {* < T- * H 'C 1 > > > +Y I* 

C05*C03+  (v'hY*C  CI4 » 

CU6*S0RT  Hr.**c  + Vrl  r ♦♦  c':* 

C07*HBS  <C06> 

C08*<C01-C02‘ ♦♦£ 

C09*C03+  <'iID4##2>  ♦ • (4.  (•♦DB> - I'1  ♦♦£ 

CO 1 0*C  0 '♦ i OP T < C 0 ? > 

coil  *»:or.'C  ai  o 

D01  * • / •'££•♦♦£■  > 

D02*< <2. 0*CD7* ♦COl 1 > /C2 
D03*l.  0*D01-D02 
D04*SC'PT‘D03> 

DOP1*1.0'D04 

D03*cVS/C2.***2 

DOC*  < (SORT  'CO'3'  >♦£.  C'^VI  > 'C2 

DO 7*1 . (•♦DO*-  • DOE-/ 1 OPT  -CO?'  » 

DOP*  «D0P1 •♦SORT ■ DO 7* 

TB*<CE*CO:Df  THEC»>  > 'CS 
Tn*flco:<Ti.> 

TB2*tTBl*130.  0>  PI 
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IF  <f:T . EC.  5.  0>  GO  TO  9998  ' 

IF  (FT.FO. 5.  0»  GO  TO 
TYF'E  9998 

9996  FOPMftT * I H . DIO  MOT  HAVE  ft  VftLUE  FOP  FT'  > 

9998  CftLL  I:T3‘TE:i,PfcL> 

60  TO  9997 

9999  CftLL  FT5<TF£»PFL> 

9997  TYPE  * 9995 . T F£ . PPL 

9995  FOPMftT  <1H  *£E10. 3> 

SI*  <£. 0*CS> - CO! D ' THEO> 

s£*ftco:  < uiMftx^co:  o • theo  • » ✓co 
S3*  ( < <THE 0*P  1 > •'  iso.  0*-i£>  ^Gi 
S4=ftcos < •cmim*i:o:d • thecd  • /csj 
S5=  <S4-:£>  •'Gl 

S6*ftC  0 S t'C D»COSD  < THE 0 * > ✓CS > 

S7*  <S4-S6> /GF 
S8*  <S3+S5*S7* ♦£ 1 

OZZ !*-<•£. (•♦ft*i8>  + <£. 0*PPL> >s 10. 0 
0"Z£=  <1 0. 0**CCZ1 * 

poi  *8.  ■ «xm**£)  ♦:  nui  <theo>  ♦co:d  .the  o • ♦ozzs 

PRSQ*F'C'l/  < < <B1 1 ♦ftf  i.  - PIS  . > .*4. 0» 

TL*  <1  0.  C^ftLOGl  '•  <PF  ZC"  .• 
riODE*'FD'  ' 

TDIF*T4'£.  ('♦VHiO>-TO-i.PO  Cc'^-'H  Cl) 

RETUF'M 

END 


Bolt  Beranek  and  Newman  Inc 


TM  W307 


TYPE  BT3.F4 

SUBROUTINE  BT3 <T££» RBL) 

COMMON  /IN'  VOX * VAY » X i A . Y I ft  * H * FS  VS  t X I S» D » RXN • BT 
COMMON  /OUT/  T»R»RP» TA*TDIF»TL»  DOP«  DE»  PHI »FRSQ»R*  -'10DE 
COMMON  /CC/  Cl  t C£»  CS»  CE  ? D£»  DS » CM  IN*  CMRX  « EM  I N»  ETA»  GS»  GB  ? RD « TD  t P 1 1 
1 61 

IF  <F.LE. 300.  0)  GO  TO  1 

IF<F.6T.300.  O.RND.F.LT.  750.  0)  60  TO  £ 

IF<F.6E.?50. O.ANB.F.LE.  1500. 0>  6D  TO  3 
IF <F . GT . 1500. 0. AND.  F. LE • £700. 0)  GO  TO  4 
IF<F.GT.£700. O.ANB.F.LE.  5000. 0)  60  TO  5 
IF  (Ft6T,500il.  0)  GO  TO  6 
TYPE  556 

556  FDRMRT <1H  » 'OUTSIDE  FREQ.  RANGE-  OF  PROGRAM'-/) 

1 IF<TB£.LE. 11.0)  60  TO  8 

GO  TO  £00 
3 RBL*0. 0 

GO  TO  100 

£00  IF  <TB£. GT. 11.0. AND.  T££. LT. 50. 0)  GO  TD  9 
GO  TO  300 

9 R£L= <<10. 0*T££) /39 . 0> -£ . 8£ 

GO  TO  100 

300  IF  <TE£. GE. 50. 0)  GO  TO  10 
GO  TO  500 

10  RBL=10.  0 
GO  TO  100 

£ IF<TB£.LE.  13.  0)  GO  TO  1 1 

GO  TO  400 
i 1 p»i  —3  n 

GO  TO  100 

400  IF  <TE£.GT.  13.0.  RND.  TE2.LT.  £0.0)  GO  TO  t£ 

GO  TO  600 

1£  RBL*  < <£. 3*TE£) /?. 0) -1 . 27 

GO  TO  100 

600  IF<TE£.GE. £0.0. RND.  TE£.lE. 35.0)  GO  TO  13 

GO  TO  700 

13  RBL*  < <3. 4*TE£> /1 5. 0) -.77 
GO  TO  100 

700  IF  <TE£. GT. 35. 0.  AND.  TE£.  LT. 5£. 0)  GO  TO  14 
GO  TO  701 

14  RBL*  < <1 . 6*T££) /1 0.  0>  +3. 1 
GO  TO  100 

701  IF  <TE£.  GE.  5£.  0>  GO  TO  15 
GO  TO  500 

15  RBL=1 1 . 0 
GO  TO  100 

3 IF<TE£.LE. 15. 0)  GO  TO  16 

GO  TO  70£ 

16  RBL* 3. 0 
GO  TO  100 

70£  IF<TB£.GT.  15.0.  AND.  TB£.LT.  45.  0>  GO  TO  17 
GO  TO  703 

17  RBL* < <8. 0»TE£)/30. 0>-l. 0 
GO  TO  100 

703  IF  <.T££.  GE.  45.  0)  GO  TO  18 

GO  TO  500 

18  PEL* 11.0 
GO  TO  100 

4 IF <T££. LE. 13. 0>  GO  TO  19 

GO  TO  704 
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19 

REL* 3. 0 

GO  TO  100 

7 04 

IF  <TE2. GT. 13.  O.AND.  1 > T.42. 0) 

GO  TO  705 

GO  TO  £0 

£0 

REL*  < <.3.  0+TE£>  /£9.  0>  59 

GO  TO  100 

705 

IF <TB2.GE.42. 0)  GO  TO  21 
GO  TO  500 

21 

RBL=1 1. 0 
GO  TO  100 

5 

IF  <TE£. LE. 2. 5>  GO  TO  22 
GO  TO  706 

22 

REL=7. 0 
TYPF  1111 

1111 

FORMRT  <1H  * 'REL  HRS  ERROR  GREATER  THRU  £ 
GO  TO  100 

706 

IF  <TE2. GT . 2. 5. AND.  TB2.  LE. 12.5) 
GO  TO  707 

GO  TO  23 

23 

REL* < <7. 0+TE2> /\ 0. 0> +5. £5 

TYPE  1111 
60  TO  100 

707 

IF  <TB2.  GT.  12.5.  AND.  TB2.  LT . 32 . 5'1 
GO  TO  703 

GO  TO  24 

24 

REL*  < < 1 . 5 ♦ T E2> / £ 0 . 0 > + 1 3 . 06 
TYPE  1111 
GO  TO  100 

IF CTBS. GF. ?2. 5'  TO  25 

GO  T0500  " 

r>- 
^ • 

REL* 15. 5 
TYPE  1111 
GO  TO  100 

6 

IF  <TB2.  C.E.  7. 5')  GO  TO  26 
GO  TO  709 

26 

RBL=S. 0 
GO  TO  100 

709 

IF <TE£. GT. 7. 5. AND. TB2. LT. 15.  0> 
GO  TO  710 

GO  TO  £7 

27 

REL*  (.  <5.  0+TE£>  ^7. 5>  >3.  0 

GO  TO  100 

710 

IF  <TB2. GE. 15. 0.  AND.  TB2.  LT. £5. 0> 
GO  TO  711 

GO  TO  23 

2S 

REL*  < k3.  0+TE2>  -'10.  0)  +3.5 
GO  TO  100 

711 

IF  <TE2. GE. 25. 0>  GO  TO  29 
GO  TO  500 

29 

PEL* 16. 0 
GO  TO  100 

500 

TYPE  555 

555  FORMAT  (1H  •.'CANT  FIND  REFLECT  I DM  LOSS'-  '> 
100  RETURN 
END 
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TYPE  BT5.F4 

SUBROUTINE  BT5<TB£»R£L> 

COMMON  /IN/  VAX  * VAY » X I A • Y I (=1 , H*  F > VS  > X I S > B > flXN * BT 
COMMON  /OUT/  T*  R*  RP*  TA»  TDIF*  TL*  DDF*  DE > PHI , PRSQ,  ft,  MODE 
COMMON  /CO/  C1*C£*C$»CB* DE* DS» CMIN* CMAX. DMIN*ETA>GS*GB»RO. TO* PI t 
1 G1 

IF  <F. LE. 300. 0)  GO  TO  1 

IF<F. 6T.300. O.AND.F.LT. 750. 0>  60  TD  £ 

IF  <F.  GE.  75>0.  0. AMD.  F.  LE.  1500.  0.)  GO  TO  3 
IF  <F.  GT.  1500.  0.  FIND.  F.  LE.  £700.  O')  GO  TO  4 
IF<F. GT. £700.0. AND. F.LE. 5000.  0>  GO  TO  3 
IF <F. GT. 5000.  0>  GO  TO  6 
TYPE  556 

556  FORMATOH  * 'OUTSIDE  RANGE  OF  FREDS.  '/> 

GO  TO  100 

I IF <T££. LE. 5. 0)  GO  TO  8 
GO  TO  £00 

8 RBL*4. 0 
GO  TO  100 

£00  IF  (TB£.  GT . 5.  0.  AND.  TB£.  LT.  £1 . 0>  GD  TO  9 
GO  TO  300 

9 R£L«<<10.  0»TB£.)/16.  0.>  + .88 
GO  TO  100 

300  IF  <TE£. GE.  £1 . 0)  60  TO  1 0 
GO  TO  500 

10  R£L*14. 0 
GO  TO  100 

£ IF <TB£. lE. 3. 0)  GO  ^0  11 

GO  TO  400 

II  REL*=8.  0 
GO  TO  100 

400  IF  <TB£. GT. 3. 0. AND. TB£. LT.  ££.  0>  60  TO  1£ 

. GO  TO  600 

1£  RBL*<:<,8.  0»TB£>/19.  0>+6.74 

GO  TO  100 

600  IF (T££. GE. ££• C>  GO  TO  13 

GO  TO  500 

13  RBL*16. 0 

GO  TO  100 

3 IF <TB£. LE. £. 5.)  GO  TO  16 
GO  TO  70£ 

16  PEL *9. 0 
GO  TO  100 

70£  IF'TB£. GT.£.  5.  AND.TEc.LT.  17.  5>  GO  TO  17 
GO  TO  703 

17  RBL*<<9. 0«TB£)/15. 0 +7.5 
GO  TO  100 

703  IF<;TB£.GE.  17.5>  GO  TO  18 
GO  TO  500 

13  P£L=13. 0 

GO  TO  100 

4 IF  <T6£.  LE.  £.  5>  GO  TCI  19 
GO  TO  704 

19  RBL*9. 0 

GO  TO  100 

704  IFaB£.GT.£.5.AND.TB£.LT.£0.  0>  GO  TO  £0 

GO  TO  705 

£0  RBL*  < <9.  0*T££>/17.5>+7.71 

GO  TO  100 

705  IF <TE£. GE. £0. 0)  GO  TO  £1 
GO  TO  500 
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£1 

5 

££ 

706 
23 

707 
£4 
703 
£5 

e 

£6 

709 
27 

710 

•*r» 

711 

£9 

500 

555 

100 


RBL-18. 0 
GO  TO  100 

IF <TB£. LE. £. 5)  GO  TO  £2 
GO  TO  706 
RBL«1 1. 0 
GO  TO  100 

IF<TB2.GT.£.5.flND.TB£.LT. 1£.5>  GO  TO  £3 
GO  TO  707 

RBL  = < <7. 1 *TB£,>  /"l 0.  0>+9.£3 
GO  TO  *100 

IF  < TB£.  GE.  12.5.  AND.  TR2.LT.  £2.5)  GO  TO  £4 
GO  TO  703 

RBL*  < <1. 9*TB£> /1 0. 0>  + 15.  73 
GO  TO  100 

IF  <TE£. GE. ££.  5>  GO  TO  £5 
GO  TO  500 
REL*£0.  0 
GO  TO  100 

IF <TB£. LE. £. 5)  GO  TO  £6 
GO  TO  709 
RBL*10.  0 
GO  TO  100 

IF <TB£. GT.£. 5. RND. TB2.LT.  10.0)  GO  TO  £7 
GO  TO  710 

RBL*  < <6. 0*TB£) /? . 5) +8. 0 
GO  TO  100 

IF<TB£.G£. 10. O.fiMD. TB2.LT.20. 0)  GO  TO  £8 
GO  TO  711 

r»f>|  r / /-5  rt*TD*is  * * ft  ft  s .»  4 <-»  ft 

'-••fc-  ' 'W*  . • I A/k/  • « Wl  V/  • V 

GO  TO  100 

IF<TE£.GE.£0.  0>  GO  TO  £9 

GO  TO  500 

RBL*19. g 

GO  TO  100 

TYPE  555 

FORMRTOH  * 'CftNT  FI  HD  REFLECTION  LOSS'/' 

RETURN 

END 


Reproduced  from 

best  availeble  copy^^y 


T 'PE 


£00 


ERFS.F4 

SUBROUTINE  ERFS  ■'hRG*  ERF ) 
S8RG*»t|RG  '8B S oiiRG) 

X«RES<ftPG> 

ER«1 . 0 

IF  <X. GT , 4 . £ > GO  TO  £00 

E*1 .0+.  07052 30784e>:+ , 04££8£ 

1 . 000152U143* <X**4> + . 00027 

EP*1 . 0-1 . O''  'E**16> 

ERF»EP*SPRG 

RETURN 

END 


01£3*X*X+. 009£705£7£e • X** 
656 7£e <X**5»  + . 0000430623* 


;*  • 4 
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